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PSErOO-RANDOM AND RANDOM TESTING 


Robert S. Norin 
Tinie/Data Corporation 
Ptilo Alto, California 


This paper describes control techniques* hardw'are implementation* and opera- 
tional aspects of current digital control systems used for pseudo-random and 
random vibration testing. 


INTRODUCTION 

The advantages of digital random vibration control 
systems* when compared with analog systems, have 
been well established since the first digital systems 
were introduced in the late 1960’s and early 1970’s 
(1* 2* 3* 4). Since this time* digital random control 
systems have evolved in their sophistication and per- 
formance. This is a trend primarily due to the in- 
creasing awareness of the potentialities of digital 
systems on the part of environmental testing person- 
nel and manufacturers of digital systems, as well as 
the rapid evolution of higher performance processing 
hardware. The impact of increasing integration and 
decreasing size and price of key electronic modules 
has enabled manufacturers to offer basic digital vi- 
bration systems which are more than competitive with 
analog control units on the basis of cost alone. 

The purpose of this paper is to describe the control 
techniques, hardware implementation, and operational 
aspects of a typical state-of-the-art digital random 
vibration control system. The first portion of the 
paper discusses a number of techniques for dealing 
with important technical asp^icts of the random vibra- 
tion control problem. These include the generation 
of pseudo-random and true random noise, the control 
spectrum estimation problem, the accuracy /speed 
tradeoff, and control conection strategies. The sec- 
ond half of the paper deals with system bar iware, the 
operator-system interface, safety features, and op- 
erational capabilities of sophisticated digital random 
vibration control systems. 


BASIC CONSIDERATIONS 

The goal of any random vibration control system is to 
excite a shaker system with a continuous random 


drive signal whose power spectrum is adjusted in such 
a way that the power spectral density (PSD) of a re- 
sponse accelerometer signal (or averaged spectrum 
of a number of signals) conforms exactly to a specified 
reference spectrum. 

The basic function requirements for digital random 
control operation are; 

1. Random noise generation. 

2. Response (control) spectrum estimation. 

3. Spactrum modification. 

Figure 1 shows a functional block diagram for a typ- 
ical digital random \ibration contra’ \stem. 

The timing flow for a typical implementation of these 
functional tasks in a digital system Is illustrated in 
Fig. 2. 

We now proceed to examine in detail the three basic 
functional requirements as they apply to a digital 
control .* ystem. 

GENERATION OF RANDOM AND PSEUDO RANDOM 
SIGNALS 

[f we are given a particular drive spectrum defined by 
a certain finite bandwidth and number of spectral lines, 
then by adding a random phase angle to each of the 
spectral lines to produce a complex Fourier spectrum, 
and taking the inverse Fourier transform, a frame of 
a time-domain signal with the required spectral char- 
acteristics is produced. Repeating this operation con- 
tinuously and concatenating the resulting frames, as 
shown in Fig. 3, produces a continuous random sig- 
nal (1, 3), In order to prevent discontinuities and to 
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RANDOM NOISE GENERATOR 


FEEDBACK CONTROL 


I 


POWER SPECTRAL DENSITY ANALYZER 


rig. 1. Functional Block Diagram for Digital Ramlom Vibration Control System 


maintain the proper correlation between frames, it is 
necessary to use u smoothing or windowing operation 
at the boundary between adjacent frames. By ac- 
counting for these and other subtleties, the random 
noise produced by this system is indistinguishable 
from the noise produced by a natural Gaussian Ran- 
dom noise source, but has the advantage that its 
spectrum can be rapidly changed to any arbitrarv' 
shape with a frequency resolution of 0.1^ or less. 

The normalized amplitude distributio n of a random 
noise signal generated by this digital process is shown 
in Fig. 4. This distribution, as a function of standard 
deviation units, is displayed both logarithmically and 
linearly. The Gaussian distribution curves are super- 
imposed on the experimental results. 
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Fig. 2. Tvpical Random Control I.oo|> Timing 


Although the original spectrum from which the digital 
random signal was derived had onl}* a finite number of 
spectral lines, the spectrum of the random signal, 
when analyzed with finer resolution and averaged over 
a large number of frames, is seen in Fig. 5 to ap- 
proach the continuous spectrum expected from a 
Gaussian random signal. 

A large number of frames of data are required to ob- 
tain a meaningful spectrum because of the statistical 
uncertainty or variance associated with any measure- 
ment of a random signal. The variance problem can 
be circumvented by generating a ^’pseudo-random'' or 
’’zero variance” signal. Basicially this is accom- 
plished by initially randomizing the phases of the 
spectral components of the drive signal and maintain- 
ing them constant throughout the test, while adjusting 
the magnitudes every L frames. Thus, the same 
frame is repeated L times. L is usually chosen to be 
large enough that a steady state condition is reached 
prior to Unalysis. If the anal>sis frame period T is 
the same as the output frame period, then a single 
analysis frame is sufficient to determine the spec- 
trum. 

Because of periodicity of the w'aveforra, energy 
will be produced only at those frequencies defined in 
the drive spectrum, i.e., the spectrum is "leakage 
free” and contains "holes” as illustrated in Fig. 6. 

Similarly, the amplitude histogram of a pseudo- 
random signal (.shown in Fig. 7), while Gaui»sian-like, 
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contains structure not seen in the random case. 

While pseudo-random noise can be used in principle 
to meet many random test specifications, it lias not 
found wide acceptance in vibratior laboratories be- 
cause of its peculiar ''sound discrete spectrum, 
discrete amplitude structure, and historical lack of 
precedence in analog random equipment. 
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Fig. 3. ('lencration of Continuous Handom Noise 
Signal 
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Fig. 5. Spectrum of Digital Random ^’oise. 
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Fig, 4. Theoretical vs Experimental Normalized 
Amplitude Distribution Curves for Digital 
Random Noise Signal 


Fig. G. Spectrum of Digital Pseudo-Random Noise 
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Fig, 7. Typical Amplitude Distribution Curves 
for I^eudo-Random Noise Signal 
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PSD ESTIMATION 

If random noise is to be used fcr testing* then because 
of the statistical uncertain^ of the spectral measure* 
ments one is faced with a tradeoff between control ac* 
curacy and system respocse time* At the heart of the 
control problem is the task of estimating the spectrum 
(PSD^ of the response signal we are trying to control* 

A convenient measure of the statistical reliability of 
a spectral estimate is its ''equivalent number of de- 
grees of freedom (k)'\ defined as: 


It can be shown that for linear spectral averaging* 
k = 2 for each frame of Gaussian distributed random 
data (5). By way of illustration, after an estimate 
with 16 degrees of freedom (d.o.f. ) there is 90% con- 
fidence that a particular spectral \^lue is within an 
interval of about 5 dB around its correct value; for 
64 d.o.f. the interval is about 2.5 dB; for 256 d.o.f. 
the inter\-al is about 1.25 dB. 

The response PSD can be estimated using either of 
two methods; direct and indirect. Both methods are 
described in the following paragraphs. 


y(t)* in the tlme*domain for M consecutive frames of- 
uniform length T* The PSD of each of the two summed 
data frames is then computed* The above process is 
repeated N times with the autospectra being linearly 
averaged. Finally the estimate |y|^ of the response 
spectrum IY|2 is given by; 


ivr == 


2 M ,„,2 


I 


where jYj^j | ^ is the PSD resulting from the 

sum of M frames of the response signal y(t), 

I^Mj I ^ the PSD resulting from the 
sum of M frames of the drive signal x(t)* 

and |x I ^ is simply the drive spectrum used when 
^Mj ^Mj obtained. 


It can be shown (7) that the variance in the PSD esti- 
mate using equation (2) is inversely proportional to 
MN, while for a fixed M and N, the variance tends to 
be greatest around resonances. 


The Direct PSD Estimate 

The traditional method of estimating the response 
PSD, illustrated in Fig, 2, is to acquire a frame of 
data and then directly compute the PSD from the 
sampled data by Fast Fourier Transform (FFT) tech- 
niques (2, 3). As PSD' s are accumulated, they are 
averaged. If the averaging is linear, a drive spec- 
trum correction is made after N PSD's are averaged, 
and then a new average is accumulated. If exponen- 
tial averaging is used, a running average with a time 
constant of N frames is maintained, and a drive cor- 
rection is made after each frame. 


The timing flow for a typical implementation of this 
indirect PSD estimation method with M = 4 and N = 2 
is shown in Fig, 8. 


ACQUtRE INPUT P 


SVK INPUT 
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ESTIMATE PSD 

COMPUTE ERROR 
AND MODIFY 
DRIVE SPECTRUM 

GENERATE 
DRIVE SIGNAL 


The variance in the direct PSD estimate is inversely 
proportional to N for both linear and exponential aver- 
aging (6). 


COMPUTE OUTPUT 

SIGNAL ><UUMUUMUMMMUUUUUUUUUUU4MHVIMMI>* 

OUTPUT DRIVE ^ , ■ ■ , , 

SIGNAL ^ ^ ^ ^ ^ f ^ ’ l-H-H-l- 


The Indirect reD Estimate 

A new algorithm has recently been introduced (7) 
which significantly improves the speed and accuracy 
with which an estimate of the response PSD is made. 
This algorithm is based on a combination of time- 
domain and frequency-domain techniques. 

Basically, *he procedure used is to first separately 
sum the drive signal, x(t), and the response signal, 


Fig. 8. Indirect Method Random Control Timing, 
N^2, M-4. 

There is a twofold advantage to the Indirect method. 
First, since a time-domain summation calculation is 
faster than a PSD calculation, and only two PSD cal- 
culations are needed for each M data frames, more 
input (response) data can be processed per unit of 
time using the indirect method than with the direct 
method. Table 1 summarizes the percentage of time 
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Table 1* Percentage of Time Devoted to Input Data Acquisition for Different 
Parameters and Hardware Configuration 


CONTROL 

bandwidth 

HZ 

PROCESSING 

HARDWARE 

(Note) 

INDIRECl 


DIRECt 1 

N=2 

M=2 

N=2 

M-4 

N=2 

M=8 

N=2 

N=4 

N=8 

1024 

1 

22 

37 

55 

18 

19 

20 

1024 

2 

30 

42 

66 

27 

35 

42 

1024 

3 

33 

52 

69 

28 

36 

43 

2048 

1 

13 

22 

35 

10 

11 

11 

2048 

2 

20 

36 

52 

18 

23 

29 

2048 

3 

25 

41 

59 

21 

27 

30 

3072 

1 

7 

13 

16 

5 

6 

6 

3072 

2 

12 

24 

37 

10 

15 

16 

3072 

3 

20 

34 

51 

16 

25 

29 


Note: 1 ^ CPU (PDPll/35) only 

2 = CPU Hardware FFT Processor 

3 5= CPU + Hardware FFT Processor Output Processor 
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that is spent acquiring input data for both methods as 
a function of M and N, control bandwidth, and typical 
processing hardware. Input acquisition is discontin- 
uous, i. e. , occurs less than 100% of the time, due to 
the additional processing functions required during 
the control loop operations including random dn .*e sig- 
nal generation, drive spectrum correction, PSD esti- 
mation, displays, safety checks, and miscellaneous 
bookkeeping. It is observed from Table 1 that by ad- 
ding hardware to perform tasks such as Fourier com- 
putations and output signal generation in parallel with 
other functions, input data can be processed more than 
50% of the time for control bandwidths up to 3000 Hz 
using the indirect method with N = 2 and M = 8. 

The second advantage is that for a given number of 
input data frames the statistical variance of the indi- 
rect method estimate is significantly lower than the 
corresponding direct method estimate. It has been 
observed empirically that the variance of a indirect 
PSD estimate, although somewhat dependent on the 
characteristics of the load being driven, is typically 
less than the variance obtained with a direct PSD 
estimate using four times as many data frames. In 
terms of statistical degrees of freedom using linear 
PSD averaging this means that for the indirect meth- 
od, k 8MN, while for the direct method, k = 2N. 

Figures 9 through 12 illustrate the reduction in var- 
iance using the direct method. Figures 9 and 10 show 
direct and indirect PSD estimates, respectively, of a 
nominally flat 2000 Hz spectrum, each using 8 input 
data frames. A tolerance interval of 4 dB (± 2 dB) is 
superimposed on the estimates. Figures 11 a.id 12 
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illustrate corresponding results with 32 data frames. 
Comparing Figures 10 and 11 it is seen that the in*- 
direct estimate with 8 frames is at least as good as 
the 32 frame direct estimate. Similarly, the variance 
in the 32 frame indirect estimate (Fig. 12), ±1 dB, is 
comparable to that expected from a 256 degree of , 
freedom estimate, requiring 128 frames using the 
direct method. 



rcMiTOR sPErmm — m control loops 

DIRECT METHOD/ N«B (bH DCF) 


Fig. 11. Direct Estimate of Flat Spectrum, 
32 Frames 



Fig. 12. Indirect Estimate of Flat Spectrum, 
32 Frames 


THE COMPLETE CONTROL LOOP 

Overall random vibration control is achieved by com- 
bining a rapid, accurate PSD estimation technique 
with a bold, yet stable, drive spectrum modification 
strategy. A typical strategy for drive spectrum mod- 
ification is the following (3): For each spectral com- 
ponent, determine the error as the ratio of reference 
to control PSD estimate. Convert the ratio to a loga- 
rithm and correct the drive by one-third of the loga- 
rithmic value, e.g, , a 3 dB error produces a 1 dB 
correction, a 6 dB error produces a 2 dB correction, 
etc. 


Figure 13 illustrates a typical control point PSD esti- 
mate using the indirect method while controlling ' 
through a peak-notch transfer function with a 36 dB 
dynamic range. The control strategy implemented is 
based on the method outlined above. Tolerance bands 
of ± 2 dB are superimposed on the plot. The corre- 
sponding drive spectrum is shown In Fig. 14. Note 
that the variance in the estimate tends to be 
greatest at the 500 Hz resonance and 1000 Hz anti- 
re'sonance. Figures 15 and 16 show corresponding 
control and drive spectra using the same control strat- 
egy and the direct PSD estimation method with N = 4. 
Notice Chat the variance using the indirect method is 
significantly less than using the direct estimate, even 
though the total acquisition and processing time was 
only 20 to 40% longer, depending on system hardware, 
for the indirect method (both estimates required 4 PSD 
calculations). The cumulative effect of a good PSD 
estimate and control strategy is Illustrated by compar- 
ing the drive spectra. Figures 14 and 16. A smoother 
drive spectrum could be obtained in the direct case by 
adopting a more conservative correction strategy, but 
at the expense of overall correction time. 



TYPICAL CONTROL POINT PSD cBTIMMTE 
INDIRECT METHOD/ N»S' 


Fig. 13. Typical Control Point l^D Estimate, 
Indirect Method, N 2, M 



TYPICAL 3b OB DRIVE SPECTRUM 
INDIRECT METHOD/ N>2/ M«B 

Fig. 14. Typical Drive Spectrum, Indirect Method, 
N-2, M=8 
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A representative digital random vibration control 
system Is shown in Fig* 17* A block diagram of tbjt 
basic system components and options Is shown in 
Fig. while Fig* 19 shows □ close-up of the con- 
trol runnel as used for random vibration control o|Jcr- 
atfon. A unkiue ndvantige of many dl-^tnl systems in 
that the same system hanlware can lie used for other 
control ami analysis applications simply by loading 
different o|x?ratlng programs. The system shown 
here, for example, uses the same control panel and 
system hnixhvaro for random acoustic contra ' , swept 
narrowband random on wideband random \>1bration 
control, and swept sine on random ^dbrutlon control 
orieraUon l\\ replacing the overlay on the right 
half of the control jianel and loading ajjfiroijrlatc soft- 
ware the svslem becomes a r inusokhl (91 or trnn- 


slenl (Id) control svsten*. 


Kig. 15, 1 ypical Control mint PSD Fstimatt 
IKrect Method, X 4 


DIGITAL SYST KM TLAinAVARK 

Baste hardware requirements for a digital random vllj- 
ratton control system include a central processing unit 
(typically a IG-bit mlnf computer) with associated mem- 
ory (Ifi to 32 thouFand words), programmable analog- 
to-dlgital and dlgltaMo-analog converter systems 
flO or 12-bits, 1 to 8 Input channels, with a maxlin^jm 
sampling rate at least 20 kHz), programmable filters 
for Inixit and output (low -pass, minimum 48 dB/octnve 
roUofO» a display device (typically 4^' x 5'" CRT), a 
program loading device (typically a hlgh-spc^d ixiper 
tarxj reader), and a system control panel and/or an 
alphanumeric keyboard, A variety of peripheral de- 
vices such as rjcirallel processing units, mass storage 
devices (magnetic tape units, disks, or cassettes), 
large screen displays, and hard-copv' devices IX- Y 
plotters, electrostatic copiers, etc*) enhance lx)th the 
performance and convenience of a basic digital control 
svstem. 


Fig. n. 


T\pical nigilal Handom Vibration (\mtnd 
Svstem 


SVSTEM OPERATION 

0()crator interaction with the system described here 
is through the control panel (Fig. 19) and aljihanu- 
merlc kevboaiY] terminal. Test ijara meters arc Ini- 
D 1 1 ally entered by the ‘nerator from the kevlKMird, A 
typical dialog lietween P operator and the system Is 


AH- 







illMStnited in Fig. 20, U|jerator responsijs are under- 
lined. The total to be eon' rolled, filter 

5 pJ»eing, and the GUMS value of the reference spee- 
tnim are computer- calculated responses. Once test 
piira meters hare Isoen entered, they can be edlteci 
through a correelions routine, listed for a printed 
record, and saved on disk or punched jmper laiie to 
allow for fast, accoiratc setup of future tests. 

a 

!f more thiin one control channel Is specified, then _ , 

the average of the estimate<l rSD's of each selected '*1 
channel is used as the spectrum to he controlled. 
Additionally, the sjiectram of any channel (''auxiliary'' 
channel^ mn be indejxjndenily displayed during and 
after lest operations. The auxiliary elrinnel definition 
ean he changed during test ofjcrations* Control panel 
switches also |iermit the o|ierator to nrbitrarilv mod- 
f fT the refe rente s pec t nmi , e ve n w h i li- the les t i s 


With test parameters entered, the operator ain ini- 
tiate the test at any time from the control panel by c 
selecting the TIi:ST function and pressing the START 
sivltch. Before the test begins, a very low level, 
short duration check is mude to confirm that a pro- 
IK'T rcsrwnse signal is I’^elng received \yy each channel 
s|iecified in the setup. If ihe 'loop-check'* is unsaiSs- 
faeiorv' (e. g, , faulty' charge amplifier, improjicr sig- 
nal connections, etc.} a diagnosljc message \s print- 
ed for the ojieralor and the test is not |XM*mitted to 
start, rpon com pinion of a satisfactory loop-cliock, 
the »est l>egins up»on receipt of a 
openitor, ** • v " *■ 
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BJTER PARAMETERS l«YES 0«MOt i 
XMPUT l«KYai>> 0«RORt i 

1 TEST tot Tai23 

2 HEADlMGt T-2123 2-AXiS 

3 aAMUVlOTHt 2000 
aAi^DWlUTKs 2048* 

4 RESOLUTION 64/ 128/2 $6/ 512 t 256 
fREOUENCY INCREMENT! 8*000 

REFERENCE SPECTRUM! 

5 INITIAL SLOPE# Da/0'**T! 12 

ALARM LIMIT DB t 3 
ABORV LIMIT 03 ! 6 

6 FREQUENCY HZ.i 50 

LEVEL GSQR/HZ*! .01 
ALARM LIMIT 03 ! 3 
abort limit OB t 6 

7 FREQUENCY HZ*! 200 

LEVEL GSQR/HZ*: .01 
alarm limit DB ! 3 
ABORT LIMIT 03 t 6 

8 FREQUUCY HZ* t 400 

LEVEL JSQR/HZ.! ..I 
ALARM LIMIT 03 : 6 
abort LIMIT 03 ! 9 

9 FREQUENCY HZ*! 600 

LEVEJL GSQR/HZ* t «1 
ALARM LIMIT OB : 6 
ABORT LIMIT DB : 9 

10 FREQUENCY HZ* ! 0 

10 FINAL SLOPE# Oa/OCT: -18 
GRMS X 6.511 

11 LOW LEVEL# -DBt -20 

12 LEVEL INCREMENT# DB: 3 

13 START-UP TIME SECi 10 

14 SHUT-OOVN time SECi 1 

15 TEST TIME HRS# MIN# SEC* 0# 5# 0 

16 AUTOMATIC INCREASE# l«Y£S 0sNO* 0 

17 CONTROL CHANNELS* 1 

18 AUXILIARY CHANNEL* 2 

19 ACCEL SENS MV/G* 

CH 1* 10 

CH 2* 10 

20 DRIVE CLIPPING 1«Y£S# 0-NO: 0 

21 ALARM LEVEL GRMS* 8 
ABORT LEVEL GRMS* 10 

22 LOOP-CHECK MAX ORIU£( VOLTS) * .*1 

CORRECTIONS 1«YES 0«NO t 0 

Fig. 20. Dialog for Entry of Typical Test 
Parameters 


Initial equalization to the reference spectrum takes 
plac3 at the low level specified during setup, e.g. , 
in Fig. 20, -20 dB with respect to the full GRMS level. 
When the test has equalized at low level, it may be In- 


creased to the full GRMS level either manually, In 
steps specified during setup (e.g. , 3 dB) or automat- 
ically, either after a specified time at low level or by 
pressing the INCREASE switch The test shuts down 
automatically if permitted to run for the specified 
duration at full level. The test also shuts down auto- 
matically if an :’^^ort condition occurs, and may be 
shut down manually at any time. Whenever a test is 
shut down, the test level is decreased smoothly in the 
specified shut down time and a post test message is 
printed, providing a record of the tebt duration at full 
level and of why it sh down, as whown in Fig. 21. 


Conditions which cause an automatic abort include: 

1. Initial open loop condition on any ?. lected channel. 

2. Loss of signal on any selected channel. 

3. Control signal exceeding the GRMS limit spec- 
ified, 

4. G-/Hz level of any line of the control spectrum 
deviating from its corresponding reference 
speculum line by more than the specified abort 
limit. 

5. An operator abort command. 

6. Closure of any of 16 external contacts. 

Z-AXIS RANDOM ON T1012 
ABORTED 

CTL SIGNAL LOSS I 
TEST TIME HRS# MIN# SEC: 0#to#22 
MONITOR FRAMES* 5 


T-2123 Z-AXIS 

ABORTED 
OPERATOR HALT I 

TEST TIME HRS# MIN# SEC: 0# 0# 58 


Z-AXIS RANDOM ON T10I2 
COMPLETED 

TEST TIME HRS# MIN# SEC* 0#L#0 

MONITOR FRAMES* 25 

Fig. 21. Examples of Post Test Documentation 

During and after a test, calibrated and annotated dis- 
plays (see Figures 5, 6, and 9 through 16) of the con- 
trol, monitor (averaged control), reference, error, 
and auxiliary spectrums are available for test ana- 
lysis and documentation. 
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CONCUSSIONS 

The control techniques, hardware implementation, 
and operation aspects of digital random vibration con- 
trol system? have been discussed. When compared to 
analog systems, today’s digital random control sj^s- 
terns offer faster test setup, more resolution, wider 
dynamic range, more flexible frequency coverage, 
repeatable test execution, more complete test docu- 
mentation, and easy, economical conversion to other 
control and anlysis tasks. 
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SINE-SWEFt* TESTING USING DIGITAL CONTROL 


Dr. Alfred G. Ratz 
Ling Electronics 
Anaheim, California 


A review is given of the theory of swept-sine vibration testing usirj“ an 
elect rodynamic shaker. The objective is to jhow how tho coirtrol rcauire- 
meats can be fully satisfied with a digital system, based on the use of 
a minicoR^Jter. The only equipment needed over that for control of either 
random or pulse testing is a dig itally~coat rolled sine oscillator. 

Autcjiatic corop»^ession speed, specimen Q, sweep rate, tracking filter 
feedback, control to the average or the maximum of a number of simtltaneous _ 
acceleration signals, specxmen protection, are among the topics covered. 
Finally, examples are given i»aowing the performance of a practical syste-^ 
with single and multiple return accelerometers. 


INTRDDCCriON 

The purpose of this paper is to describe 
the theory and practice of putting together an 
automatic sine-test control system for an 
electrodynamic exciter, based on the use of a 
minxcoQ^uter. The discussion centers around a 
particular embodiment for reference and example; 
however, the emphasis of the discussion is on 
fundamental techniques. Some of the techni- 
ques are basic to swept-sine testing for any 
type of control; other techniques are truly 
unique to digital systems. 

There are a nusd^er of important 
advantages in using a digital system for sine 
(xntrol. The system pemits te^t parameters 
to be automatically set up ahead of time with- 
out tying up the exciter and associated 
equipment. There are no switches or potentio- 
meters to be adjusted by the operator, as 
compared with the considerable number of 
critical controls that must be tuned in a 
sophisticated analog system. System setup 
errors and debugging times are eliminated. 

The level of operator skill is considerably . 
reduced. Automatic monitoring of pertinent 
test results is easily arranged. It is also 
easy to expand a system to meet new needs by 
simply modifying the software programming. 

In summary, then, the advantages 
possible with digital control are: increased 

productivity (that is, we can quickly go from 
one type of test to another) , operational 
cost savings (due to the Increaseu productivity, 
as well as the ability use operators with 
less skill) , repeatability (the identical 
program held in paper tape, on disc, OJ: using 


some other digital storage medium, can be 
used over and over so that the identical test 
is called oat), test integrity (no human error 
in setting up the test, no controls to be 
inadvertently manipulated during a test) , and 
safety and protection. In the latter case, 
the^computer can watcii over a variety of 
parameters with an unflagging attention, a 
reaction speed, and a sensitivity impossible 
with a practical analog system. 

Sine control using a digital computer 
has tended to take second place to the random 
control or the shock pulse contro) of an 
electrodynamic exciter. There are several 
reasons for this , some quite valid. Accurate 
shock rxintroT , for example, has always proven 
very difficult to achieve, prior tc the 
introduction of the digital control system. 
And, the cost and complexity of analog random 
systems have made them a target for newer 
technology. But, ir terms of the performance 
requirements and the design problems to be 
handled, the control of sine testing presents 
a level of technical difficulty that is in 
many ways more severe than that experienced 
with ra.idom control. For example, if one 
examines the various digital random control 
systems presently described in the literature, 
one finds that, in dynamic range of equaliza- 
tion, hardly any match the performance of 
high-quality analog systems. The degradation 
in performance becomes glcssed over by digital 
considerations; but with swept-sine control, 
it is difficult to achieve a similar camou- 
flage. as the sinewave sweeps over the 
freqaency bard, as each specimen resonance 
is dealt with in its turn, the quality of 
control remains at all times clear. 
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discernible# and disguised* 

Sine control capability ot a digital con- 
trol system is usually added as an adjunct to 
random control; and the system serving as the 
basis of the present discussion is indeed 
obtained by the sin?)le expansion of an auto- 
matic random control system. The only 
ancillary eouipment added to give the system 
sine-control capability is a digitally-con- 
trolled oscillator. 

THE ELECTHDDYNAMIC EXCITER 

Eiectrodynamlc vibration exciters are 
important features of the environmental test 
facilities associated with the design, 
development, and proof-testing of aircraft, 
missiles, and satellites. The principle of 
the electrodynamic exciter follows that of 
the permanent-magnet loudspeaker. But, the 
en^hasis of design is quite different. 

Frequency response is secondary. Considera- 
tions of force, power, length of stroke, etc., 
are much moire in^ortant. Consequently, the 
transfer characteristic of the exciter is far 
from flat. 

Interactions between load and exciter 
can nxxiify the frequency response, intro- 
ducing sharp peaks and notches across the 
spectrum. Resonances are often introduced 
by the properties of any fixture used to 
mount the specimen to tlie exciter, as well as 
by the transfer function of the specimen 
itself. 

Thus, the unpredictable properties of 
the frequency characteristic of the exciter 
and its load make mandatory the use of clc*ied- 
loop control, otherwise, it is not possible 
to carry out representative testing. Any 
attempts simply to drive the exciter open loop 
can only lead to chaotic and indeterminate 
results. For efficient (^ration, feedback 
is needed. 

Swept-sine is the basic sine-test mode. 
With swept-sine testing, the electrodynamic 
vibration exciter is driven by a sinusoidal 
signal, whose frequency is slowly swept back 
and forth between two limits, at a pre- 
determined (usually logarithmic) sweep rate. 

At each frequency, the vibration level must 
adhere precisely to a value established by a 
programmed input curve (level versus fre- 
quency) . 

Other forms of sine tests; dwell, 
phase control, multiple-accelerometer control, 
etc, , are simply modi fi cations and adaptations 
of swept-sine, and are best discussed as such. 

BPSIC THEORY OF SINE CONTROL 

We now discuss the parameters needed for 
a successful sine vibration system. The 
object is to apply a sweeping sinewave with 


a specified amplitude profile to the test 
object. 

Compression Control . A real-time 
amplitude feedbadc control cannot be used for 
several reasons. The complex nature of the 
response of the test item would require a 
specific cor^Lex equalizer in order to achieve 
a stable closed loop. Test objects vary 
tremendously in their respcnse. A different 
specific equalizer would be required for each 
test package. Most important is the fact that 
there is insufficient bandwidth through the 
control system (exciter, driving as^lifier, 
accelerometer electronics, etc.) to permit 
a stable reliable closed-loop instantaneous- 
signal control system. 

Amplitude control is acdiieved by apply- 
ing the te<*nique often called coc^ressor 
control, or AVC. The essence of the metiiod is 
to obtain a measure of the vibration level 
(nns or peak) at the exciter, and then to 
attenuate or amplify the amplitude of the sine 
drive to eliminate any discrepancy between 
the measured and the required level. 

TL -2 required attenuation gain range 
over «diich the sine level must be adjustable 
is of importance. A typical test c^ject can 
be expected to have resonances in its trans- 
fer-function response . with values of Q up 
to and exceeding 60 dB. Even so, the typical 
required accuracy is 2% of specified level. 

The control must thus cover over a 1000:1 
range with 2% accuracy of point at every 
attenuation level in the active dynamic 
range. 

Compression With Random Control . The 
range of an^litude control required to handle 
spectral fluctuations in the response of the 
specimen is demonstrably much greater with 

sine testing than for random. With wide- 
band random, resonances with extremely sharp 
values of Q become invisible to the control 
system, and do not need to be dealt with at 
all. Also, one finds the limitations of 
control performance obscured behind such 
factors as the statistical jitter on the out- 
put spectral plot. 

With random control, the numbers 
representing the spectrum deal with it 
basically as a set of linear quantities. One 
finds little awareness of how seriously the 
spectral dynamic range is limited in 
practice by the significant figures actually 
available to repre ent the individual spectral 
lines. Even so, marginally successful 
random control systems are reported. But, 
there is no such lucky situation for the 
designer of the sine system. The sinewave 
stands out clear and naked. If its amplitude 
were represented linearly, three decimal 
digits would be needed to express accuracy 
and three to handle ranging, for a total of 
six decimal digits. Tlu excessive number of 
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decinuil oiglts causes the designer to search 
for an alternate method of representing signal 
strength, to the linear method that passes as 
acceptable for random control. 

A logarithmic format can handle a large 
d/namic range vith constant percentage 
accuracy. Hence, a log control scheme is the 
one used in all successful sine*^control systems. 
A logarithm control schere is also essential to 
solve the problems inherent with compressor 
control stability, as it tarns out. 

Servo Control Loop . Figure 1 shows the 
basic elements of a coepressor control loop, 
whether digital or analog. In Figure 1, the 
reference (ref) is the required vibration 
asplitude on the test package at the present 
value cf frequency. The difference is obtained 
between the reference (ref) and the level at the 
test object. It is used to generate an error 
signal. The integrator smooths out the error 
signal and applies a gain command to the 
attenuation anplifier. The anplifier 
adjusts the level of the drive signal and then 
applies it to the test package. The detector 
rectifies the sinewave representing the 
resultant acceleration, and generates the new 
control signal, thereby closing the feedback 
loop. 



Figure 1 - Basic Servo Loop 


Compressor Speed . An important para* 
ineter of the compressor control loop is its 
speed of response. The requirements for wide 
dynamic range and high accuracy force the 
feedback loop gain to be high. The complexity 
of the test specimen requires that the com- 
pressor-loop time constant be the predominant 
time constant in the servo. 

The speed of response of the control 
loop is set by two considerations. The control 
loop must be fast enough to satisfactorily 
hold the desired acceleration level as the 
drive frequency is swept through the test 


resonances. The control loop must be slow 
enough to permit smoothing of the ripple from 
the rms detector, thereby permitting the 
control to accurately follow the average level 
of the sinewave. (If the control tends to 
follow some instantaneous component, it 
introduces distortion into the sine drive 
signal.) 

The goal of the control loop is to 
obtain a 2% control accuracy. For design, 
let us allow a maximum of 1% detector ripple. 

Wc can now calculate the compressor speed for 
a \rticular frequency. The integrator 
used to reject ripple does so at the rate of 
6 dB/octave. If we assume that the detector 
is a full-wave rectifier or a square-law 
detector, then the ripple frequency is twice 
the fundamental frequency. The average 
detected value is roughly equal to the 
air^litude of the unsmoothed ripple. Therefore, 
the compressor loc: bandwidtli must be 
decreased to the point where the drive 
frequency is larger than the loop bandwidth 
by orders of magnitude. For example, for a 
fundamental drive frequency of 1000 Hz, the 
compressor feedback loop bandwidth must be 
no wider than 5 Hz to 10 Hz. 

Feedback loop bandwidth of the com- 
pressor system serves as a measure of the 
ability of the compressor to control the 
anplitude of the accelerometer signal under 
the many fluctuating dynamic conditions that 
occur when the frequency is swept. It is ex- 
pressed in radians per second (rad/s) , 

The expected error in control as the 
sine drive is swept through a resonance can 
be predicted by calculating the control error 
E(dB) , 

E = (i) 

The quantity v is the amplitude slewing-rate 
capability of the sys**<^m, measured in dB/sec. 
Anplitude slewing rate is also called ”com- 
pressor speed." From Equation (1) we see 
it also is a measure of the bemdwidth of the 
control loop as a servomechanism. The 
relationship to the servo bandwidth can be 
deduced as follows. Suppose that f(t) 
renresents the output signal, and that f(t) 
has its envelope changing with time (A is a 
constant) : 

f(t) » A • e"“c^ (2) 

Taking logarithms and differentiating yields 

(w^ = 2TTf ) : 
c c 

H, , dL ( t ) . a. 68 u « 54.5 f (3) 
dt c = 

There are two main operational para- 
meters influencing the required range of 
values for the slewing rate. These parameters 
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are (1) the Q of the resonanv:e that the drive 
frequency is presently sweeping through and 
(2) the freqxaency sweep rate. The effect of 
these in establishing an acceptable sle%r-rate 
is discussed below* 

SYSTEM EieODIMENT 

The above discussion concerning the 
theory is applicable to either analog or 
digital systems (the latter our subject here) . 

Block Diagram * in the digital applica- 
tion, it is efficient that the portion of the 
diagram cf Figure 1 whicdi consists of the 
reference, summer, and integrator, should be in 
software in the computer* The sine generator 
and the variable-gain amplifier are inplement- 
ed with devices handling the signal in analog 
form, but controlled digitally. The frequency 
of the sinewave locks to the nuid>er digitally 
supplied to the generator* The gains of the 
variable-gain amplifiers are set by numbers 
sui^lied by the computer* The detector can 
be handled either in an analog or in digited 
form* In the present s^'stem, it is handled 
digitally* 



Figure 2 - Compressor Loop - Digital System 

Figure 2 shows a simplified block 
diagram. Shown are the sampling and quan- 
tizing devices required by the con^uter: an 

A/D converter and multiplexer (MUX) , as well 
as a D/A converter. The detector and its 
filter are achieved in software. 


signal must have their sasgile rates smoothly 
and Infinitely programmable without losing 
timing accuracy. In view of the complexity 
of synchronous san^ling, it is worth exploring 
the non-syn^ronous san^ling of the cycle of 
duration t, Fig\ure 3: suppose we are 

attempting to find the rms value of the sine- 
wave by taking the mean square value of N 
samples of the sinewave over one period. The 
sang>ling rate is sucdi that there are (N+m) 
samples per cycle of the wave (m < 1) * The 
first sample occurs %ihen the argument of the 
sinewave has the value 6. 



Figure 3 - Sampling: MUX and DATA 

The computer forms the average, , 
given by : * 

NS^ = sin^O + 2i^N+m))+ sin^ (6 + 4TT^+m) ) 

+ *+sin^ (6 + 2w (N-l^(N+m) ) 


(4) 


By trigonometric manipulation, the equation 
for NS“ can be rewritten: 


2NS^ = N cos 26 


im*h) * 4ir 
sin N + m 



+ sin 20 

. 2irtm+l) 
sin • sin 

N+m 

2Trm 

N+m 


“"‘to ’ 


(5) 

If m “ 0, 9 =0, becomes the true rms, - 

Also, 2NS » N 
o 


Non-Synchronous Sampling * The method 
of "detecting” the return signal is to take in 
a frame of data, t seconds long; t is very 
nearly equal tn the period of the fundanental 
of the sinusoid. If p is the sampling rate, 
it is desirous to keep p synchronous with f, 
the frequency of the sinewave. Since f shifts 
in very fine steps, this weans a continuous 
and ^ery fine updating of p; the A/D con- 
verter and multiplexer dealing with the input 


and so S = 0.707 
o 

Examination of Equation (5) with N >> 1 
indicates that it is certainly bounded by, 
for all m: 

2NS^ a N ± 1 (6) 

a 
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flence, « S^(l ± **N) (7) 

rc control to better than ±0*3 dB, we miist keep 

N > 16 - 2^ (8) 

Thus# as f shifts, we do not need to shift p, 
u^itil Equation (B) is violated. When p is 
I'hanqed, we must change the cut-off frequency, 
associated with the anti-aliasing filter at 
. le HUX input. We keep the -0.3 dB (cut-off) 
requency, F, of a given filter, as the lowest 
/alue possible, satisfying 

f j< F (9) 

When f starts to exceed F, the filter is 
stepped to a larger value of F. 

Suppose we are engaged in multiple- 
cceleration control; i.e., suppose we are 
ontrolXing to the rms average or extremal of 
acceleration signals. We mu«' bring in N 
ine histories during the t*' x. Then, the 
ultiplex rate is raised to . A typical 
.aximiun allowable value for Mp set by hardware 
j)erfonnance is 125,000 H^. Hence, for eight 
i^nput channels, the highest values permitted 
for f is given by setting 

128 F = 125,000 

Thtjs f £ 3000 Hz (10) 

If we need to test at frequencies higher than 
loco Hz, 6 channels, we singly use more than 
one cycle of the sinewave. 

number of input channels that can be 
handled by (MUX A-D) in K cycles of the 
operating frequency, f, is given by: 

Mj^ F - 8000 K (11) 

Frxn the diagram of Figure 3 we note that the 
cc. .trr.)l algorithm is handled in the interval 
of time, Tj^, minus the time x. x is now equal 

to ( K/f) seconds. Since the need for using 
multiple cycles occurs only at the higher 
f rer. uencies , little restriction is imposed 
whe; we are forced to go to multiple cycles. 

i »ata Samples . Fc-. each input frame of 
the return signal tl jrefore, we compute an 
rms value: for t! purposes of control, the 

ms value so obtained constitutes a single 
data point. . assessing the return signal 
fiom the t. .ter, the actual time history is 
of no valuo. The information sought is the 
rms le^ 1. The servo speed must be geared to 
tlie • 'te of change of the information in the 
sxquai, not that inherent in its time history. 

H ce, the servo speed is actually slow 
enough so that we do not nave the computer 
‘»xamir ) the return sig'iil continuously. It 
can saiple one cycl'i of the wave (for the 
time duration, x) allowing a time, T^, be- 
tween frames to carry out the subsidiary 


confutations in the computer (Figure 3). Thus, 
the control is a sampled-servo system with 
a sanfling period, T^: 

T « X + T (12) 

L O 

During each time, T^, the coBfuter (CPU) com- 
putes the level of the return signal from the 
sanf les derived during the previous ingesting 
period, x. It then derives a new control 
signal to be sent to the gain control register 
of tn output variable-gain amplifier to 
correct for any deviation in the return level 
from that specified by the program. 

The output amplifier, under CPU control, 
carries out the logarithmic gain-control 
function discussed above. 

Specimen Q . One parameter that has a 
large effect cn the control performance is 
the maximum value of the specimen Q expected. 

In analog systems, the operator must set in 
via a front panel switch the expected specimen 
Q. With the digital system he mtist likewise 
indicate the Q to be expected via instructions 
to the CPU. 

It is to the operator's advantage never 
to call for a larger Q than actually occurs. 

If the specimen has no sharp resonances, the 
operator calls for the iow-Q parameter for 
the control. For higher Q values, he uses the 
medium-Q parameter and for very sharp 
resonances, he picks high-Q. The low-Q con- 
dition is defined as having Q- values between 
zero and 50; the medium-Q condition is de- 
fined as having Q values between 50 and 200; 
the high-Q case 2 is having Q values larger 
than 200. 

The digital system user has one import- 
ant advantage over the user of an analog 
system. If he selects the wrong value of Q, 
the control automatically corrects his 
selection for the actual value of Q encounter- 
ed. This feature is discussed belc^. 

Control Loop Action . Figure 4 
illustrates the control loop in further detail. 
x(t) is the return sigi al received from the 
accelerometer system and its signal condition- 
er. The gain of the system input amplifier, 
denoted by G^, is set to a value consistent 

with the desired value of acceleration at the 
table, and the need to keep the rms level of 
x(t) near to a full-scale value. We must 
allow for the peak/rms ratio of the signal, 
plus a need to recover from a large, sudden 
signal at the table. It is indicated that 
the A/D must therefore be set to handle 
signal peaks at least twice the rms. The A/D 
converter has 12 bits, with 1 bit for sign, 

1 bit for the average/peak allowance just 
discussed, X bit to allow for the fact that 
can only be set in 3 dB steps. The rms 
level, therefore, approximates the 9t}i bit 
in the A/D converter. (NOTE: For what follows, 

Figure 10 provides a logic flow diagram.) 
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The first step after ^A> conversion is to 
cotnpute the mean square of the signal; i.e., 
to compute over the cycle rf interest, 

Y = = 1/n E X. (13) 

1 ^ 

Y is compared with the reference value, 
auid an error number, e, produced. 

d = y - Yj^ (14) 

The reference level, Y_, roust be calculated by 
R 

the computer for eacdi point on the test pro- 
file, as determined by the control program. 

For each anplitude value required by 
the program, there is an optirouro value for the 
gain G. , bringing the rros controlled level at 
point as close to the optimum value as the 
gain steps of the amplifier allow. To com- 
plete the calibration of the input circuitry 
to 1%, it is necessary for the program to 
adjust the reference level, Yj^, to compensate 
for the size of the minimum increment of gain 
in the amplifier G^. But, changing the value 
of Yj^ has the effect of altering the loop 
gain, and this effect must be removed: a 

gain-normalizing factor, G, is introduced; 

e = G d (15) 

where G = K/Y 

R 

The quantity, e, is the final normalized 
error number. Error numbers e are summed, 
one number for each iteration of time. 

The sum produces the number, e' i 

e' = E e (16) 


The number is then scaled by a factor R as 
required by the compression rate being used. 
The combined effect of the loop time, T , and 

the loop gain (as finally set by R) determines 
the con^^ression rate. After scaling by the 
quantity R, the quantity E' is converted to 
the final asq>lltude level connand (output 
signal) . 

D « re' (17) 

There is then a final word conversion to 
convert D to a gain-conmand word to set the 
gain of the output an^lifier. 

The output an^lifier is essential ii. 
performing the compression. It is in5>ossible 
to set gains merely by shifting data words 
digitally in the CPU. There is no way to do 
this and maintain control accuracy over the 
wide dynamic range needed. There is only one 
>ray to achieve coo^ression feedback having 
4_he sort of characteristics proven essential 
over the years, for effective sine-servo 
control, and that way employs a wide-ranging 
an^lifier. 

Compressor. The compression rate is 
affected by the Q of the specimen. If a high- 
Q specimen is being tested, it is necessary 
to slow the control loop down- Otherwise, 
the loop becomes unstable. In setting up a 
test, we set in a value of Q assumed to be 
the largest that we can expect across the 
frequency band of interest. Thus , we are in 
effect assuming a fixed Q across the band; 
and so the conpression rate (i.e., it. 

Equation (3)) to be permitted at each present 
value of the frequency, f, must be proportion- 
al to f; the proportionality constant must 
obviously be inversely proportional to Q. 

Thus, the compression rate, can be written: 

il; — f/Q (18) 

The anplitude control loop as a servo- 
mechar.ism has a frequency response, f^, that 
is related to ^ by Equation (3). Obviously, 
f^ is affected by T^ ; 



The gain R also affects f : 
c 

f ^ R (20) 

c 

Thus, at any given frequency, 5n going from 
low-Q to high-Q conditions, the control-loop 
speed must be reduced; this can be accomplish- 
ed by a combination of two techniques: 

(a) Increase the "data sample" time, T^. 

(b) Decrease the loop gain, R, so that it 
tal’.es more control iterations to correct 
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By a mutual interaction of and n, 

%ie force if, throu^ f^, to follow the curves 
of Figure 5. It is convenient with a binary 
CPU to alter the value of t for every octave 
dtange in f. thus# for frequencies between 


f^ and f « 2fj# Tj^ 


T . ; then for f 


between f « 2f^ and f • 4f^, we dtange Tj^ 
to Tl * *j Tj^j. At the lower end of the 
frequency range, can be kept constant 

for a given Q; a typical relationship is; 


df/dt • q « bf (2: 

Hence, log <f/f ) « bt (2: 

e o 

If W is the sweep rate in octaves per minute. 


If Tg = the s%n:ep time in minutes (the time to 
sweep from the lowest frequency, f^, to the 
highest, f^) 


In moving from low*Q to medium-Q control, R 
is reduced by the factor 2"^; to go to high-Q 
control, R is reduced further by a factor of 
2 -\ 



Figure 5 - Automatic Compression Rates 


then, WT = log,(f /f ) (25) 

s Z u z 

A digital system sweeps by incrementing 
the data in the frequency register of the 
oscillator at the correct rate. Thus, to 
carry out a given frequency sweep, the 
oscillator is addressed from time to time by 
the con^uter and fed new "frequency" words. It 
locks the frequency of the output sinewave to 
the latest "frequency" used. In order that no 
resonances are missed, we must keep 

Df/f < a/Q (26) 

where a is a nuinber less than unity. Let Df 
be the frequency increment. Combining 
equatiors, we can compute the corresponding 
time increment, Dt, A suitable value for K 
is 0.3. Hence, 


At the higher frequencies, relationship 
(21) cannot be followed beyond a certain 

frequency, since T C2Uinot be reduced below 
L 

the time duration needed for the computer to 
carry out the basic work that must be done 
between "data samples." Thus, T^^ at the high 

end of the spectrum is fixed and R must be 
adjusted to compensate. 


As mentioned above, the digital system 
can out-perform an analog sine control system 
in one important respect. The user selects 
whether he intends to run at the low-Q, 
medium-Q, or high-Q condition: if he chooses 

wrongly (i.e., if he selects too low a Q) , 
the computer notes the instability of the 
control, and adjusts automatically to 
compensate: the adjustment is done in steps 

of 2”^, until stability is reached. (Of 
course, an upper limit for Q is also implied 
by the sweep rate selected, so that a sweep 
slow-^down is required as well as a compres- 
sion-rate slow-down; if such is indeed the 
case, though it could adjust the sweep rate, 
the system is not presently programmed to 
do so, but to abort! It would seem that 
changing the sweep rate is something the user 
should do.) 

Logarithmic Frequency Sweep . The law 
usually followed in sine-sweep testing starts 
slowly at low frequencies and increases the 
rate of frequency change at higher frequen- 
cies; it follows the logarithmic law: 


Dt = 25/WQ 


It is reasonable to set 


An acceptable value for Dt is 0.1 seconds. 
Based on this value, we can set up the follow- 
ing table for the maximum value of W (i.e.. 


TABLE 1 

Sweep Rate Vs Specimen Q 


Condition 


Lo Q. 
Medium Q 
High Q 


Assumed 
Maximum Q 


5 oct/minute 
1.25 oct/minute 
0.25 oct/minute 


The operator of the system must restrain 
to fit the above table. Once he has chosen 

High Q, Lo Q or Medium Q, he must not then 
select cun incompatible W^. The above table 
summarizes basic physical limits imposed by 
the servo, and does not arise because either 
a digital or analog system is used. The size 
of the corresponding frequency increment is 
given by: 
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Keeping Dt at 100 milliseconds independent of 
the present value of puts Df ^t 


Df « 


Wf 

8?0 


(30) 


Thus# the computer arranges for finer fre* 
quency increments# %fhen W is less. This is 
as it should be# since a reduced value of W 
can be taken as in^lying the possibility of 
encountering a resonance with a high Q and 
its corresponding narro«rer bandwidth. The 
resonance therefore needs smaller frequency 
increments to be excited prcperly. 


Tracking Filter Mode . An important 
method of operation for any swept-sine 
vibration control system is to control on the 
level of the fundamental of the return 
acceleration signal# rather than on the 
conposite nns. With analog systems# a tracking 
filter must be added to the control loop to 
filter the acceleration signal. With the 
digital syste«n# no extra hardware is required# 
the filtering being done in software. The 
only hardware cdiange is to use one channel 
of the input MUX to bring in a frame of 
the full-scale oscillator drive signal# before 
it is «q>plied to the output attenuator, in 
other words, the oscillator output# F^^Ct) , is 
treated like another return signal# and is 
multiplexed with the accelerometer signal or 
signals (for multiple accelerometer control) . 
Let 


F^(t) = A sin wt (31) 

From the input sanqples taken during the time# 
T# the con^uter develops from each sample of 
Pj^(t), a second signal P 2 (t)- 

Fjtt) = ^ (32) 

T^it) and p 2 (t) are used as reference signals. 
The phase difference between F^(t) and the 
data samples is not important provided the MUX 
has low jitter. Each sample of the input 
signal or signals being tracked is multiplied 
by the corresponding value of Fj^(t) and of 
its derived value F 2 (t). For each input 
signal# x(t)# therefore# the coc^uter conroutes 
t)ie data gathered over each period# t# two 
data values# Y and z^s 


two variables Y and Z. T^ seconds later# the 
(j-fl)-st sample pair obtained. 


Y^ and are fed to two identical low- 

pass filters within the software of the 
coR^uter. These low-pass filters are obtained 
by digital filter techniques. Their outputs# 

Y and Z # resp>ectively « are# of course# 
oj oj 

derived from the stream of values Y^ and Z^. 

They are combined to obtain a single value# 

X. : 

3 

+ Z^. (34! 

3 03 03 


replaces 


in the servo block diagram. 


In working out conditions for the 
tracking filter# the filter 3-dB bandwidth# 

B# must be wide enough to be useful. Very 
narrow values for B are much easier to design# 
but their utility and efficienc/ are limited 
to very special situations# since they force 
the user to restrict the sweep speed that can 
be permitted to very low values. 

Since vibration phenomena are usually 
constant-Q in nature# it is ^-ry useful to 
have B proportional to f. Thus# as f increases 
B must be widened. This cannot be achieved 
indefinitely as f increases: there is a 

value of f beyond which it is not practical 
to widen B. 

For sejn^o-loop stability# we must 
assume that 


BT » 1. (35) 

L 

This is a basic servo law that can never be 
violated by digital or analog equipment. 

Hence# corresponding to the conditions for 
which Equation (21) applies# et seq.# we can 
draw up a similar set of conditions. We set 
Bj^ as the absolutely maximum value of B that 
can be permitted. Smaller values of B can 
always be obtained. 

The value of B is influenced by R and 
H 

Tj^. The following table gives values of B 
that are representative of what can be outain- 
ed with acceptable values of Tj^ and R(T^ and 
R# of course# are affected by f) : 


Y^ « / x(t) • P^(t) dt 
Z^ - x(t) • F 2 (t) dt 


(33) 


‘j 


and Z^# of course# are coirputed by sunma- 

P(t) 


tion.of the appropriate x(t) * F(t) products 
fo> the data sanples gathered over the time # 
Thus# Y and Z are the j-th samples of 


TABLE 2 

Acceptable Filter Bandwidths 


f 


B 

250 

Hz 

2 

Hz 

500 

Hz 

4 

Hz 

1000 

Hz 

8 

Hz 

2000 

Hz 

15 

Hz 
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Rftferring bade to the discussion on com- 
pression speedy the CPU can use the largest 
value of the ratio (f/B) to be used in a ^iven 
test# to set the Q value for determining \{>. 

MULTIPlJB INPUTS 

For large fragile specimens# controlling 
to a signal computed from the outputs of a 
nunber of control accelerometers is very 
desirable. Two methods of control are: 

(a) Averaging# (b) Extremal selection. With 
method (a)# as has been discussed above# all 
the return signals are multiplexed over the 
frame duration# t* The sum of the sequence 
of all KUX sauries# normalized by the number 
of input diannels# gives a mean square value 
for control. 

With method fb) # control is switched to 
the return signal having the largest (or small- 
est) value. Either the rms or peak value can 
be used. In either case# each input must have 
its level confuted separately during the time# 
T# and the decision carried out as to ^e 
controlling channel immediately after the 
frame has passed. 

The speed of selection is a popular (and 
prob 2 d>ly critical) criterion in the specifying 
of a selection system for vibration control. 

It has been a big factor in the specification 
of analog selection circuits and should be 
equally is^rtant with digital systems. It is 
important# therefore# to establish the maxim*im 
useful selection speed. It turns out the 
selection speed depends on the speed of 
response (conqpression rate) of the control 
loop. It is no use selecting faster than the 
loop can respond# and extremely fast selection 
does not improve system performance. An 
estimate of the switching speed needed is 
important. Suppose we have two input signals 
and are operating in a selection mode. 

Suppose the first signal is in control# with 
the second signal initially set to zero 
amplitude Suppose now the second signal 
increases xn amplitude# the change following 
a step function: the level of the controlling 

signal is passed instantaneously. The 
pattern of change just described represents 
the worst possible case. 

Assume# now# that the selection is 
instantaneous by some miracle of design* The 
control loop commences to reduce the error 
introduced at the co*i^ression rate. The 
loop reduces at the maximum rate of 8.68 dB/ 
time constant. Three loop time constants 
reduce the error to 5% of its Initial value. 

Now suppose instead of the step-function 
change in the selection# we had the selection 
process follow a romp function so as to take 
a time equal to one loop time- constant to 
reach the new error signal. Then the loop 
would require four time constants to reduce 
the error# instead of three for the case of 


instantaneous response# (a slow-down of 33.3%). 
This seems a small decrease in performance# 
even for the extreme operating conditions 
surmised. In practice# with slower sigaal 
changes and slower congress ion speeds# the 
slow-down is probably never larger than 10%. 

Besides speed# signal weighting is 
another in^rtant factor affecting selection 
design. Signal weighting can be dealt with 
in a number of ways. The easiest method in- 
volves har<hfare only: it is sinply to adjust 

the accelerometer sensitivity as fed into 
the input normalizing amplifier# if such exists 
Otherwise# weighting can be done in software 
by introducing scale factors to the control 
loop via the CHI. 

ALARM/ABORT IMPLEMENTATION 

Alarm/abort limits are very useful 
features# easy to program into digital- control 
systems. When the acceleration level 
deviates from the desired by more than the 
alarm limit# the operator is notified. If it 
deviates by more than the ah t limit# the 
test is t<‘*-^\nated. 

It I ccessary uo eliminate false 
alarms oc 32 aborts. This can be done by 
having alarm <xnd abort controls programmed 
to act in some manner as the following: 

(a) A running count# n# is to be kept 
separately for alarm and for abort. 

(b) Each period of time# T^# the return 
signal is examined for exceeding the 
alarm and the abort levels. If it does 
exceed a level# the relevant count# n# 
is increased by unity. 

(c) Each period of time# T # the return 
signal does not exceed a level# the 
relevant count, n# is decreased by unity. 

(d) The val\:ie of n can never go negative. 

(e) When n reaches the value N^# the relevant 
level is considered to have been exceeded# 
and abort action is initiated. 

N^ largely depends on experience. A value of 
two or three for N^ has proven satisfactory. 

The sine control system is uniqiie in 
one way. Its chief cause of consistent 
alaxm/abort conditions is the encountering of 
a hlgh-Q situation with too fast a conpression 
speed set in by the operator. Thus, when 
the control loop sees an alarm condition# in 
addition to notifying the operator# it alters 
the con^ression Icqp gain by decreasing the 
value of R# as is described above. A 6 dB 
change is made. 

After an alarm condition has been 
reached# and is reset to zero# the action can 
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commence again. It the alarm is again exceed- 
ed, R is again adjusted by 6 dB. Thus, ic is 
possible for the system to eventually correct 
for a serious error in the original specifica- 
tion of Q. Thus, unlike any analog system, 
the system has the ability to adaptively 
adjust the control parameters, depending on 
the specimen Q. 

SPECIMEN PROTECTION AND SAFETY 

In common with ra’idom control, the 
specimen must be guarded against malfunction- 
ing of CPU, amplifiers, MUX's, etc. Most 
malfunctions are caught singly by their 
effects of driving the signal outside alarm 
and abort limits- 

One iijportant item to be uniquely 
guarded against is failure of the accelero- 
meter. Several algorithms are possible, 
but it has been proven that the best scheme 
here actually copies the ana^g protection 
device; it measures the rate of ch^mge of 
the signal level received from the accelero- 
meter. If the increment in the time incre- 
ment T is negative and exceeds a previously 
established value, A^, the excessive 

demand is assumed to be due to lack of 2 Ui 
accelerometer (return) signal. The serious- 
ness in setting the value of A depends on 

M 

the fragility of the specimen, and the 
program must permit the user to select A 
at will. Whenever an abort is called for? 
the con^uter attenuates the drive signal to 
the shaker smoothly and gracefully to zero, 
preventing any possible sericus acceleration 
transients from occurring. 

MIXED SINE/RANDOM 

Mixed sine/random testing involves 
superimposing sine and random tests on a 
specimen simultai\eously. The tests are 
programmed individually, and the first task 
of the control is to separate the return 
sine signal from the return random so that 
each type of signal can be used in its own 
control loop. 

The sine component is easily pulled 
out of the random by using the tracking- fil- 
ter mode. Narrowing down the bandwidth, 6, 
of the tracking "filter,” is easy to do with 
a digital system, and makes mixed sine/random 
easy to accomplish from the sine poxnt of 
view. 

The randcai signal is much more 
difficult to extract from the sine with a 
digital control system: the best way to ex- 

tract che random seems to be as follows: 
il) The composil:e return signal is fed to 
the FFT process; and a Fourier transform 
of the composite signal is produced, (2) At 
the present value of the frequency, f, of 
tiie sinewave, the FFT line is much higher in 


level because the energy of the sinewave is 
superio^osed, (3) However, the computer is 
aware of the lines affected, since it has 
control of f, and for the line or lines in 
the neighborhood of f, the random closed 
loop control ceases to function; instead, the 
PSD level of the output signal is set locked 
to tliat worked out before f swept into the 
region of the lines* Once f sweeps past the 
line(s) in question, the closed-loop control 
is again free to operate for them. 

The principle problem to be worked out 
with mixed sine/random is the adjustment of 
the signal levels through both the sine and 
random systems (the latter on a PSD basis) , 
to allow for the presence of large, inter- 
f erring signals. Fortunately, it has been 
proven practical to get excellent control for 
a considerable range (35 dB or more) in the 
value of the ratio of (random rms)/(sine 
rms) , the range being centered at 0 dB. 

EXAMPLES 

Out of a multitude of test results, a 
few have been selected to illustrate the sort 
of equalization achievable binder sine control. 

The first test shows the equalization 
of a peak/notch filter used to simulate an 
actual shaker and specimen. See Figure 6. 



Figure 6 - Peak Notch Correction 

A notch xs set close to 500 Hz, and a peak 
close to 990 Hz: the Q of the notch is 61, 

and the Q of the peak is 72, for a total 
p/n spread of 73 dB. The sweep rate is 0.5 
oct/min. The resulting control shows a 
maximum error of 0*6 dB. Further reduction 
of the error is possible, If system para- 
meters are adjusted. For example, the MUX 
sample rate is set by Equation (8) ; if it 
is increased, relative to (1 /t), the error 
can easily be cut in half. However, the 
present values for the parameters would seem 
adequate for practical control. 
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NOW we turn to a more interesting 
example. In this example, eight accelerometers 
are placed along a resonant beam, and mounted 
on a shaker. Control is affected as follows: 
control to the maximum level of inputs #A, 

#3, »4, #5, #7, 18; input «A is the 

average of acceleration signals #1, #2; 
inputs #3 through #6 are the signals out of 
the corresponding accelerometers, weighted 
by a factor 0.25; inputs #7, #8, weighted by 
a factor 0.125. 

The program specified is: 0.75 G, 

10 Hz to 100 Hs; constant DA from 5 Hr to 
10 Hz; sweep rate: 4 oct/min. Figure 7 

shows the input program, with channel #? 
signal superisqposed; we see at What frequency 
channel #2 in effect takes over all by itself# 
channel #1 yielding too weak a signal at those 
frequencies to contribute significantly to #A. 

We also see where the signal from 
channel #2 does not control alone; either the 
average or one of the other signals controls 
in the regions where signal #2 drops away. 
Figure 8 shows the return signal from 
accelerometer #8; it must reach 6 G to take 
control, and only achieves this level when 
the drive frequency is 36 Hz approximately. 
Figure 9 plots all t ight input acceleration 
signals, weighted. 



Figure 7 - Channel No. 2 


CONCLUSION 

Swept- sine vibration testing provid^'s 
many interesting and unique aspects of 
interest in the techniques of shaker control. 
Control is essentially non-stationary in 
nature, as opposed to random control or 
shock pulse control. Adaptive servo 
techniques have proven practical and useful. 

The control of swept-sine vibration 
testing involves handling of much wider 
dynamic ranges than are required with random 


or pulse control. And there is no opportunity 
here of glossing over the dynamic- range per- 
formance. The sinewave stands out clear 
and unencumbered. 

The versatility of digital sine control 
is obvious. It will only be a matter of time 
when its potjntial will be fully realized 
for handling testing situations where analog 
devices have performed marginally at best; 
for exan^le, it will not be long before it 
will be reported that a practical digital 
system has been developed for the simultaneous 
phase and amplitude control of multiple 
shakers connected by a common resonant load, 
the goal of considerable desire and effort 
over the past years. 




Figure 9 - All Channels - Eight-Channel Control 
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TIME HISTORY SYNTHESIS FOR 
SHOCK TESTING ON SHAKERS 


D. 0. Smallwood 
Sandia Laboratories 
Albuquerque, New Mexico 


Several digital methods now ■available for matching shock 
spectra with oscillatory type transients are reviewed and 
the merits of each discussed. The methods include WAVSYN. 
SHOC, sums of decaying sinusoids, fast-sine sweeps, modu- 
lated pseudo random noise, classical pulses, and the method 
of least favorable response. While the last method is no<* 
strictly a method for matching shock spectra, it is dost . 
related. The discussion is limited to procedures for sy? 
thesizing the ivaveform. Additional parameters (in addit 
to the shock spectrum) which can be specified to limit th 
classes of functions available for performing a particulai 
I test are discussed. 


INTRODUCTION 

The need for an osci 1 1 atory-type 
pulse to represent certain field en- 
vironments has been recognized for 
several years. During this period, 
it was also recognized that electro- 
dynamic or electrohydraulic shaker 
systems were well -suited for repro- 
ducing these pulses. With the devel- 
opment of digital control techniques 
it became possible to reproduce very 
complex time histories. Therefore, 
considerable effort has been expended 
in developing these methods. 

A comi.«on tool for measuring the 
character of a transient has been the 
shock response spectrum of the event. 

As a result, many current methods are 
based on producing a transient whose 
shock spectrum matches a specified shock 
spectrum. Although considerable con- 
troversy sti"" exists regarding the 
value of the shock spectrum, its use 
is very common and will continue for 
years regardless of the outcome of the 
controversy. Because of its projected 
use and the fact that procedures which 
exist can result in quite different 
waveforms, it was felt that a review 
of the available techniques would be 
useful . 

Both analog and digital methods 
are being used. However, the purpose 
of this paper is to review the digital 


techniques. The problem addressed in 
this paper is: given a specified shock 

spectrum how can a waveform be genera- 
ted which will have the same (within 
some tolerance) shock spectrum. The 
companion problem of how to reproduce 
the synthesized time history is left 
for another paper. 

Before the time history synthesis 
can be discussed, the limitations placed 
on the waveform by the shaker system 
will be reviewed. 

The types of transient vibration or 
shock pulses which can be accurately 
reproduced on both electrodynamic and 
el ectrohydraul i c exciters are very 
much dependent on the physical limita- 
tions of the exciters. These limita- 
tions are listed in Table 1 and are 
briefly discussed here. 


TABLE 1 

Exciter Limitations 


Limi ta tion 
Number 

Initial 

Final 

Maximum 

1 

x(0) = 0 

)i(T) = 0 

Limited 

2 

x(0) = 0 

x(T) = 0 

Limi ted 

3 

x(0) = 0 

x(T) = 0 
(Electro- 
dynamic) 

Limi ted 


The initial and final acceleration 
and velocity of a transient must be 
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zero for both electrodynamic and elec- 
trohydraullc exciters. As with any 
type of testing machine, maximum 
attainable values of acceleration 
and velocity are limited. Accelera- 
tion Is actually limited by the force 
capabilities of the exciter. 

Flexures In electrodynamic exciters 
generate restoring forces which 
return the exciter table to Its 
originating position (defined as zero), 
and limitation 3 holds. This is not a 
requirement for el ectrohydraul 1 c sys- 
tems; however, by imposing this limita- 
tion one can take advantage of both the 
forward and return portion of the 
stroke to generate the required tran- 
lent. This effectively results in 
doubling the displacement capacity of 
the exciter for generating transients. 
Thus for the purposes of this discus- 
s’on limitation 3 will also be consid- 
ered a limitation for an electro- 
hydraulic exciter. 

The ' .itial slope of a lightly 
damped shock spectrum is related to the 
velocity and displacement changes re- 
quired. This is illustrated in Figure 
1. Here is an example of waveforms for 
different values of the initial slope 
(i.e., the slope of the shock spectrum 
as the frequency approaches zero on a 
log-log plot). Note that the shaker 
limitations will require that only 
pulsf:s of the last type M>12) can be 
reproduced on a shaker system. If the 
initial slope of the specified shock 
spectrum is not greater than 12 dB/ 
octave there will exist a lower fre- 
quency below which the shaker displace- 
ment and velocity limits will not per- 
mit the spectrum to be matched. If the 
test item responds like a rigid body at 
all frequencies less than this lower 
limit an argument can be made that it 
is not important to match the spectrum 
in this range. 

Techniques using oscillatory 
pulses whose shock spectra have an 
initial slope of S = 1 , as for example 
decaying sinusoids, have been moderately 
successful because the shaker system 
will act as a high- pass filter removing 
the low frequency energy from the wave- 
form. This, combined with the flexure 
restoring force, will force the velo- 
city and displacement to return to 
zero. The manner in which the accelera- 
tion-time waveform is distorted to re- 
move the velocity and displacement 
change will be characteristic of the 
individual shaker system used. This 
makes it difficult to predict the 
velocity and displacement waveforms 
until after the test is run. The 
velocity and displacement waveforms 


TYP£ INITIAL SLOPE ACaiLRATION VELOCITY PIS PLACEMENT 

M > dB^Ct 
S * No<)dim. 





TIME 

Fig. 1 - Initial Slope of the Shock 
Spectrum for Some Common 
Waveforms 


are important as the shaker system 
places limits on the magnitudes which 
can be reproduced. 

The shock spectrum at very high 
f-equencies is identical to the peak 
amplitude of the input time history. 
Again the shaker limitations will 
determine if this level can be repro- 
duced. In many real specifications 
the shock spectrum is not specified 
to a high enough frequency to deterniin* 
the required peak input. In this case 
some flexibility is available in 
synthesizing the pulse. 

With th.s introduction to the 
limitations of shaker systems and to 
the characteristics of a shock spectrum, 
the methods used to generate a time 
history whose shock spectrum will match 
a specified spectrum can be discussed. 

PARALLEL FILTER METHODS 

The first method discussed is the 
parallel filter method. This technique 
is a direct digital implementation of 
older analog methods [ij. Using this 
technique tfie waveform synthesis and 
shaker equalization are not separated 
but accomplished together with the test 
item mounted on the shaker. Using this 
method the required shock spectrum is 


if. 
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broken into regions (frequency ranges, 
typically 1/3 octaves). For each region 
the peak response is specified, along 
with a basic waveform whose energy is 
concentrated in the same frequency range. 
Any of the basic waveforms discussed 
later could be used. The most common 
ones are WAVSYN and decaying sinusoids. 

In the earlier analog methods the basic 
waveform was the response of a bandpass 
filter to an impulse. An initial ampli- 
tude for each basic waveform (one for 
each region or frequency band) is chosen 
and a composite waveform is formed by 
summing the waveforms. This time his- 
tory is then applied to the power ampli- 
fier and the control accelerometer re- 
sponse is measured. The shock spectrum 
of the response is calculated and com- 
pared with the required spectrum. This 
information is then used to modify the 
amplitudes of the basic waveforms and 
the process repeated as illustrated 
in Figure 2. 



Fig. 2 - A Flow Chart for a Parallel 
Filter Iteration Method 


A variation of the procedure out- 
lined in Figure 2 is shown in Figure 
3. fhe principal difference between 
this variation and the previous one is 
that the frequency components are added 
one at a time starting with the lowest 
frequency and working up to the highest 
frequency. When all the components 
have been added the complete spectrum is 



Fig. 3 - A Modified Flow Chart for a 
Parallel Filter Iteration 
Method 

checked. If the spectrum is not within 
spe.cifications the procedure is started 
over at the lowest frequency using the 
previously determined amplitudes as 
starting values. Using this technique 
only one frequency component at a time 
is modified. As a result, the method 
is more stable than the previous method 
(more likely to converge) as the inter- 
action between components is less like- 
ly to cause problems. The advantage of 
the previous method is that if a stable 
solution is found, fewer pulses are re- 
quired :o equalize the system. 

Th? only difference between the 
digital implementation and the older 
analog methods is that a larger variety 
of basic waveforms can be used, the 
shock spectrum is computed digitally, 
and more sophisticated and automated 
an.litude modification methods can be 
used. The advantages of these methods 
include: The methods are relatively 

easy to implement and use, as details 
of the time history are not controlled. 
Quite large values of the shock re- 
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sponse spectrum can be generated as 
the shaker is not required to be equal- 
ized over sharp notches in the frequency 
response function. This tends to con- 
centrate the pulse energy in those fre- 
quencies where the shaker responds well. 
Also since equalization is good only in 
a broad sense, a large variety of shock 
spectra can be matched. The disadvan- 
tages include: Since details of the 

time history are not controlled, little 
is known about the velocity and dis- 
placement requirements. Because the 
shaker is not equalized with a fine 
frequency resolution, all frequencies 
may not be adequately tested. If the 
shock spectrum is analyzed with a finer 
resolution than the initial range, 
deviations from the required spectrum 
may be found. 

TIME HISTORY SYNTHESIS 

Several methods have been developed 
to synthesize a time history which will 
match a specified shock spectrum. The 
synthesis is usually done on a digital 
computer. It is assumed that this time 
history can be reproduced on a shaker 
system using methods originally devel- 
oped by Favour and LeBrun [2] and will 
not be discussed in this paper. Several 
of the methods have checks Included to 
determine the suitability of the tran- 
sient for shaker reproduction. Re- 
producing the pulse on the shaker system 
is independent of the procedure to 
develop the waveform. 

O'le method will be discussed in 
detail. The remaining methods will be 
discussed more briefly as the procedures 
are similar. The first method discussed 
is sums of decaying sinusoids (OS). The 
other methods discussed include HAVSYN, 
Shaker Optimized Cosines (SHOC), fast- 
sine sweeps, modulated random noise, 
modification of field time histories, 
classical pulses, and least favorable 
responses. 

Sums of Decaying Sinusoids 

It has been recognized for years 
that many field environments can be 
adequately represented by sums of 
decaying sinusoids, and as a lesult 
several authors Ts, 4, S] have sug- 
gested their use to match shock spectra. 
The usual basic waveform is given by 

g_(t) « A e"^m^m^sin w„t taO 

mm m 

“ 0 t<0 

The basic waveform is shown as 
Figure 4, and the normalized shock 
spectra, for several values of c are 
shown as Figure 5. Several waveform:; 


are added to synthesize a time history 
to match a complex shock spectra. 


i 



Il^« «i 


Fig. 4 - Decaying Sinusoid Acceleration, 
Velocity and Displacement Char- 
acteristics 



Fig. 5 - Normalized Shock Spectra for a 
Single Frequency Decaying 
Sinusoid 
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However» this waveform does not 
meet the requirements of a zero velocity 
and displacement change. If an attempt 
is made to reproduce this waveform, it 
will be slightly distorted by the shaker 
system removing the velocity and dis«^ 
placement change. Since the exact dis- 
tortion will be a function of the shaker 
system, it becomes difficult to predict 
the velocity and displacement waveforms. 

Two modifications of the basic wave- 
form have been suggested. The author 
[ 3 ] suggested that a time-shifted, highly 
damped decaying sinusoid be added for 
velocity and displacement compensation. 

It is also possible to introduce time 
shifts to Improve the appearance of the 
waveform. Allowing for time shifts of 
the individual components, the waveform 
becomes 

X(t) - ^U(t-T )A 

• sin 






•sin ^^(t+x) 


(c! +1) S 


T . » y Vm_ 

, i t 

I)'-*' 

The magnitude (A^) and the shift 
(t) of the velocity and displacement 
compensating pulse are fixed by the other 
parameters. A plot of a typical waveform 
is shown as Figure 6. 

Prasthofer and Nelson [ 4 ] suggested 
velocity and displacement compensation 
by adding two exponential pulses and a 
phase shift to the decaying sinusoid to 
give 


X(t) - E 


sln(u t ♦ e)> 

i m » 


“c " “n/ VI- 

ut 


a * w^/2if 
c 




(a-b) [(c-a)^ 


(a-b) ^(c-b)*' + 
I(b-;)^+ (a-c)^+w^ 


0 • tan~^ tan~^ 

m 

-1 *^ifi 

The first two terms are added for 
velocity and displacement compensation. 
The phase shift, 6, is added to force 
the initial value to zero. 

A plot of a typical modified pulse 
is shown as Figure 7. 

lOi — I 1 1 



Fig, 6 - Modified Decaying Sinusoid 
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Fig. 7 - Modified Decaying Sinusoid 


A set of curves can also be gener- 
ated which gives the peak response 
ratios (the ratio of the peak response 
to the input amplitude at the frequency 
of highest response) as a function of 
the decay rate (c) of the decaying 
sinusoid and the damping factor (n) of 
the shock spectrum (Figure 8). 

Example Using Decaying Sinusoids 

Let us assume a test requires that 
a pulse be reproduced on an electro- 
hydraulic shaker system. The required 
shock spectrum is given in Figure 9. A 
typical field time history which resulted 
in the required composite shock spectrum 
is also included (Figure 10) for compari- 
son purposes. The shaker limitations 
are: peak acceleration 50 g's o-p, peak 


velocity 120 In/sec o-p» and peak dis- 
placement 8 in p-p. The pulse will be 
composed of sums of decaying sinusoids 
compensated using the method suggested 
by the author. 

The spectrum cannot be matched at 
the lowest frequency because the slope 
of the spectrum is less than 12 db/oct. 

In actual fact, a 16 ft/sec velocity 
change is associated with the spectrum. 

It will be assumeo that the test item 
is not likely to have a..y resonances 
below 6 Hz or above 115 Hz and the spec- 
trum will be matched from 4 Hz to 115 
Hz. The first step is to assume all 
the decaying cinusoids will act in- 
dependently. Using Figures 5 and 8 the 
amplitudes and decay rates for the 
components can be picked. A 4 Hz 
pulse is picked to raise the low end of 
the spectrum. A large decay rate is 
picked to keep the spectrum reasonably 
flat in this region. Further fre- 
quencies are picked at the peaks in the 
spectrum. 

The initial estimates (Trial 0) are 
shown in Figure 11 and Table 2. From 
displacement considerations the compen- 
sating frequency is picked at 1 Hz with 
a decay rate of 1.0. The shock spectrum 
of the composite pulse is shown on 
Figure 12. The spectrum is low at 4 Hz 
due to the interaction of the 1 Hz pulse. 
The spectrum is also low from 4-12 Hz, 
high from 12-90 Hz, and low above 90 Hz. 
As expected the assumption of independ- 
ence was only approximate. An iterative 
procedure is then used to modify the 
component amplitudes to match the spec- 
trum at the frequencies of the decaying 



Fig. 8 - Peak Response Ratio for a Decaying Sinusoid 
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Fig. 9 - Example Specified Shock Spectrum 



TIME (sec) 


Fig. 10 - A Sample Time History 

sinusoids. One possible iteration 
procedure is outlined in Figure 13. 

The results are shown in Figures 12 and 
Table 3 as Trial 1. Note that the solu- 
tion does not converge; i.e., the ampli- 
tudes of the 36 and 47 Hz components 
were reduced to zero, but the spectrum 
was still too high. The high level is 
caused by the interaction of the lower 
frequency components. For the next 
trial (Trial 2), the decay rates of the 
26 Hz, 36 Hz, 47 Hz, and 60 Hz compo- 
nents were lowered to reduce the inter- 
action. Components at 8 Hz and 80 Hz 
were added to raise the level near those 
frequencies. The level of the shock 
spectrum at 112 Hz was also raised to 
improve the match at that frequency. 
Again an iteration procedure was used 
to determine the component amplitudes 
which would result in matching the re- 
quired spectrum at the component fre- 
quencies. As before the amplitude of 
the 47 Hz component was reduced to zero. 
This time the spectrum was about S% 
high at 47 Hz, The results are shewn 
In Table 4 and Figure 12. At this point 
it was determined that the match was 


suf:*icient and the pulse parameters were 
within the shaker capabilities. 



Fig. 11 - Composite .»l*uck Spectrum 

Formed by Several Decaying 
Sinusoids 



Fig. 12 - Shock Spectrum for Several 

Composite Waveforms Composed 
of Sums of Decaying Sinusoids 

In this example an acceptable match 
was achieved In two iterations (Trial 0 
is not counted, as it is usually skipped, 
but was Included for illustrative pur- 
poses.) The time histories of the 
acceleration, velocity and displacement 
are included as Figures 14-16. The 
parameters in Table 4 can now be used to 
define a time history to be reproduced 
on the shaker system. 
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Fig. 13 « Flow Diagram for Picking 

Decaying Sinusoids to Match 
a Given Shock Sp^'^trum 


TABLE 1 - Parameters for Trial 0 



FREQ 

ZETA 

AMP 

TAU 

I 

(Hz) 


(g) _ 

(sec) 

1 

1 

1.0 

-5.37 

-0.141 

2 

4 

0.5 

7.1 

0 

3 

14 

0.2 

8.8 

0 

4 

26 

0.2 

8.8 

0 

5 

36 

0.1 

5.7 

0 

6 

47 

.05 

3.4 

0 

7 

60 

.01 

2.3 

0 

8 

112 

.02 

8.1 

0 


TABLE 

3 - Parameters 

for Trial 1 


FREQ 

ZETA 

AMP 

TAU 

I 

(Hz) 


(9) 

(sec) 

1 

1 

1.0 

-5.71 

-0.138 

2 

4 

0.5 

8.67 

0 

3 

14 

0.2 

10.1 

0 

4 

26 

0.2 

8.52 

0 

5 

36 

0.1 

0 

0 

6 

47 

.05 

0 

0 

7 

60 

.01 

2.22 

0 

8 

112 

.02 

8.47 

0 


TABLE 4 - Parameters for Trial 2 


I 

FREQ 

ZETA 

ANP 

TAU 

(Hz) 


(9) 

(sec) 

1 

1 

1 

-5.74 - 

0.141 

•} 

4 

0.5 

6.47 

0 

3 

8 

0.3 

4.45 

0 

4 

14 

0.2 

11.1 

0 

5 

26 

.08 

3.51 

0 

6 

36 

.05 

0.646 

0 

7 

47 

.03 

0 

0 

8 

60 

.005 

1.78 

0 

9 

80 

.01 

2.73 

0 

10 

112 

.02 

9.93 

0 

Acceleration Range 

+18.1, -9.0 g 

Velocity 

Range 

+63, -75 in/s 

Displacement Range 

+2, 

1. -4.9 in 



TiMt 


Fig. 14 - Acceleration of a Composite 
(Trial 2) Composed of a Sum 
of Decaying Sinusoids 



Fig. 15 - Velocity of Three Composite 
Waveforms 
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Fig. 16 - Displacement of Three 
Composite Waveforms 

WAVSYN 

A method developed for the U. S. 
Army is described by Yang, [6* 7] This 
method is called WAVSYN» The basic 
waveform is given by 


g (t) * A^sin(2irb t)sln(2irf 1*) 
uj mm m 

for OStsTji^, 

-0 for t>T„, 

where ■ N b , 
m mm 

T • l/( 2 b_), 
m m ' 

and N is an odd integer, 
m 


A specified shock spectrum is then 
matched by summing several of the basic 
waveforms. 




The basic adjustable parameters 
are the number of pulses summed L, the 
frequencies f«, the number of half 
cycles ana the amplitudes of each 
component A^. 

The delay t,- is used to improve 
the appearance or the time history 
(i.e., make the time history look more 
like field data) and does not seriously 
change the shock spectrum. Normalized 
shock spectra for N » 3, 5, 7, and 9 
are shown in Figure 18. Note that the 
peak in the shock spectra occurs near 

^m* 

The shape or magnification can be 
modified by changing N^. The frequency 
of the peak in the shock spectrum can 
be controlled with f^ and the magnitude 
of the whole curve can be moved up and 
down with Aj^. Detailed procedures [6, 

7] are described for matching a wide 
variety of shock spectra using the tech- 
nique. An advantage of the technique 
is that the velocity and the displace- 
ment changes are zero. This is essen- 
tial for accurate reproduction on a 
shaker system. 


01 damping 


The waveform for N = 5 is shown as 
Figure 17. 


x(t) ■ A sin Ci>t sin Nut 





0.1 0,2 0.3 0.5 1 2 3 5 

NORMAL IZLD FREQUENCY (V ) 


Fig. 17 - WAVSYN Time History, N « 5 


Fig. 18 - WAVSYN Shock Spectra 
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For comparison» the same shock spec- 
trum as In the previous example was 
matched using MAVSYN. Eight components 
were used as listed In Table 5. 


TABLE S - WAVSYN Components 


No 





1 

4 

3 

1755 “ 

0 

2 

6 

3 

4.82 

0 

3 

9 

3 

7.37 

0 

4 

12.5 

7 

3. 

0 

5 

27 

5 

4.3 

0 

6 

33 

7 

0.94 

0 

7 

70 

11 

2.44 

0 

8 

112 

19 

8.86 

0 


Maximum Acceleration 20. 

9 g*s 


Maximum Velocity 

1 07 

In/sec 


Maximum Displacement 7. 

10 In 


The acceleration* velocity* and dis- 
placement time histories are shown in 
Figures 19, 15, and 16. The shock spec- 
trum is shown in Figure 20. 



Fig. 19 - Acceleration Time History of 
WAVSYN Composite Pulse 



Fig. 20 - Shock Spectra Comparing Decay- 
ing Sinusoids, SHOC* and WAVSYN 


SHOC 


A method developed by the author [8] 
for matching shock spectra with an 
oscillating transient Is called the 
SHOC (Shaker Optimized Cosines) tech- 
nique. The basic waveform for this 
technique Is given by 

c: (t) = A e*^ni“m'^cosw^ t -B_cos^ irt 
m m r» ^ — 

0<tST/2 

=0 t>T/2 


g„,(-t) . 
where B 

m 


ejt), 


ro m ' 

the pulse duration 


The first term is an exponential 
decaying cosine and the second term 
(a cosine bell) is added to force the 
velocity and the displacement changes 
to zero. The waveform is symmetric 
about the origin and builds up and then 
decays. 


The basic waveform is illustrated 
as Figure 21. The normalized shock 
spectra for several values of c are 
shown as Figure 22. As in the case 
of WAVSYN, a parameter is available for 
modifying the shape of the curve (c„,)> 
the frequency of the peak of the shock 
spectrum (to ) and the level of the curve 


As in the case of WAVSYN, several 
of these components are added to syn- 
thesize a waveform which will match a 
complex shock spectrum 


L 

x(t) * E 

m=l ^ 


Also, as in WAVSYN, J.he velo- 
city and the displacement changes are 
zero. 

For comparison, the same shock 
spectrum as used in the two previous 
examples was matched using SHOC. Four 
components were used and are listed in 
Table 6. 
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Fig. 21 - A SHOC Pulse 



Fig. 22 - SHOC Shock Spectrum 


TABLE 6 - SHOC Components 


FW 

(Hz) 


Amplitude 

x(sec) 

4.0 

1.0 

10.0 

0.54 

13.0 

0.1 

2,5 


32.0 

0.2 

3.0 


li 2 . 0 

.03 

9.8 


Maximum Acceleration 23,7 


Maximum Velocity 69.5 

in/sec 

Maximum Displacement 6.5 

in 


The acceleratlont velocity* and dis^ 
placement time histories are shown In 
Figures 23» 15* and 16. The shock spec- 
trum Is shown In Figure 20. The three 
previous methods are also compared in 
Table 7. 

Note that for this example the 
characteristics of all three pulses are 
quite similar* ai7d the three methods 
would be expected to have a similar 
damage potential* although this has not 
been firmly established. 





00 ai a? a) a 4 ns ao 

TtM 


Fig. 23 - Acceleration of a Composite 
SHOC Pulse 


Fast-Sine Sweeps 


The first analytical procedure 
developed for matching shock spectra 
was a technique using fast-sine sweeps* 
which lasted from a few hundred milli- 
seconds to several seconds. The early 
work In this area was done by Crum and 
Grant [9] and later expanded by several 
authors Including Roundtree and Fre- 
berg tlO]. The basic procedure is to 
assume an acceleration time history of 
the form 


x(t) * A(t) sin G(t) 


A particular form is then chosen for 
for A(t) and 0(t) which includes a 
limited number of variable parameters. 

The derivative of 0(t) is the Instanta- 
neous frequency of the time history. 
Procedures are then given for picking 
the variable parameters to match the 
specified spectrum. Crum and Grant 
chose 

A(t) = A (a constant) 


0(t) * -2 ttN’ 
0(t) * -2 itN» 


In^l- 


ln|l+ 



Increasing 

frequency 



decreasing 

frequency 


3 ’ 
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TABLE 7 - Comparison of Decaying Sinusoids* WAVSYN, and SHOC 

Peak Ratio of Peak Velocity Displacement 

Input Duration'*' RMS Duratlon^^^ Response to Reoulrement Requirement 
(e) (ms) (ms) the Peak Input (In/sec) (in) 


Method 

Peak 

Input 

(fi) 

Duration 

(ns) 

HAVSYN 

21 

300 

SHOC 

24 

120 

Sums of 

Decaying 

Sinusoids 

19 

210 


(1) The duration Is defined as the time beginning when the waveform first reaches 
10-percent of Its peak value and ending when it decays to 10-percent of its 
peak value for the last time. 

(2) See Appendix A 


where the variable parameters are the 
constant for A(t), the starting fre- 
quency fg, and the effective number 
of cycles at each frequency N'. The 
effective number of cycles at each 
frequency (N*) is based on the ratio 
of the shock spectra with a Q of 5 
and 25> fp was chosen as the lowest 
frequency of interest, and A was 
determined by the required amplitude 
of the shock spectrum. The form chosen 
by Roundtree and Freberg was a more 
general formulation given by 

« Rf(t)^, 

= 2rrf(t), e(0)-0. 

Where the variable parameters are B, a, 

fo. R. Y. 


As mentioned earlier, the deriva- 
tive of 0(t) is the instantaneous 
frequency f(t). B represents the rate 
of change of the amplitude function, 
A(t), with respect to the instantaneous 
frequency, f(t). The starting fre- 
quency Is fq. The Initial value for 
the amplitude function A(t) Is a. R and 
Y are used to control the rate of fre- 
quency change and the sweep duration. 

For example, if B = 0, the sweep ampli- 
tude is held coniitant, and If y = 0 , 
the sweep rate will be linear. The 
sweep will be exponential for y - 2. 

Roundtree and Freberg present 
graphs for picking values for the para- 
meters based on the requested shock 
spectra computed at two different values 
of damping. 


The methods of Crum and Grant have 
been used for several years. The 
methods of Roundtree and Freberg are too 
new to be well established. The ad- 
vantage of these methods is that they 
produce pulses which are well suited 
for reproduction on shaker systems. It 
has also proved practical to match 
shock spectra at two different values 
of damping. The principle disadvantage 
is that the techniq*'es produce pulses 
which do not resemble (in the time do- 
main) many (if not most) field events. 
Another limitation is that the tech- 
niques have been highly developed for 
matching spectra which can be repre- 
sented by a straight line on a log-log 
shock spectrum plot and adapt poorly 
for matching spectra with other shapes. 

A typical time history for these methods 
is shown in Figure 24. 



Fig, 24 - Time History of a Fast-Sine 
Sweep 
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Modulated Random Noise 

Authors In the seismic field have 
long recognized the somewhat random na- 
ture of seismic records. As a result, 
numerous proposals [11. 12] have been 
made to derive a random process which, 

when multiplied by a suitable window, 
will provide an adequate simulation of 
seismic events. The attempt Is to de- 
rive a waveform which will have the 
statistical characteristics of a seismic 
event. Since one of the tools used to 
characterize the seismic event Is the 
shock spectrum the event, the modu- 
lated random noi:.e methods produce 
pulses which match a specified shock 
spectrum. However, it Is Important to 
recognize that these methods only match 
the shock spectrum in a probabilistic 
sense. 

Bucclarelll and Askinaza [13] 
proposed that these same Ideas could be 
used to simulate pyrotechnic shocks by 
using an exponential window. 


x(t) = g(t)n(t) 


where g(t) Is a deterministic function 
of time which characterizes the trans- 
ient nature of the event. 



The function n(t) is a representative 
of a stationary, broadband, zero-mean, 
noise process with a spectral density 

It Is then proposed that the func- 
tion S||(u) be chosen such that the aver- 
age Fourier amplitude spectrum of x(t) 
will be equal to the average Fourier 
amplitude spectrum of the field environ- 
ment. It Is shown that the following 
approximation will accomplish this 
purpose. 


ElP(u)P*(u)] aS^(u)/2B 

where 

, E[P(u)P*(ui)] 


is the expected value of the field 
Fourier amplitude spectrum. 


The above equations seem reasonable 
for the following reason. Sn(u>) Is a 
stationary process with a fixed amount 
of energy per unit of time. As g becomes 
larger, the duration of the transient 
becomes shorter and must be In- 

creased to maintain a given energy In 
the transient. 

Tsai [14] then proposes that the 
random variation In the shock spectrum 
be removed In the following manner. 

Pick a sample x(t). The shock spectrum 
of j((t) Is then computed. In areas 
where the shock spectrum Is low, add 
energy to the waveform by adding sine 
waves to n(t). In areas where the shock 
spectrum is high, filter n(t) with a 
narrow-band notch filter. Compute the 
shock spectrum of the modified x(t) and 
repeat the process until the desired 
shock spectrum Is achieved. The final 
modified x(t) will then be used as the 
test input. A procedure could probably 
be worked out to accomplish the same 
goal by modifying the Fourier amplitude 
of n(t) without resorting to adding sine 
waves and filtering [15,16]. 

This procedure looks quite Inter- 
esting, but as far as the author knows. 
It has not been tried in the laboratory. 

Modification of Field Time History 

Workers In the seismic field have 
suggested a method very similar to 
the modulated random noise technique. 

For this method the original sample 
(x(t)) 1s an actual field time history 
(In the seismic field an earthquake 
record). The Fourier content of the 
waveform Is then modified to force the 
shock :>pectrum of the time history to 
match a specified curve. 

Classical Pulses 

It may be desirable to reproduce 
a classical pulse (half-sine, terminal 
peak sawtooth, etc.) on a shaker system. 
In general these pulses are of Type 3 
(Fig. 1) and cannot be reproduced 
without modification. The pulses can 
be easily corrected to a Type 5 pulse 
by adding a bias pulse of equal and 
opposite area at the end of the pulse. 
Pulses of this type can be reproduced 
with only the distortion Introduced 
by the improper displacement boundary 
conditions. It is always possible to 
Introduce a bias pulse which will also 
meet the displacement boundary condi- 
tions. For example, consider a terminal 
peak sawtooth pulse biased with a square 
wave as shown in Figure 25. 
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Least Favorable Response Techniques [17] 

The basic assumption of the least 
favorable response techniques Is that 
the magnitude of the Fourier spect<*um 
has been specified. Since the un^'amped 
residual shock spectrum Is related to 
the Fourier modulus. It could also be . 
specified. It Is then shown that If the 
frequency transfer function between the 
mechanical Input and the response of a 
critical point on the test Item (not the 
transfer function of the shaker system) 
can be characterized by 


Fig. 25 - A Modified Terminal Peak 
Sawtooth 

It Is easy to see that the velocity 
change will be zero If 


H(w) • |H(w)| ^ 


BT = 1/2 A Tg. 


If an expression is written for the 
acceleration of this waveform, and 
it is double Integrated to solve 
for the displacement. The displace- 
ment boundary conditions can then be 
used to show that a zero displace- 
ment change can be achieved if 

T -T 2T2 

■ I T" 


then the peak response of the critical 
point on the structure will be maximized 
by the input 


X(t«) * X 


Thus two parameters are available 
Ti) to control the velocity 
and displacement charge of he compo- 
site waveform. A siiilar analysis 
can be done for other bias pulses 
and classical wavefo’-ms. If both 
the bias pulse and the classical 
waveform are symmetrical the solution 
is always to place the classical wave- 
form in the center of the bias pulse. 

A bias pulse which has been found 
to work quite sati sfactorily is a cosine 
bell or Hanning pulse. The advantages 
of this bias pulse are: The waveform 

is zero and smooth at the ends, whi^.h 
usually improves the reproduct ' ^n. 

The energy of the pulse is concentrated 
at the lower frequencies with a well- 
defined upper frequency and a sharp 
roll-off characteristic (18 dB/oct). 

This means that the effects of the 
bias pulse will be seen only at the 
lower frequencies. The frequency where 
the effects will start to be seen can 
be controlled with the duration of the 
bias pulse. 


where Xe(o)} is the specified Fourier 
modulus, and 


dm. 


where x(t) is the required time history 
at the input to the test Item. 

The Fourier transform of the re- 
spon e or the cr tical point {Y(w)) is 
given by 


Y{o)) = X{w)H(w) 


X^(ca) I H(u)) 
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The function Y(w) Is real. This respon^^e 
can be seen to produce a maximum response 
by thinking of a transient as being com- 
posed of a sum of many sinusoids. If the 
phase angle of all the sinusoids Is zero 
at some po1n;t In time all the components 
will constructively add» and produce the 
largest possible peak. It Is also shown 
In Appendix A that this response will 
have a minimum rms duration. 

Modern methods make the computation 
of the above equatic s relatively easy. 

It 1$ proposed that che phase angle (9) 
of the frequency transfer function be 
measured In the laboratory. This func- 
tion, together with the specified Fourier 
modulus (X (u)) Is tnen used to compute 
the test input x(t). Note that the 
above technique considers the actual 
test ite 1 response, while the shock 
spectrum aoproach considers only the 
response ot single degree of freedom 
systems with known damping. The least 
favorable response method does assume a 
linear system with a well-defined criti- 
cal response. 

The method does guarantee that the 
largest possible peak response will be 
achieved (hence a guarantee of a con- 
servative test). Shock spectrum tech- 
niques cannoc make this guarantee for 
multiple degree of freedom systems. 

Several examples [17, 18] using the LFR 
method have indicated peak responses on 
the order of 1 to 2.5 times a typical 
peak response to a field event. 

This method has been used on one 
known test series [18] and the results 
are compared with the shock spectrum 
methods. As the method requires consider- 
able computations and the reproduction 
of a complicated computed waveform, digit- 
al methods are required for its applica- 
tion . 


A spin-off, from the least favor- 
able response techniques. Is to assume 
a unity transfer function (i.e., 

H(w) * 1). This means that the peak in- 
put will be maximized. If we think of 
a transient as being a finite amount of 
energy with a specified frequency distri- 
bution (Fourier modulus), the procedure 
w)ll produce a transient which will 
transfer this energy to the test item 
with the largest possible peak input and 
in a minimum amount of time. 

The input to the test item is given 
by 


::(t} 


1 A* / \ lojt j 

do,. 


Note that Xe(o)) Is a real positive func- 
tion. Therefore, x(t) will be a real 
even function. As a result, a typical 
Input defined In this manner will look 
similar to a SHOC pulse. The resulting 
Input Is Independent of the test Item 
characteristics, and hence eliminates 
the need to define the transfer function 
H(w). The only required parameter Is 
the Fourier modulus (or the equivalent 
undamped residual shock spectrum). For 
example, the least favorable Input with 
a unity transfer function for a tran- 
sient witit the same Fourier modulus as 
the transient shown In Figure 14 Is 
given as Figure 26. The duration of 
this transient Is 64 msec. One test 
series [18] Indicates that this may be 
a rt^sonable approach. 


2 > 

g 10 
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Fig. 26 - A Least Favorable Input with 
a Unity Transfer Function 

ADDITIONAL PARAMETERS (IN ADDITION TO THE 
SHOCK SPECTRUM) WHICH MAY BE SPECIFIED 

Since there are many methods for 
matching the shock spectrum, it would 
be desirable to specify additional para- 
meters limiting the class of functions 
which may be used for a particular test. 
Several parameters have been suggested 
and a short discussion of several methods 
will follow. 

L 1 mi t the Duration of the Transient 


It has been suggested that limits 
could be placed on the minimum and 
maximum allowable durations for the 
transients. It is felt that if the 
shock spectrum is matched and the dura- 
tion is comparable, the "damage*' should 
be naarly the same. For complex wave- 
forms, careful attention should be given 
as to how the duration is defined. The 
author is not aware of any test program 
where this specification has actually 
been used. 
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Require the Shock Spectrum at Two 
OlTferent Values of Damping to be 
Hatched 

Since clamping Is the poorest 
known parameter In many systems* It 1$ 
felt that If the shock spectrum Is 
matched at two different values of 
damping (for example* a Q of 5 and 25) 
the resulting transient should be a 
reasonable simulation for all values 
of damping. However, It is a very 
difficult problem to find a transient 
which will match a shock spectrum 
at /o values of damping except for 
a limited class of functions. In fact 
It Is not even clear that a solution 
always exists. It is true that a 
solution can exist* for example, a 
set of shock spectra with different 
damping values for a known time 
history. But if those spectra are 
modified (for example, smoothed, raised 
In level, or enveloped), it is not 
clear that a solution will still exist. 


Specify the Allowable Ratios of 
tne~Peak Shock ResponVe and the 
Peak input Level 

The attempt is to prevent or 
encourage the use of oscillatory- 
type input as opposed to a single- 
pulse (for example, a half-sine) in- 
put. Note that, if the shock spec- 
trum is plotted at a high enough 
frequency, the maximax shock spectrum 
reflects the peak input level and 
this specification is redundant. 


In esse*^:e, the specification will 
require a certain time history be used. 

It should be remembered that if this 
is done, only those tesc^ng laborato- 
ries which have implemented the particu- i 

lar method can perform the test. 

Thought must also be given to require- < 

ments for a successful test. One Is 
tempted to place requirements on the ^ 

fidelity of the time history repro- 
duction. However, m ..t available 
techniques for reproducing a time 
history assume a linear system. If 
the test Item is non-linear, the repro- 
duction of a given time history may be 
poor through no fault of the testing 
laboratory. 

conclusions 

It is apparent from the preceding 
discussion that a large variety of 
methods are now a'^ailable for matching 
a shock spectrum on a shaker system. 

These methods can produce quite differ- 
ent waveforms, and it is not clear 
that they will all produce equivalent 
results. It is important that a shock 
spectrum not be considered as a suffi- 
cient specification to define an en- 
vironment. rh> "characteristics" of 
the field waveform should be considered. 

The personnel who write test specifica- 
tions must be aware that several pro- I 

cedures could be used to match a speci- 
fied shock spectrum, and careful con- 
sideration should be given as to which 
procedure will (or will not) produce 
an acceptable test. 


This specification would be use- 
ful to confirm the value of the maxi- 
max spectrum at very high frequencies. 
TRW has included this requirement in 
a recent specification which required 
that the amplification o^ the response 
spectrum over the peak-g input would 
be in the range of 2.5 to 5. This 
requirement was added to prevent a 
single pulse or a very long transient 
from being used. 

Specifically Exclude Certain Methods 

If the test requester is aware 
that certain methods are not likely 
to produce a good simulation, a very 
effective measure is to simply exclude 
their use. The limit of this pro- 
cedure is to exclude all the methods 
but one. In fact, a knowledgeable 
requester can specify not only the 
method, but the parameters needed to 
define the transient. For example, 
if decaying sinusoids are specified, 
then the component frequencies, decay 
rates, and amplitudes can be listed. 


Nomenclature 

A - ampl i tude 

a - constant 

B - a constant 

b - constant, also a frequency 

c - constant 

E - expected value 

f - frequency 

g - 3 basic waveform from which a 

composite is formed, also an 
ex ^cncnt^al window 

1 - 

K - constant 

L - the number of basic waveforms 

wh*’:i" have been added to form 
a ’'iposite waveform 

M - s X. on a log-log plot in dB/ 
f . . e 

N - er of half cycles, an odd 
. -eger 
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N* - an effective number of cycles 
at each frequency in a fast 
sine sweep 

Q - quality factor, 

R > a constant 

S - nondimensional slope on a log-log 

plot, also auto (power) spectral 

density 

T - a pulse duration 
t -time 

U - a unit step function 
X - Fourier transform of x 

X - a specified Fourier modulus 
e 

X - acceleration 

X - velocity 

X - displacement 

8 - a constant 

r - a constant 

C - a decay rate 

n - damping coefficient for a shock 

response spectrum 

Q - phase angle 

T - a time delay or shift 

w - circular frequency 

Subscripts 

c - a velocity and displacement 
compensating pulse 

i - an index 

m - an index 


Superscripts 
• - complex conjugate 
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APPENDIX A 

RMS Duration of a Transient 

A transient f(t) is defined with 
a Fourier transform F((i>). The rms 
duration (D) of the transient is 
defined as, 

CO 

2 \ f 2 2 

D = £ / t |f(t)| dt (A-1) 

-oo 


where 


E =y* |f(t)|^dt. (A-2) 

E is sometimes referred to as the 
energy of the pulse. It will be 
required that E is finite, requiring 
that f(t) approach zero faster than 
1/t^, as t approaches both positive 


The rms duration is a measure of 
the central tendency of a transient. 
For example, consider a transient of 
some finite energy composed of all 
frequencies in equal amounts. An 
impulse (or delta function) will re- 
present the transient in this class 
with a minimum duration. A long-dura- 
tion, low-level random waveform will 
represent the transient with a maximum 
duration. The rms duration of several 
common transients is given in Table 
A-1. 

It can be shown fl9] that the rms 
duration is ciso given by, 

"‘'■lit /|si/ * }- 

where 

F(io) = 

If A(<d) Is specified the minimum 
rms duration Is given by 



or 

4((d) = constant 


The constant can be zero. Eq. A-4 
implies that the rms duration Is 
related to the smoothness of the Fourier 
spectrum, both the magnitude and the 
phase. The smoother the Fourier spec- 
trum, the shorter the rms duration. 
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Function 
Square Wave 

half sine 


terminal peak 
sawtooth 


triangle 


haversine 


parabolic cusp 


exponential 
single side 


exponential 


exponential 

decaying 

sinusoid 


’ TABLE A- 


The rms Duration of Some Common Transients 
Equation 
f{t) = 1 0<t<T 

* 0 elsewhere 

f(t) « sin^ 0<t<T 

» 0 elsewhere 


f(t) = t/T 0<t<T 

» 0 elsewhere 

f(t.) = 1 - 2|t|/T - J<t<J 

- 0 elsewhere 

f(t) = ^ (l ♦ cos 0<t<T 

s 0 elsewhere 

f(t) = (|iti - 1)2 - |<t<i 

* 0 elsewhere 


f(t) = 

= 0 

f(t) = 


t>0 

t<0 


f(t) = e’®“^sin «t t>0 

t<0 

= 0 

for a<<l 


rms Duration 
0.29 T 

0.23 T 
0.19 T 
0.16 T 
0.14 T 

0.11 T 
1 

5ir 

1 

v^a 

?aw 

i 
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MODAL ANALYSIS USING DIGITAL 
TEST SYSTEMS 


Mark Richardson 

Hewlett Packard Company 
Santa C1ara« California 


Modal Analysis methods based upon the measurement and post test processing 
of transfer functions in digital form are discussed. The analysis which 
shows how modal data can be identified from transfer function measurements 
is reviewed. Then various techniques of transfer function measurement 
using a Fourier ana^yze»- are presented. This is followed by a discussion 
of alternative methods for identifying modal data from transfer function 
measurements. 


1. INTRODUCTION 

An analog tracking filter can be viewed as a 
device which gives the spectral content of 
a time domain signal » one frequency at a time. 
By comparison« the digital Fourier transfonn 
can be viewed as a parallel processor which 
gives the entire frequency spectnw of a time 
signal in a selected bandwidth. 

Since the digital Fourier Analyzer provides 
a broad band frequency spectrum very quickly 
(e.g. = 100 ms to give 512 spectral lines), 
it can be used for obtaining broad band 
response spectrums from a structure which 
is excited by a broao band input signal. 

Furthermore, if the input and response 
time signals are measured simultaneously, 
Fourier transformed, and the transform 
of the response is divided by the trans- 
form of the input, a transfer function 
between the input and response points on 
the structure is measured. (Division of 
one frequency domain function by another 
is straightforward when the data is in 
digital form.) Hence, a digital Fourier 
Analyzer which can simultaneously measure 
two (or more) signals is en ioeal tool for 
measuring transfer functions quickly and 
accurately. Furthermore, since the modes 
of vibration of an elastic structure can 
be identified from transfer function mea- 
surements, a multichannel Fourier Analyzer 
with additional processing capability can 
also be used for identifying modal par- 
ameters from test data. 


Modal Analysis is the process of character- 
izing the dynamic properties of an elastic 


* 

structure by identifying its modes of vib- 
ration. A mode of vibration is a global 
property of an elastic structure. That is, 
a mode has a specific natural frequency and 
damping factor which can be identified from 
response data at practically any point on a 
structure, and it has a characteristic 
“mode shape” which identifies the mode spa- 
tially over the entire structure. 

Nodal testing is performed on mechanical 
structures in an effort to learn more about 
their elastic behavior. Once the dynamic 
properties of a structure are known its 
behavior can be predicted, and therefore 
controlled or corrected. 

Resonant frequencies, damping factors and mode 
shape data can be used directly by a mechan- 
ical designer to pinpoint weak spots in a 
structure design, or this data can also be 
used to confirm or synthesize equations of 
motion for the elastic structure. These 
differential equations can be used to simu- 
late structural response to known input 
forces and to examine the effects of pertur- 
bations in the distributed mass, stiffness, 
and damping properties of the structure in 
more detail . 

In this paper the measurement of transfer 
functions in digital form, and then the 
application of digital parameter identifi- 
cation techniques to identify modal para- 
meters from the measurer transfer function 
data are discussed. It is first shown that 
the transfer matr. , vnich is a complete 
dynamic model of an rlastic structure can be 
written in terms of the structure's modes of 


f" . . • •,«- 
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vibration. This special mathematical form 
allows one tc identify the complete dynamics 
of the structure from a much reduced set of 
test data, and is the essence of the modal 
approach to identifying the dynamics of a 
structure. 

Next, various considerations which apply to 
broad band testing using a Fourier Analyzer 
are discussed. Ways in which spectral re- 
sol utiont measurement noise, and multipoint 
excitation can be handled with a digital 
machine are also discussed. 

Finally, the application of transfer function 
models and Identification techniques for ob- 
taining modal parameters from the transfer 
function data are discussed. Several alter- 
native single degree- of- freedom and multi- 
degree-of-freedom methods are covered. 

The results presented in this paper are not 
Intended to be a survey of digital modal 
testing, but rather a sunnary of ex- 
perience gained at Hewlett Packard in the 
development of an all-digital modal an- 
alysis test system. 

II. THE STRUCTURE DYNAMIC MODEL 


In the majority of present day modal analysis 
practice, the motion of the physical system 
Is assuned to be adequately described by a 
set of simultaneous second-order linear dif- 
ferential equations of the form 

Mx{t) + Cx(t) +Kx(t) = f(t) (1) 

where 

f(t) - applied force vector 
x(t) * resulting displacement vector 
x(t) = resulting velocity vector 
x(t) = resulting acceleration vector 

and M, C, and K are called the mass, damping, 
and stiffness matrices. If the system has 
n-dtmenslons (n - degrees-of-freedom) then 
the above vectors are n-dimensional and the 
matrices are (ny-i). 

For a two dimensional system, two differential 
equations of motion would be written: 

C22jl_X2(t)J 

•'ll •'12 ^ (2) 

kg, kgg X2(t)J 

or performing the indicated matrix multipli- 
cation: 

m^^ *x^(t) + m ^2 ^2^^^ ^ ^11 ^i^^^ ^ 

Ct |2 X 2 (t) + x^(t) + k ^2 X 2 (t) * f^(t) 
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n^^x^(t) ♦ m22X2(t) + C2|X^(t) + C22X2(t) 
1^21 Xi(t) + k22X2(^) “ f2(^) (3) 


Each equation is a statement of Newton's second 
law which balances Inertia, damping and re- 
storing forces against each of the applied In- 
put forces. The mass, stiffness and damping 
matrices contain the necessary mass, stiffness 
and damping coefficients so that the equations 
of motion y^eld the correct responses when 
input forces are applied. 

The mass, stiffness and damning matrices are 
usually assumed to be real valued and sym- 
metric. 

A considerable amount of modal analysis 
practice has developed around the so called 
normal mode testing. In the theory of normal 
modes. It is assumed that the damping matrix 
Is either proportional to the mass and stiff- 
ness matrices (proportional damping), or can 
be replaced by a complex valued stiffness 
matrix (structural danv)1ng), or doesri't 
exist at all. It has been shown recently, 
(refs. 1 and 2) that without these assumptions 
a meaningful transforufation of equation (1) 
to modal coordinates is possible using con- 
jugate pairs of complex modal vectors. Fur- 
thermore, a generalized modal transformation 
of equation (1) has recently been developed 
(ref. (3)) which does not require the mass, 
stiffness and damping matrices to be real 
valued or syninetrfc. 


In the approach to modal testing presented 
here, modal parameters are Identified from 
transfer function measurements. Therefore, 
it is convenient to write equations (1) in 
their equivalent transfer function form. 

Taking the Laplace transform of the system 
equations gives: 

B(s)X(s) = F(s) (4) 

where 

F(s) = Laplace transform of applied force 
vector 

X(s) = Laplace transform of resulting 
displacement vector 

B(s) = Ms 2 + Cs + K 

s = Laplace variable (a complex number) 

B(s) is referred to as the system matrix . 

The transfer matrix H(s) is defined as the 
inverse of the system matrix, that is: 

H(s) => B{s)*l (5) 

Hence, H(s) satisfies the following equation: 

X(s) - H(s) F(s) (6) 

which is equivalent tc equation (4). 
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Each element of the transfer matrix is a 
transfer function. For an n-dimensional 
system. H(s) is an nxn matrix, which can be 
written: 


H(s) 


h|^(s) . • • h^^(s) 


h„i(s) . . . h^^(s) 


(7) 


where h.j(s) is the transfer function in the 

i— row and j— column of the transfer matrix. 
Note that transfer functions are defined in 
terms of the Laplace variable s. The transfer 
function is a complex valued function, i.e. 
it has a real and imaginary part. Figure 1 
shows the real part of a typical transfer 
function plotted as a function of the s- var- 
iable. 

Note that the Fourier transform gives the 
value of the transfer function along the 
imaginary axis in the s-plane, referred to as 
the frequency axis. The transfer function 
evaluated along the frequency axis is also 
referred to as the frequency response function . 
The real axis in the s-plane is referred to as 
the damping axis. The reason for this will be 
made clear later on. 


It will be shown in the next section that it 
is possible from measured frequency response 
functions to identify all the parameters nec- 
essary to completely specify the matrix of 
transfer functions H(s), for an elastic struc- 
ture. The transfer matrix contains all the 
information necessary to completely specify 
the dynamics of the structure and it has 
been shown that the mass, stiffness and damp- 
ing matrices can be recovered once the trans- 
fer matrix H(s) is known. 

III. MOOES OF VIBRATION 

In this section, modes of vibration are de- 
fined in mathemetical f^rms and it Is shown 
that the transfer matrix can be written in 
terms of modal frequency, damping and modal 
vectors. This special mathematical form is 
the key to modal testing using transfer fun- 
ctions and has implications both in the mea- 
surement and post test analysis portions of 
this approach. 

Recall that the elements of the system matrix 
B(s) are quadratic functions of the Laplace 
variable (s). And, since the transfer mat- 
rix H(s) is defined as the inverse of the 
system matrix, it follows that the elements 
of H(s) are ratios of polynomials In s, with 
the determinant of B(s) in each denominator. 

That is the 1j^^ element of h(s) can be 
written: 


ntAL PAnr op the laplace thanseoam 

OP A SINGLE OEOAEE-OP-PIIEEOOM 
WtTH POLES AT • - -1 tfS 



hij{s) = 


b^s^""^ +b2S^”*^ +. 


+b2n_iS +b2^_2 


det (B(s)) 


( 8 ) 


where the b‘s are polynomial coefficients. 

If the system is n- dimensional then det (B(s)) 
is always a polynomial of order 2n, which has 
2n roots, i.e. values of s for which 
det (B(s)) = 0. 


If the roots of det (B(s)) are assumed to be 
distinct , then H(s) can always be written in 
the partial fraction form: 


= Z fep, (9) 

k=l 

where 

P|^ = root of det (B(s)) = 0 

a. = residue matrix for the k^^ root. 

{nxn for an n-dimensional system) 

The roots are referred to as poles of the 

transfer function. When the system Is 
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subcritical ly damped > the polos are complex 
numbers and occur in complex conjugate pairs. 
Figure 1 shows a pair of complex conjugate 
poles locations at s = -1 + j5 Hertz. When 
the structure is critically or supercriticaMy 
damped , the poles are real valued and lie 
along the real axis in the s-plane. 

At the locations (s*Pj^) it is clear from ex- 
pression (9) that the^transfer function value 
goes to infinity, i.e. its denominator goes 
to zero. 

In physical terus, this is where the struct- 
ure is said to resonate or become resonant. 
Refer to Figure 1 again. For the example 
plotted there, when the structure is subjected 
to an input frequency of 5 Hertz, it will 
resonate because a pole is located at 
s = -1 + j5 Hertz. And the lighter the 
damping is, i.e. the closer the pole is to 
the frequency axis, the more the structure 
will resonate after external forces are re- 
moved. In the extreme case where the poles 
are located right on the frequency axis the 
structure would resonate forever when ex- 
cited by an external force. 

Each complex conjugate pair of poles cor- 
responds to a mode of vibr ation in the 
structure. They are complex numbers, 
written as 


Pfc " •'’k * ’“k ’ Pk ” -°k • ’“k 

where • denotes the conjugate, is the 

modal uampinq coefficient , (assumed to 
be positive in value) and is the 

natural frequency . These parameters, which 
are the coordinates of the poles in the 
s- plane are shown in Figure 2. 

An alternative set of coordinates can also 
be defined for describing the pole locations. 
The resonant frequency is given by: 

<"> 


and the damping factor or percent of 
critical damping is given by: 



When = l,mode (k) is critically damped; 
and when <1, mode (k) is subcritically 



Modal vectors (U|^) are defined as solutions 
to the homogeneous equation 

B(P|^) u^ = 0 (13) 

Uk = n-dimensional complex vector. 

It has been shown (ref. 2) that when they are 
defined in this way, the transfer matrix can 
be written in the form: 


H(s) 


2n 


k=l 


^“k^ 

s-P. 


.14) 


where the superscript t denotes the trans- 
pose and A^ is a scaling constant. In other 

words, the modal residue matrix a^ for 

mode (k) can be written in terms of the 
modal vector as 

«k = Vk“k^ (’5) 


Furthermore, it is shown in ref. (2) that the 
transfer matrix can be written as a summation 
of (n) conjugate pairs: 



where u^^* is the conjugate of u^ and p^* is 
the conjugate of p^. (A^ = i without loss 


damped. 


of generality.) Each of the terms In the 
summation in expression (16) is an nxn mat- 
rix which represents the contribution of 
each mode (k) to the transfer matrix. 


For a 2-dimensional system, there can only 
be two modes of vibration. Therefore, there 
are two sets of complex conjugate pairs of 
modal vectors and poles, i.e. 
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p,. u, = 
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“21. 


.“21. 
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“12 
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“22 
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and the transfer matrix is written as the 
sum of submatrices 

H{s) = 



-• 


— 


“ii“ii 

^1^21 


“ii*“ii* 

s-Pi 

s-Pi 


s-Pl* 

S-P,* 

“2l“ll 

“2l“21 

+ 

“ 21 *“ 11 * 

u .u * 
21 21 

s-Pi 

s-P,_ 


s-p,* 

s-p,* 



+ 

_ 

d 

“12"12 

“l2“22 


“12*“12* 

u, *11 * 

12 22 

S-P 2 

S-P2 


S-Pj* 

S-P2* 

“22“l2 

“22«22 


“22*“l2* 

“22*“22* 

S-P2 

• 

S-P2 

mm 


S-P2* 

S-P2* 


It is Important to note that each row and 
column of the numerators in the above matrices 
contains the same modal vector, multiplied by 
a component of itself. In other words, the 
modal vectors and pole locations can be ident- 
ified from any^ row or column of the trans- 
fer matrix except those corresponding to 
components of a modal vector which are zero, 
known as node points . This leads to perhaps 
the most significant premise of modal test- 
ing using this approach. 

Only one row or column of the transfer matrix 
needs to be measured in order to identify all 
the modal parameters of a structure provided 
the following assumptions are met 

1. The motion is linear, described by the 
linear second - order equations. 

2. The symmetry of motion or reciprocity 


property is valid (i.e.* B and H matrices 
are symmetric). 

3. No more than one mode exists at each 
pole location of the system transfer 
matrix. 


In practice, the pole locations Pj^ and one 
row or column of the residue matrix are 

actually identified from one row or column 
of measured transfer function data. Sup- 
pose that the row or column of the tran- 

fer matrix was measured. Then the q^^ row 
or column of aj^, call it r^j^, would be ident- 
ified from the data. It is straightforward 
to show that the entire residue matrix can 

be constructed from the identified q^^ row 
or column using the formula 


a = 

^ ''qqk 


(19) 


where r^^j^ is the q^^ component of the 

residue vector. If a single exciter is used 

to measure the q^^ column of the transfer 
matrix, then is the residue at the 

driving point . 


COMPLEX MODE SHAPE 

The structural model used in the foregoing 
development required that the damping matrix 
be symmetric and real valued. Without further 
restrictions on damping, such as the pro- 
portional damping assumption, the mode vectors 
can in general be complex valued. 

When the mode vectors are real valued, then 
they are the equivalent of the mode shape , 
i.e. the obrerved spatial motion of the mode. 
In the case of complex mode vectors, the 
interpretation of the mode shape is slightly 
different. Recall that the transfer matrix 
for a single mode (k) can be written 

H|.(s) » (20) 


where = (nxn) complex residue matrix, 
pj^ = pole location of mode (k). 
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A component of H|^{s) 
h(s) * — 5 

21(s-p,^) 


can therefore be written 


* 



2i(s.pJ) 


( 21 ) 


• L 

where ^ = complex residue of mode (k). 

The Inverse Laplace transform of the trans- 
fer function of expression (21) Is the 
Impulse response of mode (k). In other 
words. If only mode (k) of the structure 
could be excited with a unit Impulse 
function. Its time domain response would be 


x,^{t) = Sin{u>,^t + a,^) (22) 


where lrj^| = magnitude of the residue 
aj^ = phase angle of the residue 

= arctan In’aglnary part (r ,^) 

Real oart (r^) 

Note: It is convenient to remove a factor 

of 21 from the residue ir. expression (21) 
in order to simplify the above result. 

Expression (22) is a sinusoidal waveform with 
frequency multiplied by the decaying ex- 
ponential function as shown in Figure 3. 

Hence, the response is a decaying sinusoid 
where the amount of decay (damping) is con- 
trolled by the parameter wj^. A phase shift 

in the impulse response is introduced by 
the phase angle of the complex residue 

(modal vectoO* Note that when the residue 
vector is real valued = 0®) there 

is no phase delay in the impulse response of 
mode (k). 



Figure 3. Impulse Response of a Mode (k) 

meters necessary to characterize the modes 
of vibration of the structure. 

In this section, various methods of mea- 
suring transfer functions with a digital 
machine are discussed. In addition, digital 
techniques for reducing measurement noise and 
Increasing frequency resolution are discussed. 


1. TWO CHANNEL MEASUREMENT WITH BROADBAND 
EXCITATION 

Sine wave excitation can be used to measure 
transfer functions using a Fourier Analyzer, 
but it has some clear disadvantages cojnpared 
to other types of signals. Sine testing is 
slow since it only excites the structure at 
one frequency at a time, while the digital 
FFT is capable of generating all frequency 
lines in a band of interest simultaneously. 
Therefore, a broadband excitation signal such 
as an impulse or random noise source which 
excites .the structure over a spectrum of 
frequencies simultaneously makes transfer 
function measurement faster with a Fourier 
Analyzer. 

Secondly, since modal analysis assumes that 
the dymanic behavior of the structure is 
linear, any measured non-linear behavior will 
tend to invalidate the modal parameter est- 
imates. A broadband random signal is act- 
ually helpful in removing the effects of non- 
linear behavior from the measured transfer 
function data. 


It is this phase delay which is represented 
by the complex mode shape. 


III. MODAL TESTING CONSIDERATIONS 


In the last section, 
one row of colunn of 
be measured in order 


it was shown that only 
the transfer matrix need 
to identify the para- 


TRANSIENT TESTING 

A fast method of performing transient test- 
ing is to use a hand held hammer with a load 
cell attached to impact the structure. The 
load cell measures the i.iput force, and 
accelerometers mounted on the structure mea- 
sure the response signals. The setup time 
required to use this method is negligible 
compared to the mounting and alignment of 
shakers. A typical te?‘ setup is shown in 
Figure 4. 
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Whether a hammer or shaker is used to Impact 
the structure, the signal measuring equip- 
ment must contain a triggering mechanism 
capable of capturing the majority of this 
short duration signal. 

The pi iTidry drawback of the transient test- 
ing method is that the power spectrum of the 
Impulse cannot be controlled as well as a 
random spectrum. Figure 5. A shows a typical 
Impulse signal and its corresponding power 
spectrum. This particular signal prov^* ^es 
energy from DC (0 frequency) to 5 KHz, but 
the energy density may not be high enough 
to excite an entire structure. The ^'nly 
way of providing more energy is to hit the 
structure harder, whict may destroy It! 




Figure SB. Impulse Signal and Power Spectrum 

On the other hand, a softer head can be 
placed on the hanmer or shaker head to soften 
the impact. This has the effect of changing 
the impulse to a wider time width pulse, and 
the bandwidth of the corresponding power 
spectrum will in general be less than that 
of the Impulse. Figure 5.B illustrates this 
case. This method can be used to concentrate 
more energy at lower frequencies but still 
may not provide the energy levels necessary 
to test a large, heavily damped structure. 

RANDOM TESTING 

Two primary advantages of using a shaker 
driven by a broadband random signal are 
that it is more controllable than a tran- 
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slent signal and when used In conjunction 
with power spectrum ensemble averaging 
(described later), noise and non*1 Inear 
components of the structural dynamics can be 
removed from the measured data. 

A random signal can be generated In several 
different ways: 

METHOD #1. 

A straightforward approach Is to use a random 
signal generator and pass the signal through 
a band pass filter to concentrate the energy 
In a band of interest. Figure 6. Illustrates 
this case. G<»nerally, the signal spectrum 
will be flat except for filter roll-off at 
the ends of the spectrum and you can only 
control the overall level of the signal 
(I.e. the overall spectrum level). One draw- 
back of this approach Is that although the 
shak'Ar 1$ being driven with a flat spectrum 
signal, the structure Is excited by a signal 
with a different spectnmi due to the Impedance 
match between the structure and the shaker 
head. To overcome this, a signal with a 
shaped spectrum must be output to the shaker. 

Another more serious drawback of this method 
Is that the measured signals (Input and re- 
sponse) are not periodic 1n the measurement 
time window. A key assumption of digital 
Fourier Analysis Is that a measured time 
waveform be periodic In the measurement time 
window. If It Is not, the corresponding fre- 
quency spectrum will contain so called "leak- 
age”. That is, energy from the non-periodic 


parts of the signal will “leak" Into the 
periodic parts of the spectrum, this giving 
an Inaccurate result. Figure 7. Illustrates 
the difference between a signal which is per- 
iodic In the window and one which Is non- 
periodic. Figure 7. A. shows a sine wave of 
period T. The measurement window covers the 
first time Interval (0, T) but note that if 
the measurement were repeated again In the 
Internal (T, ZT), the exact same waveform 
would be measured. When a signal repeats In 
each successive time period (T), It Is said 
to be periodic In the window (0, T). 

Note that the non-periodic signal In Figure 
'’.B. does not repeat In each successive time 
period (T). In summary a signal x(t) Is 
periodic In a window (T) If and only If 

x(t + T) » x(t) 

Typical Input and response signal pairs for 
transient and random measurements are shown 
In Figure 8. The transient waveform In 
Figure 8. A. are periodic In the window,' 
However, the random signals are non-periodic 
In the window (0, T) since they do not re- 
peat In successive windows of length T. 

The Hanning Window 

The Hanning window Is normally used, to re- 
duce leakage In the spectrum of a non-per- 
lodlc random signal. 

The Hanning window Is shown In Figure 9. 

When a non-periodic time waveform Is multl- 
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Figure 6. Random Test Setup 
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plied by this wlridow. the values of the 
signal In the measurement window r.wre closely 
satisfy the requirements for a periodic 
signal. 



Figure 7A. Signal Ferioaie in Window 


corresponding Fourier transforms. Hence, 
although multiplication of a non-periodic 
signal by the Hanning window reduces leakage, 
the spectrum of the signal Is still distorted 
due to convolution with the Fourier trans- 
form of the Hanning window, I.e. Its line 
shape. 





Figurt 7B. Signal Non-pertocMc in tha Window 

Therefore, leakage in the spectrum of a signal 
which has been multiplied by the Hanning win- 
dow Is greatly reduced. 

However, multiplication of two time waveforms, 
I.e. the non-periodic signal and the Hanning 
window. Is equivalent to convolution of their 


Figura 9. Hanning Window 

Therefore, when a ncn-periodic random noise 
source Is used to excite the structure, a 
Hanning wir.low must be applied to the mea- 
sured time-domain signals and the effect of 
this window must be deconvolved from the 
transfer function. This Is easily done 
using digital techniques provided that the 
spectrum of the Input signal is flat. 

METHOD #2 

To avoid the problems of using a non-periodic 
random signal. It Is pout’tie with a digital 
computer to generate a random s'^gnal 
and then output It to the shaker through a 
digital to analog converter (DAC). To do 
this, enough random data Is synthesized In 
the computer to drive the shaker for T sec- 




RESPONSE 


Figure 8A. Traniient Input and Response Signets 




RESPONSE 


Figure 8B. Random Input and Response Signals 






V 




ondst and this data Is then repeatedly out- 
put through the ADC every T seconds. The 
random signal Is periodic since It repeats 
every T seconds, and Is referred to as a 
pseudo random signal. This method of ex- 
cftatlon not only removes windowing problems, 
since the measuted signals are periodic, but 
allows you to shape the spectrum of the ran- 
dom signal line by line to compensate for 
loading of the shaker by the structure. 

However, because It Is periodic, a pseudo 
random signal will not effectively remove 
nonlinear distortion components from the 
measurement. 

METHOD # 3 

The best possible random signal Is one that 
Is both periodic In the window, I.e. satis- 
fies the conditions for a periodic signal 
and yet changes with time so that It excites 
the structure In a truly random manner. 

This type of signal is referred to as a 
periodic random signal. It starts out 
exciting the structure In the same manner 
as a pseudo random signal. Then, after the 
structure has reached a steady state condi- 
tion, I.e. the transient part of Its motion 
has dlea out, a measurement Is taken. Then, 
Instead of continuing with the same pseudo 
random sr^nal, a different one Is synthesized 
and output to the structure. This new ex- 
citation will excite the structure In a dif- 
ferent steady state manner than the previous 
one, and when the power spectrums of these 
two measurements and many others are averaged 
together, non-linear distortion components 
are removed from the measured result. Ref- 
erence (4) explains these methods In more 
detail. 


2. REDUCING MEASUREMENT NOISE 

One of the problems of any modal testing 
system is that extraneous noise from various 
sources in the measurement equipment or on 
the structure Itself Is measured along with 
the desired signals. Figure 10 below depicts 
this situation for the case involving trans- 
fer function measurements. 

where 

F(w) * Fourier transform of measured input 
signal 



> vu> 

— 


N(w) 


Figtir* 10. Mtwurtmtnt of Signal piui Noim 


Y(u) » Fourier transform of measured re- 
sponse signal 

N(w) * Fourier transform of noise 

X(w) « Fourier transform of desired re- 
sponse measurement 

H(w) * Transfer function of structure 

Since we are Interested In Identifying modal 
parameters from measured transfer functions, 
the va^^lance on the parameter estimates Is 
reduced In proportion to the amount of noise 
reduction in the measurements. Two methods 
of noise reduction which take advantage of 
the digital data form and the data proces- 
sing capabilities of a Fourier Analyzer are 
discussed here. One is an averaging tech- 
nique applied at measurement time, and the 
other Is a smoothing technique which Is ap- 
plied after the measurements have been made. 

COMPUTATION OF THE TRANSFER FUNCTION IN 
THE PRESENCE OF NOISE 

In Figure 10. aoove, we wish to compute 
H{o)) from the measured and transformed sig- 
nals F(o)) and Y{w). Recall that H(w) Is 
defined as the ratio X{w)/F{iu). Dropping 
the frequency dependence from the notation. 
Figure 10. yields the expression 

Y « HF + N (23) 

All quantities in equation (23) are scalers. 
The input (or auto) power spectrum is defined 
as 


Off = FF* (24) 

and the cross-power spectrum as 

- YF* • (HF+N)F* = + NF* (25) 

where • denotes the complex conjugate of the 
transform. 

Now consider averaging the quantity 
The average value of Sf from n different 
measurement s Is 


n 

V ; 2 (!6) 

1=1 

where Gy^(1) Is the 1^^ measurement taken. 
Therefore, 


Gy^ = HG^^ + PF* 


HGff * 


(27) 
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where is assumed to be invariant (i.e. 


THE COHERENCE FUNCTION 


not changing from measurement to meabureinent) 
and and 6^^ are defined similarly tj 

Now if the noise is assumed to h.^‘/e zero mean 
value and to be incoherent with the input 
signal, then 

H = ^ ^ (28) 

®ff °ff 


As the number of averages grows larger, the 

noise term is reduced and the ratio 

more accurately estimates the true transfer 
function. 


Whenever transfer functions are measured on 
a digital Fourier Analyzer in the manner 
just discussed, the coherence function can 
also be easily computed. The coherence 
function denoted y^(w) is defined as the 
ratio 

2 = (response power caused by applied input ) 
^ " {measured response /«q- 


The coherence function is also dependent upon 
frequency. Recall from Figure 10. that the 
measured response power spectrum coiitains the 
power caused by the input and the power due 
to extraneous noise sources. The measured 
response power spectrum is written 


Gyy = (HF+N) (HF+N)* 

= |HF|^+ H*NF*+ HFN* + NN* 


(30) 


Figure 11. shows che effect of averaging on 
a typical transfer function measurement. 


or 


‘5yy ' |Hp + H + G^I^ 


(31) 
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Figure 11 A. Measi*renient with 3 Averages 


The response power caused by the 
input is 

Sx “ ^ff ^32) 

While the cross oower spectrum between mea- 
sured Input and response is 

Gyf = HGf^ (33) 

Substituting these into the definition for 
coherence and taking average values from a 
large number of measurements, as explained 
with the previous transfer function calcu- 
lation, gives 


v2 = 




(34) 
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0 frequency 1 kHi 

Figure 11B. Measurement with 50 Averages 


From the definition, it is clear that when 
the measured response power is in fact caused 
ty the measured input power, the coherence 
value is one, for all freque":ies. But when 
the measured response power is greater than 
the measured input power because some extra- 
eous noise source is contributing to the out- 
put power, then the coherence value will be 
less than one (but greater than zero) for 
those frequencies where the noise source 
adds power to the response signal. 

Hence, tne coherence function is used to in- 
dicate the degree of noise contamipation in 
a transfer function measurement and it can 
be used in a qualitative way to determine 
hew much averaging is necessary to remove 
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noise front the measurement. Figure 12. 
shows coherence functions corresponding to 
the transfer function measurement In Fig- 
ure 11. Note that with more averaging the 
estimate of coherence contains less vari- 
ance, thus giving a better estimate of the 
noise energy in the measured signal. 





Figure 12A. Cdherenoe o* Measurement with Z Averages 
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Figure 128. Coherence of Measurement with 50 Averages 

Figure 13. contains a block diagram of the 
power spectrum averaging process used in a 
digital machine to reduce noise in t»*ansfer 
function measurements. 

EXPONENTIAL SMOOTHING 

If after a reasonable amount of averaging, 
the transfer function measure»«ents are 
still relatively noisy, the noise can be 
further removed by means of an exponential 
smoothing process. Many different types of 
schemes couV be used to smooth data, 
but the exponential funccion is advantageous 
in thi , case because its effect on the data 
1 » <nown. 

Recall from section II that the inverse 
Fourier transform of a transfer function is 
the sum of the impulse responses of the modes 
in the measurement bandwidth, i.e. 


m 

x(t) = 2 !r,^|e ( 35 ) 

k=1 


An example of this sum of decaying sinusoids 
is given in Figure 14. A. Measurement noise 
adds uniformly to this impulse response and 
therefore gives tie result shown in Figure 


14, B. 



Since the noise is distributed evenly through- 
out the block and the signal decays exponent- 
ially, the signal- to-noise ratio of the com- 
bined signal also decays exponentially. 

Now if the measurement plus noise is multi- 
plier by a decreasing exponential function 
(called an exponential window), as shown in 
Figure 15., the noise at the right hand end 
of the block is truncated, while the signal 
toward the left hand end is preserved. The 
overall signal to noise ratio of the data is 
increased since the combined signal plus noise 
is weighted more heavily in favor of the 
signal and less in favor of the noise. When 
the resulting waveform is transformed bxck to 
the frequency domain, the corresponding 
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transfer function has been smoothed, as 
Illustrated in Figure 16. Furthermore, 
the width of the modal resonances, which are 
governed by the amount of damping In each 
mode have Increased by a known amount. If 
the exponential window Is represented by 

«(t) = e-*»‘ (36) 

b = a known constant, then the Impulse re- 
sponse after windowing Is 


«(t) = Z + a^) (37) 

k«l 




Figure 15. Expofientizl Window 
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Figure 14B. impulse Response plus Random Noise 

In other words, a known amount of damping (b) 
Is added to each mode with each multiplication 
of the Impulse response by the exponential 
window. 


SOO Mc 

Figure 16A. Transfer Function Before Smoothing 
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Figure 16B. Transfer Function After Smoothing 


The smoothed transfer function data can there- 
fore be used to Identify modal parameteis and 
the correct damping coefficient can be recov- 
ered by subtracting the amount of damping oue 
to the smoothing process from the identified 
damping value. 


The drawback of this approach to noise removal 
is that if modes are closely spaced in fre- 
quency, exponential smoothing will smear them 
together so that they are no longer discern- 
able as two modes. 
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3. GAINING FREQUENCY RESCLUTION 

Normally, the Fourier transform Is computed 
In a frequency range from zero frequency 
(D.C.) to some maximum frequency 

digital Fourier transform is spread over a 
fixed number of frequency lines (typically 
1024) which therefore limits the frequency re- 
solution between lines. 

Band-selectable Fourier Analysis, the so cal- 
led ZOOM transform, is a measurement tech- 
nique in which Fourier transform based dig- 
ital spectrum analysis is performed over a 
frequency band whose upper and lower fre- 
4 uencies are independently selectable. A 
comparison of Band- selectable Fourier Ana- 
lysis (BSFA) and standard (Baseband Fourier 
Analysis is shown in Figure 17. 

BSFA can provide an Improvement In frequency 
resolution of more than a factor of 100, as 
well as a 10 dB Increase in dynamic range, 
compared to Baseband Fourier Analysis. It 
Is also called the ZOOM transform since it 
zooms In on a portion of ^ Baseband spectrun 
and "magnifies** it with more lines of defini- 
tion, much as a camera zoom lens magnifies a 
picture. The topics of resolution and dy- 
namic range are discussed separately below. 
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Figure 17. Bum Band vs. Zoom Transform 


INCREASED RESOLUTION 

In any measurement technique, the resolution 
achievable in the frequency domain is deter- 
mined by the length of time that the time- 
domain signal is observed. Specifically, the 
frequency domain resolution is the reciprocal 
of the time length of the measurement 
(4f = l/T). 

Baseband Fourier Analysis provides uniform 
frequency resolution from D.C. to F^^ (one- 

naif the sampling frequency). Thus, the 
frequency resolution can also be expressed as 

af = F ^ /(N/2), where N is the block size; 
max 

the number of samples describing the real 
time function. There are N/2 complex (mag- 
nitude and phase) samples In the frequency 
domain. 

In actual practice, F Is fixed by the 

frequencies of major experimental Interest, 
and by aliasing considerations. Thus, the 
only way to improve frequency resolution 
in Baseband Fourier Analysis is to Increase 
the block size. There are two reasons why 
this Is an Inefficient way to increase fre- 
quency resolution: 

A. Digital processing times Increase with 
block size. 

B. The maximum block size is limited to 
some -datively small number of sam- 
ples, based on Computer memory size. 

BSFA solves these problems by providing 
greatly Increased resolution about points 
of Interest in the frequency domain, with- 
out requiring an Increase In the system's 
block size. 

This is done by digitally filtering the 
Incoming time-domain data, and storing only 
the filtered time-domain data, corresponding 
to the frequency domain band of Interest. 

Since the frequency resolution is still 
the reci,>rocal of the time length of the 
incoming signal, the digital filters must 
process T seconds of data to obtain a fre- 
quency resolution in the analysis band of 
Af * 1/T. The resolution obtained in the 
frequency band of interest is approximately 

if = bandwidth/(N/2) (38) 

Thus, by restricting attention to a narrow 
region of interest below an increase 

in frequency resolution proportional to 

(where BW is the BSFA measurement 
max 

bandwidth) can be obtained. 

increased dynamic range 

BSFA can provide increased dynamic range re- 
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lative to Baseband Fourier Analysis* This 
is due to the increased ** processor gain” on 
the analog to digital converter (AOC) quant- 
ization noise. However, a BSFA system can 
take advantage of this processor gain only 
to the extent that is taade possible by the 
noise level and out-of-band rejection of 
the pre-processing digital filters. Certain 
BSFA filters can provide more than 90 dB of 
signal-to-noise and out-of-band rejection 
(Ref 5.). 

Processor gain refers to the effect of in- 
creased frequency resolution on white noise 
in the presence of a narrow-band signal, 
such as a sine wave. The sine wave energy 
exists at a single frequency. Its peak value 
is, therefore, independent of the frequency 
resolution. The white noise peak amplitude, 
however, is reduced 3 dB in power each 
time frequency resolution is increased by a 
factor of 2. BSFA provides increased pro- 
cessor gain by increasing frequency re- 
solution in the analysis band (relative to 
the baseband measurement) . 

4. SINGLE POINT EXCITATION 

Having discussed how single transfer functions 
are measured, this section returns to the 
fundamental modal testing problem, measurement 
of one row ot column of the transfer matrix. 
Consider the following two-dimensional ex- 
ample which defines the transfer matrix as: 


m m 

x,(s) 

h 

X2(s) 

h 


h'ii(s) h|«(s) 


12' 

22< 



F,(S) 

1 

F2(s) 


(39) 


In order to measure the first row of transfer 
functions along the frequency axis in the 
s-plane, that is to measure the frequency 
response functions 


hii(w) “ h<ji 


(s)|, h^ 2 (u,) = h^ 2 (s) (40) 


s^iti 


S=lu» 


the structure should first be excited at 
point #1 and its response simultaneously 
measured at point #1. Then the first trans- 
fer function h,,(u)) is computed by forming 
the ratio of 

Xi(w) ^ 

, F2(«) = 0 «1) 

F, (u)) = Fourier transform of the Input at 
point <1 

X,(u>) = Fourier transform of the response 
at point #1 


Next, the structure would be excited at point 
#2 and the response at point #1 simultaneously 
measured. Then the transfer function 

is computed by forming the ratio 


X|(w) 

**12^“^ * F'^Cw) • = 0 (42) 

F^(w) = Fourier transform of the input 
at point #2 


Note that the first row of the transfer 
matrix is the same as the first column 
since H(s) is synmetric, i.e. ’ 


To measure one column of the transfer ir.atrix, 
the exciter is placed at one point on the 
structure (point #1 to measure column II; 
point 12 to measure column 12) and the re- 
sponses at points #1 and #2 are measured. 

For Instance, column #2 would be measured 
by forming the ratios: 

• This measurement Is actually made using 
expression (28) but for simplicity it is 
represented here as the ratio of transformed 
response over input. 

as defined before, and 

X^(oi) 

‘'22<“> ' f|q * ° 

where X 2 (w) = Fourier transform of response 
at point #2. 

The method just described requires that the 
structure be excited at only one point at a 
time, i.e. single point excitation . 

5. MULTI-POINT EXCITATION 

In many cases a structure may be so large 
or heavily damped that a single point ex- 
citer does not put enough energy into the 
structure to excite the modes of interest. 

In addition, the structure may have multi- 
ple modes at the same frequency. In these 
cases, two or more exciters may be required 
to excite the structure. All of these 
Inputs and the corresponding responses must 
of course be measured and because the data 
Is in digital form, the desired row or 
column of the transfer matrix can still be 


formed using matrix algebra. 

In general, the transfer matrix is defined by 

x(s) = H(s) F(s) («) 

Post multiplying by the conjugate transpose 
(denoted by T) of the transformed input 
vector gives 

X(s) F^{s) = H(s) F{s) F(S)^ (46) 

then solving expression v46) for H(s) gives 

H(s) = X(s) F^(S) [F(s) F(s)^r’ (47) 

If iiH inputs are used to excite an n-dimensional 
system (m < n), then the corresponding m coi- 
unns of H(s) can be computed using expression 
(47). Note that the inputs cannot be com- 
pletely correlated or the indicated inverse 
will not exist. 

Consider the following case where two inputs, 
located at points i and j on the structure 
arc used 



'x,(s)' 


• 

- 


★ * 


★ 

H(s) = 


[F^.F.] 

[F.,F.] 

i’ j** 


ys)_ 


• 



Dropping the s-variable dependence from 
the notation in expression (48) gives 


H(s) = 


X,F* 

X,F* 




w V* 

w V* 



F.F* 

Vi 



J jJ 


-1 


(49) 


= QP 


,-l 


Q is an n X 2 matrix and P is a 2 x 2 matrix. 
The inverse of P can written in closed form 
using Cramar's rule, aiid the indicated mult* 
iplication of expression (ag) can be carried 
out using the Fourier transforms of input and 
response signals to yield frequency response 

functions corresponding to the i^^ and 
columns of the transfer matrix. 

This procedure can be extended to larger 
numbers of inputs but the required algebra 
involving the input and response Fourier 
transforms becomes more complex. 


IV. MODAL PARAMETER IDENTIFICATION 

When a structure is excited by a broad band 
input force many of its modes of vibration 
(degrees-of-freedoml are excited simultaneously. 
Since the structure is assumed to behave in a 
linear manner its transfer functions are 
really the suns of the resonance curves for 
each of its modes of vibration, as shown in 
Figure 18. 

In other words, at any given frequency the 
transfer function represents the sum of 
motion of all the modes of vibration which 
have been excited. However, near the nat- 
ural frequency of a particular mode its con- 
tribution to the overall motion is generally 
the greatest. The degree of mode overlap , 
i.e. the contribution of the tails of ad- 
jacent modal resonance curves to the transfer 
function magnitude at a modes' natural fre- 
quency, is governed by the amount of damping 
of the modes and their frequency separation. 
Figure 19 illustrates light and heavy modal 
overlap. Figure 19. A. shows modes with 
light damping and sufficient separation 
so that there is little modal overlap. 

Figure 19. B. shows modes with heavy damping 
and/or sufficiently high modal density 



Figure 18. Magnitude of a Multi Degree-of* Freedom System Transfer Function 
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(little frequency separation) such that there 
is plenty of modal overlap. This condition 
is also referred to as modal coupling . 


This complex resonance curve generated by 
expression (50) along the frequency axis 
(s=iu) is shown in Figure 20. 


In cases where modal overlap is light, the 
transfer function data can be considered in 
the vicinity of each modal resonance as if 
it were a single deqree-of- freedom system. 

In other words it is asstjned that the con- 
tribution of the tails of adjacent modes near 
each modal resonince i> negligibly small. 

On the other hard when modal overlap is 
heavy a single aegree-of- freedom approach to 
modal parameter identification will not work; 
the parameters of all the modes must be id- 
entified simultaneously. 

1. Single Degree-of-Freedom Hodels 

The nk ity of modal parameter identifi- 
cation techniques used today are based upon 
single degree-of- freedom models. This is 
so because the presently popular normal mode 
multi -shaker testing techniques re pre- 
dicated upon the excitation of one mode of 
vibration at a time and hence the measured 
transfer function primarily reflects the 
motion of a single mode of vibration. 



Figure 20. Resonance Curve for Singte-Degree-of-Freedom 

The resonance curve is generated by summing 
two terms together, one that describes the 
motion primarily around the positive fre- 
quency u)|^, and one that describes the motioi. 


Recall from expression (21) that the transfer 
function of a single mode (k) of a linear 
system can be represented as 




* 



21(s-P|^) 2i(s-pj^*) 


(50) 


primarily around the negative frequency 


Normally the transfer function is only mea- 
sured for positive frequencies, i.e. along the 
+iw axis, and it is clear that for positive 
frequencies the majority of the modal reson- 
ance is described by the formula 


where 

p|^ * pole location in the s-plane 


H,(S) = 


2i(s-p,j) 


(51) 


- complex residue 


The frequency response function Hj^(w) can be 


f - 






( 
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written In real and imaginary parts as 


Re[H^(..)] = I 

where 

"k = ■'’k ^ 

"k " -'^k * ’“k 



s = iw 


(52) 


Now at modal resonance, i.e., w = wj^, it Is 
clear from the expressions above that 


where 

|rj^| = magnitude of complex residue 
a|^ = phase angle of residue 
= arctan ( Im(r|^)/Re( r^^) ) 

If the residue is real valued and of unit 
magnitude then 



It is straightforward to show that 


Re 

Im [H|^(uj|^) ] s- rl|^/2aj^ 


(53} 


MODE SHAPE FROM QUADRATURE RESPONSE 

Ti.'' real part of H|^(cj) is also referred to 

as the coincident or in-phase response, and 
the imaginary part of H|^(w) is called the 

quadrature or out-of-phase response. The 
so 'called co-quad plot is simply the trans- 
fer function plotted in rectangular (real- 
i mag inary) coordinates along the frequency 
axis. 


In cases of light damping the residues are 
found to be almost always real valued. There- 
fore r2=0 and Re(Hj^((j))=0. When this assump- 
tion is valid the mode shape can be identified 
by merely picking the quadrature response of 
each measurement at the modal resonant fre- 
quency. Figure 21 illustrates this case. 

MODE SHAPE BY CIRCLE FITTING 

Another single degree-of- freedom approach to 
identifying the mode shape is based upon use 
of the Nyquist (or real vs. imaginary) plot 
of the transfer function. Expression (51) 
for a single mode (k) can be rewritten 



Re[H,]2 . (In[H,] * (56) 

k k 

That is, the modal resonance curve .s a circle 
in the Nyquist plane with radius 1/4 and 

centered at -i/4 as shown in Figure 22. 

The complex residue vector merely 



MOOf SHA« - OSTAINtD AT SAWI FREOUtNCV FROM Atl 

measurement points 


Figure 21. Mode Shape from Quadrature Re$pon$e 
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expands the radius of the circle and rotates 
the circle In the Nyquist plane by the clock- 
wise angle away from the negative Imag- 
inary axis, as shown In Figure 23. 


MODE SHAPE BY COMPLEX DIVISION 

This technique Is also based upon the reson- 
ance curve for the positive frequency pole 
of a single mode, which Is described bv 



REAL 


IMAG. 



In genera] circle fitting involves the fitting 
of a circle in parametric form 

(x + a)^ + (y + b)2 = (57) 

by a least squared error procedure to a nuni- 
ber of measured transfer function data points 
in the vicinity of each modal resonant fre- 
quency. In the process the center of the 
circle (a, b,) and its radius (c) are ident- 
ified. The radius of the circle is then used 
as the magnitude of a component of the mode 
shape and the phase angle is computed as 

= arctari (j) (58) 

The complete mode shape is obtained by fitting 
ci»^cles to each measurement of an entire row 
or column of the transfer matrix. This method, 
which is referred to as the Kennedy-Pancu 
method (Ref. 6), is generally more accurate 
than simply using the quadrature response. 




(59) 


Since p, the pole location. Is assumed to be 
constant no matter where a transfer function 
measurement Is made on the structure, then 
the ratio of two transfer functions (n) and 
(m) is 



(60) 


In other words dividing one transfer function 
by another yields a complex scalor at all 
frequencies . Therefore If a single transfer 
function (m) is divided into all other mea- 
surements from a row or column of the trans- 
fer matrix, and several values of the re- 
sultant quotients in the vicinity of each 
modal resonance are averaged together after 
each divide operation, the complex vector 




n m 


(61) 


is identified, which is proportional to the 
mode shape. 

This method is also better than the quadrature 
response method since, like the circle fitting 
method, it uses a number of data points in the 
vicinity of the modal resonance. However, it 
is sensitive to shifts in the modal fre- 
quency from measurement to measurement. When 
these shifts do occur, the resulting mode 
shape has small phase changes in it. 

MODE SHAPE BY A DIFFERENCE FORMULA 

This formula is based upon the differencing 
of transfer function data between adjacent 
frequency lines. In fact both residues and 
pole locations can be determined by this dif- 
ferencing approach. 

Let the transfer function for a single mode 
(k) at frequency uj be defined by 



i 
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where ' wj, 

and at the previous frequency by 



The product of these two values Is 


"A-.-l > S^HV. * <4 

and the backward difference {aHj = 
is 


■» 4 

^"j=2 ^ 1c,) 

u 


(65) 


Dividing the difference into the product gives 
the desired residue value 



Hence by forming the product and difference 
of adjacent frequency lines of data the 
residue is estimated using the above formula, 
which also includes the frequency resolution 
between lines (w.- This formula ap- 

j j” * 

plies at all frequencies and, as in the 
previous case, is best applied to data in 
the vicinity of the modal resonance and then 
an average of the results computed. 

POLE LOCATION A DIFFERENCE FORMULA 

Suppose that expression (62) is first multi- 
plied by ojj to give 



(67) 


and then the difference formed 




[ S-“i-l)^"k ^ ''kl 


Then the desired parameters and are 
obtained by forming the ratio 


"k + i-k = 



(69) 


and this formula also applies at each fre- 
quency Wj. 

Both of the above differencing formulas can 
be derived using higher order differences 
also. 


2. MULTIPLE DEGREE-OF-FREEDOM MODELS 

All of the above methods work reasonable well 
when the amount of overlap between modes is 
small. However, in cases where the amount 
of modal overlap is sufficient to cause 
significant errors by the single degree-of- 
freedom techniques a multi-degree-of- freedom 
identification technique must be used. 

In general some type of curve fitting tech- 
nique which matches the following summation 
expression evaluated along the frequency 
axis (s=jo)) to measured frequency response 
function data is necessary. 



For each mode (k) the four parameters fre- 
quency, damping and complex residue (two 
parameters) must be estimated simultaneously. 

LEAST SQUARES ESTIMATION 

Oi.e technique that has been successfully ap- 
plied at Hewlett-Packard is the application 
of equation (70) by an iterative least squared 
error technique to the complex transfer 
function data. The iterative equations re- 
sult from a minimization of the error func- 
tion 


N 

e = («j - (71) 

j=l 

where 

H. = the measured data at frequency u> ■ 

J 

H((j. ) - the analytical model a' frequency 
N = data block size 


( 68 ) 



with each Iteration of the algorithm new 
values of all four unknowns for each mode 
are estimated which further reduce the error 
between the measured and analytical waveforms. 
Figure 24, Illustrates a mul t1- degree- of- 
freedom curve fit using this process. 
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beforehand and Initial estimates of their pole 
locations given as starting values. Modes 
are generally Identified from their re- 
sonance peaks on the transfer function mag- 
nitude curve. This can be done either by 
a computerized peak picking procedure, or 
by hand. Starting estimates of modal damp- 
ing and residue can usually be satisfact- 
orily Identified using one of the previously 
described single deg ree-of- freedom methods 
or some other technique. 

This apparent disadvantage, I.e. having to 
Identify the modes before the least squared 
fitting process begins, turns out to be an 
advantage from an operational standpoint. 

In a test situation the machine operator Is 
able to maintain control over the modal par- 
ameter Identification process by specifying 
modes at certain frequencies even though 
from particular measurements their presense 
may not be apparent. 



Figure 24. Magnitude Curve - Multi Degree<of-Freedom Curve Fit 


An Alternate Approach 


From inspection of expression (70), it 
Is clear that the transfer function Is a 
linear function of the residues 

Therefore, If the pole locations were 

known, the residues could be Identified 

using a least squared error procedure that 
involves solving a set of simultaneous lin- 
ear equations. Hence, a simultaneous lin- 
ear equation solution algorithm coniilned 
with a search procedure which iterates to- 
ward an optimum estimate of pj^ Is another 

approach to Identifying modal parameters 
with a m»ilt1- degree- of- freedom model. 

With either of the Iterative techniques just 
mentioned the modes must somehow be Identified 


THE COMPLEX EXPONENTIAL ALGORITHM 

A direct solution of the modal parameter 
Identification problem underwent a consider- 
able amount of development for the U.S. Navy 
several years ago (Ref. 7). It 1$ refer- 
red to as the complex exponential algorithm. 
The approach Is direct In the sense that no 
starting modal estimates are necessary. The 
algorithm works with the Inverse transform 
of the transfer function i.e., the Impulse 
response. Expression (22) the characteristic 
damping sinusoidal response, can also be 
written using complex exponential functions 
as 

k-1 

Using the complex exponential algorithm 
this expression is curve fit to one mea- 
sured row or column of the Impulse responses 
(corresponding to a row or column of the 
transfer matrix) and the parameters pj^ and 

are Identified in each measurement. 

The algorithm Is a two step process which 
solves fo** all the polynomial roots Pj^ first 

(recall that the poles are roots of the 

polynomial c.et (B(s) = 0) and then the re- 
sidues r|^. Originally the algorithm would 

find (n/4) modes in a impulse response that 
Is (n) data points long. However, it has 
been Impv'oved recently so that parameters 
for a prespecified number of modes are Id- 
entified. Nevertheless, the weakness of 
this approach is that the machine operator 
has no control over the root solving pro- 
cess and hence the algorithm will put modes 
wherever necessary to fit the data. There- 
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fore the same mode niiy not be Identified 
In different transfer function measurements, 
thus making It difficult to Identify a mode 
shape by this process. 

V. CONCLUSIONS 

The development of al1-d1g1tal systems to 
collect and process dynamic vibration data' 
from mechanical structures Is beginning to 
mature more quickly as test engineers real- 
ize the power of the digital Fast Fourier 
transform as a tool for transfer function 
measurement. The testing and analysis tech- 
niques presented here are but some of the 
approaches that can be used with measure- 
ments that are in digital form. 

Although analog based single mode testing 
techniques are currently more widely used 
than the digital techniques presented here, 
there are some Inherent advantages to broad 
band testing e.g., speed, accuracy, resol- 
ution, noise rejection, and repeatability, 
which indicate that more Fourier transform 
based systems will be used for modal test- 
ing and analysis in the future. 
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HARDWARE/SOFTWARE DEVELOPMENT TECHNIQUES TO SUPPORT 
COMPUTER CONTROLLED TEST SYSTEMS 


Laurie R. Burrow, Jr. 

Spectral Dynamics Corporation of San Die^o 
San Diego, California 


An economical systems approach is described for controlling a desired 
random spectrum shape at a reference point. Hardware concepts for performing 
a "Power of 2" inverse FFl’ from a 500-point real time spectrum and software 
considerations for control and simple user dialog are presented . Examples are 
given of a complete test sequence including system calibration, establishing 
safety criteria, and displaying calibrated results. The extension of this concept 
to shock analysis and control is also given. 


INTRODUCTION 

The process of ucing the Fast Fourier 
Transform (FFT) algorithm coupled to a computer 
for generating cu«d controlling a random vibration 
test has evolved into a practical and realistic 
concept today . 

The basic approach is that of measuring 
a response from a control accelerometer using 
a forward transform, then comparing this with 
an operator-entered reference spectrum. Next, 
creating a spectrum that when converted to a 
random time function by use of the inverse FFT 
and used to excite the shaker system , would 
result in a response spectrum that matched the 
reference spectrum within prescri^ ed limits. 

An investigation into the available systems indi- 
cated that a less complex and less costly approach 
was needed to really make digital random control 
and its inherent benefits available on a practical 
basis . 

The purpose of this paper is to describe 
the hardware and software design considera- 
tions for the Spectral Dynamics Digital Random 
Control System . Many system details not con- 
sidered pertinent are omitted . As in most cases , 
we were somewhat biased in our approach due 
mainly to our years of experience in the environ- 
mental field, as well as our practice of making 
minimum use of other manufacturers' hardware 
in SDC systems . 


spectral DYNAMICS* APPROACH 

O nginally , it was intended that a version 
of the SD Digital Signal Processor (DSP) , 

Fig. 1, would be used to perform forward and 
inverse Fourier Transforms because of its high 
processing speeds . The basic cost of FFT hard- 
ware caused a re-evaluation of plans for a Digital 
Random Control System . After a review of the 
analyzers available , a high speed 500-line Real 
Time Spectrum Analyzer, Fig. 2, was chosen 
as the Power Spect* al Density Analyzer since 
it had the required Malysis ?*peed, economy, 
and reliability. Having made this decision, 
a basic approach was formulated which needed 
to be unique since an analyzer was being used 
that was not geared to "Powers ol as are typical 
FFT implementations . Thus , the RTA would be 
used to perform the PSD respo;?se measurements 
from the control accelerometer and software used 
tc perform the inverse FFT in the computer. 

MARKET INPUTS GUIDE DESIGN PHILOSOPHY 

The Basic System requirements as viewed 
by SDC Marketing after observations and dis- 
cussions with various customers were: 

1) Personnel cutbacks in Aerospace and 
Government laboratories resulted in very high 
workloads on system operators , and many new 
and less experienced people were called upon 
to operate complex systems . As a result, operator 
errors and setup errors of systems were on the 
increase . 
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SOFTWARE CONSIDERATIONS 


Figure I * SD30^j Digital Signal Processor 


2) Commensurate witli simple operation 
. ..d setup is that package protection considerations 
must be given high priority to remove as much 
operator decision making as practical* q 


The Analyzer/ Computer system uses a* 
Digital Equipment PDP-it 16-bit minicomputer 
with the following features which make it ideal 
for use in a digital control system: 

a ° 

a> Asynchronous operation * System com- 
ponents run at their highest possible speed* 
replacement with faster devices means foster 
operations without other hardware or softvfare 
changes * 


WP- 


I b) Modular component design resuu^ng 
ease and nexibillty in configuring systems* 


3) Experience with our Computerized 
Spectrum Processing System and Campbell Dia- 
gram Systems indicated that use of a te?eiype 
or other ke ''board input device is very difficult^^*^ 
for the average operator: therefore, every .jffori 
should be made to minimize or eliminate depin- 
dence on a teletype. 


Direct Memory Access capability. 


d) Automatic priority interrupt system . 
of 4 levels permits grouping ' interrupt lin \ ^ 
according to response reqvu,cments , 


4) To give test conductors , inspectors , 
and design engineers total '^visibility" during 
a test by visual observation of the Response 
Spectrum on a real time basis * It was noted that 
many customers owning Random Control systems , 
either analog or digital, were using real time 
analyzers for on independent verification of 
equalization due to the lack of a real time display 
from existing systems . 


e) Vectored inloi-rupts ^ yit^lding fast 
interrupt response to system periphe'*nls * 


Automatic power fail and restart opera - 


g) A single com n uni cations path CE^NIBUS) 
structured for easy communication with a wide 
variety of peiipheroJs and special interfaces* 


5) To make a system as inexpensive as 
practical so as to be attractive as a replacement 
for an older analog system . 




In addition to these hardware features * 
PDP-H's are used in SDC systems because of 
their software support , accepted reliability and 
performance, and a worldwide service system* 


^^^6) To make a system w^th equalization 
speeds equal * or foster . than the analog systems * 
to allow test pro files to be equalized conservative- 
ly a: 1 *5 dB , over a dynamic range of greater 

than 60 dB * 


Software can he roughly categorized as 
high level (FORTRAN, FOCAL, BASIC, etc.) 
and low lev'ei (machine language, assembly lan- 
guage). Relative to low* level languages* high 
level languages generally: 


/ 7) To retain expansion capability for shock 
synthesis, and consider sine and, eventually, 
modal analysis by adding additional hardware. 


Require iess program development 
time; 


Are easier to modiiV; 


To make the system such that those 
wishing to do R ii D evaluations of fixtures , etc* , 
, using random or psuedo raudom could do so * 


9) To use the stand-alone hardware ap^ 
p roach, enabling separate use of the Analyzer 
as a data reduction tool srd the computer for 
calculations. Also, hardware approach would 
allow independent dual real time disph*/ of the 
averaged and unaveraged Response Spectrum. 
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c) Require more memory sTorugo^ and 

d) Execute slower. 

An optimized blend of a high level DEC 
language called FOCAL , and low level assembly 
language is used in Ihd system. This blend yields 
_ a software package Iha^ , while designed to perform 
a specific task, is extremely flexible and fast. 


FOCAL (Formula CALculalor) is an easily 
learned, highpowered, compact programming 
language . It is an interpreter level ionguage , ^ 'V| 
which means that a program mav he written, ' 
debugged, modified, and execute.: "on-line** 
without the usuSuy lengthy "compile" process 
required of compiler level languages , such as 
FORTRAN • 

“ n 

FOCAL is similar to BASIC in its use and 
structure with two si.^nificant exceptions . FOCAL 
permits line and charactcr-by -character editing 
(modifying) of program stalements. This edit 
feature, along with the "interpretive" nature 
of FOCAL, reduces program development tim^f 
signiticanlly , Secondly * the FOCAL interpreter 
only occupies appi\>ximately 2400 words of core C 


os opposed to BASIC which occupies aim'^^st 4000 
words of core . D 


FOCAL is used to perform tasks that do 
not require optimum speed, but require greatest 
flexibility’ These include: 


a) User dialog/ com mumCdUon with the 
system - FOCAL provides a simple , efndem 
means f>r transferring teat information into the 
compute^ - 

b) Report headings and data tabulations 
Assuming the system output device is a low to 
medium speed device such as a teletype or printer 
FOCAL is sufficiently fast to handle output 
This frees tne programmer from having to write 
simple, yet tedious and cumbersome data formal 'O 
conversion rouUnes , output routines , etc 

I A c) Overall system co'Srol and sequencing - 
FOCAL can call on assembly language routines 
to perform high speed input /output aperations* 
and high-speed "number crunching” fsuch as^ 
an FFT> routines. Since FOCAL can pass infot 
mat ion to these routines na parameters . and 
con be easily and quickly modified to call on 
these routines in different sequence, overall 
system operation remains flexible , yet fiist* 

d) "Modeling" of program tasks - In the 
early stages of program development , many pro 
grams can be wrllten and tested in FOi'aL, and 
once proved feasible , convened to the faster 
assembly language version ^ For exomple, a 
new algorithm for rr^mbining narrowband data 
into 1/3 octave data can be tested and debugged 
almost completely ai the FOCAL level whore pro 
gram creation, modificatign, and evaluation is 
desired. O 


Figure 3. SD30LD-C Real Time Analyzer 
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The assembly language portions of the 
software are the workhorses of the system. 
Although usually most difficult to write and more 
time consuming to debug than the FOCAL routines . 
these program sections handle the high-speea 
input/output functions and data manipulations 
that are beyond the capabilities of most high 
level languages . These routines are usually 
modular in form, and are made as general and 
flexible as possible without sacrificing speed 
and memory allocation. This usually results 
in assembly language modules that may be used 
in many systems and/or applications with little 
or no modification necessary . 

Figu’“e 3 shows a typical group of FCCAL 
statements whereby dialog is set up with the 
operator of the Systei.i. 

4.05 S D=FN(PAi;,12);S DN=FX(1 ,FN (DST .10)) 
-D^200 

4.10 D4.95 

4.15 T "DEFINE?^*;D 5; I(X-0Y)4.2,4.5 
4.20 T "UST?^’;D 5; I(X-0Y)4. 25 ,4.6 
4.25 T "EDIT?”;D 5;I(X-0Y)4.3,4.7 
4.30 A "PLOT RANGE (DB): " .X; I(X)4.45 . 

4.45; X FN(PLCT8.X) 

4.45 T "DISABLE TEST?"; D 5; I(X-0Y)4.9 , 

4. 8. 4. 9 
4.5C D 20;G 4.9 
4.60 D 21;G 4.9 
4.70 D 22;C 4.9 
4.80 FN(PAN,2.D);D 20.4 
4.90 X FN(PAN.8) 

4.95 1 •/. 2 , 1 . "TEST NUMBER: " ,D , 1 

Figure 3 . Typical Group of FOCAL Statements 

3.10 X FN(PAN.6,0);X FN(PAiM0,4) 

Figure 4. Typical i.'^rogram Control Transfer 

Figure 4 shows typical "calls" or transfer 
of program control to high-speed assembly lan- 
guage routines . As an example , the routine , 

FN (PAN. 6. 4) has one parameter. The FN tells 
the program that it is a user-defined function; 
the PAN, 6 tells the program which function; 
and the four is the parameter. This particular 
function defines at what pretest level the system 
switches from psuedo random to random drive . 
When the assembly language program has been 
executed, the routine returns to FOCAL and pro- 
ceeds to execute the next sequential program 
statement . 

With a compact interpreter language, opera- 
tor dialog can be rapidly constructed fcr entering 


test profiles and limits, as well as listing and 
editing test profiles . Everything connected with 
the use of the teletype is as simple and short 
as possible for the average user , but for the 
more sophisticated user . there is the ability to 
change certain parameters as desired for special 
applications and conditions . For example . the 
dialog can be modified so as to be in Spanish . 
German. Frencn, etc. .averaging times can be 
changed, psuedo random can be used, etc. 

SYSTEM CONTROL PANEL 

Since much of the effort put into the SDC 
Digital Random System centered around the re- 
quirement that the system be safe and easy to 
operate, a new instrument was designed to ser\'e 
as the system control panel . as well as to house 
certain system pe»i|/herals . such as programmable 
amplifiers, output low pass filter, and system 
timing. As a matter of fact , oac*^ the test profiles 
are entered into the system . there is little need 
for the operator to concern himself with anything 
other than the control panel and the scope display . 

By referring to Figure 5 . it can be seen 
that a different approach than most of the other 
systems that were available on the market w'as 
taken. The test operators may enter tests and 
recall tests by simply assig^iing a number to 
a particular test and then using this number, 
recall and run those tests usin^' the "Test No" 
selector. Provisions for inserting ten different 
test pre files and then recalli; g. running the 
tests, plotting the results, etc. . without any 
further intervention of the teletv’pe were made. 
Thus , the teletype can be removed from the syster-. 
and stored until additional test profiles or new 
profiles are desired to be entered. The control 
panel is sectioned in such a manner as to make 
its use as simple as practical . 

The teletype button, which is on *he lefthand 
side of the control unit in the first section . is 
used to initiate dialog with the teletype connectec 
to the .^vstem . Below that . is a Post Test Data 
switch whereby different responses and records 
made during a test can be recalled , displayed , 
or plotted. The next section to the right is for 
calibration of the system and the actual process 
of starting and running a test . The user can 
define up to four levels below full test so that 
equalization can be accomplished in steps. The 
present pretest-level condition of a test is indi- 
cated by the LED indicators. This is particu- 
larly helpful if many non-linearities are encoun- 
tered or if there is oncern over equalization 
dynamic range. 
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The toggle switch in this section is used 
for several purposes; 

a) During calibration of a test* the toggle 
switch can be used to adjust the full test level 

if it is desired to ran the same profile shape 
at a djGTereni full test level; 

b) After starting a test* the toggle switch 
can be used to manually control the lest level 
through tliC various pretest levels; and 

c) After a test is terminated, the toggle 
switch can be used to set the plotter Y-axis span 
in dB and to select whether or not the plots are 
'auioranged'* or plotted on the same scale as 

the standard spectrum . 

Also in tills section is the test number 
switch whereby □ particular lest profile can be 
selected* recalled and run . Below the test 
number switch is the system STOP button. At 
any time* the test can be lerminated by the STOP 
button . 

The •"Status*’ section in the center of the 
panel tells vfhether the test has aborted, whether 
the system is not equalized within tolerance, 
or whether additional output gain is needed to 
reach the full test level. 

The righthand section gives the operator 
a block diagram of how the sy* era is set up and 
the condition of the test at any particular time. 

The LKP '^ttdout o the lower righthand side 
of the panel is utilized for a multitude of functions. 


During a lest, they indicate the actual ov^crall 
of lest * but at any time , this 
^rms ean be over-ridden by pressing 

the time button and the actual elapsed lime of 
a lest will be displ^Li. 

During a run , if for any reason the lest 
^aborts or is stopped, the LED indicators provide 
IS different stop codes to inform the operator 
Vfhy the test stopped* These codes can be seen 
below the control panel on an ttuxiliary panel 
that is labeled, *’Stc^ Codes Also, on the fame 
panel on the left* are the Post Test Data corthtions 
and the data that can be retrieved from the system 
after a test has been terminated . When plotting 
data* the LED indicator tells the operator what 
g^.'Hz value to label the upper decade of his 
plot, 

SYSTEM CONFIGLRATIOX - Conceptual Block 
Diagram 

Figure S illustrates the conceptual block 
diagram of the system, A brief explanation of 
signal flow through the system follows. 

The output of the signal conditioning ampli 
fiers is fed directly into a programmable attenu 
ator. This programmable attenuator, which is 
under control of the minicomputer, has a primary 
function to optimize the input signal le’vel of the 
Real Time Analyzer so that the s>'siem is always 
operating in an optimized dynamic range condition 
for the analysis of the response pickup. This 
is a very important consideration since dy.mmic 
range , when analyzing random datu. is very 
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Fi.’ure 6. Systerr. Block Uia^rarr. 


important. One of the most common mistiikes 
made when doing random data analysis is not 
to have the input amplifiers of the Real Time 
Analyzer up to their masdmum level. 

The outputs from the Real Time Analyzer 
go in two directions . One output is fed to the 
spectrui.i averager and then to the display oscil- 
loscope. The Real Time Analyzer used has the 
ability to simultaneously display both the averaged 
and the unaveraged spectrum data on the oscillo- 
scope, thus alerting the operator for sudden and 
unexpected changes in the spectrum which . due 
to the averaging process . would show .some time 
later on the screen. This unaveraged display 
capability , while new to most operators , is an 
extremely valuable tool for observing sudden 
changes in response spectrum . 

Another output from the Real Time Analyzer 
is from the Analog-to* Digital Converter that con- 
verts the spectrum amplitude into digital data 
output to the minicomputer. This data is fed 
into the computer memory , averaged . and used 


as the basis for doing the follow-on calculations. 

.*\fler the input data i.s compared to the 
reference spectrum, limit checked, etc. , a random 
time function is fed from ihe minicomputer through 
a direct memor\" access channel and a digital 
to- analog convertor. This random time function 
is fed through an additional programmable attenu- 
ator whose output is passed through a low pass 
filter to remove harmonic components that might 
be generated in the output and then is fod on 
to the power amplifier and into the bhakor. The 
programmable output attenuator adiusts the cor- 
rect overall rms level into the power amplifier. 

.Additional inputs and outputs from the 
System Control unit are for sensing and con 
trolling the indicator light.s of the panel, and 
for supplying power to the panel. .Also contained 
in the System Control unit is a hardware test 
timer and provisions for external hardware abort 
connections. .All the essential design require 
ments for the system arc contained in this con 
ceptual block diagr im 
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THEORETICAL CONSIDERATIONS 

The equalization soeed of any random shaker 
control system is determined by two main factors: 

X. The time it takes to attempt one equali- 
zation; 

2. The number of equalization attempts 
necessary to achieve equalization. 

A certain amount of trade-off can be exer- 
cised be ween these two points. If one attempt 
of equalization is based on data with low statistical 
confidence, the speed of one equalization attempt 
is high , but a higher number of attempts is neces- 
sary* lo achieve equalization: plus the fact that 
tti® resulting equalized spectrum will have a 
distributed mean which is very* undesirable. 

The Random Shaker Control system minimizes 
the number cf attempts; therefore, input data 
is needed with a high statistical confidence. 

We know from various references the statis- 
tical accuracy* of a PSD estimate is a function of 
the effective noise bandwidth of the analysis 
filter and the time over which the output of the 
filter is averaged. Mathemotically , this can 
be expressed as: 

% € = ^ X 100 . where B is the . . 

nBT Noise Bandwidth of the 
analysis filter and 
T is the total averaging 
time in seconds . 

Or, 

c dB = 20 log / 1 + _J \ . (2) 

' sBT / 

The above for.nula is valid for <: < ± 20 

Since we decided we wanted an f of 
approximately 12% or t 1 dB , we could use the 
above formula. 

In the Random Shaker Control system , 
the Real Time Analyzer produces a 500- line PSD 
estimate of the accelerometer signal every 50 
msec, or twenty 500-line estimates per second, 
independent of rate at which the input signal 
is sampled; that is. independent of the bandwidth 
of analysis. The RTA samples the input signal 
at 11 times the Nyquist rate (equal to 3 times the 
full scale frequency range); therefore, the 500- 
line PSD ertimate is based on the 1500 sampled 
values . As long as it takes less than 50 msec 
to collect the 1500 sampled values (one Memory 
Period) , the PSD estimates produced at 50 msec 


rate are independent of each other. However, 
in practical cases the memoiy^ period is longer 
than 50 msec. At 2 kHz for instance, the memory 
period is 250 msec, resulting in 5 PSD estimates 
per me nory period or about 20% independence 
of successive PSD estimates . As each PSD esti- 
mate is produced . it is averaged with previous 
estimates . We have termed this redundant aver- 
aging. 

Although redundant averaging does 
not contribute to the accuracy of the PSD esti * 
mate faster than non - redundant averaging, it 
does reduce the variance of the estimate faster 
by reducing the measurement noise due to instru- 
ment noise and the beating of nearby frequency 
components. The accuraty* of the PSD estimates 
in the Random Shaker Control system can be 
further improved by choosing 1> 2 or 1 '4 of the 
measurement resolution for control resolution . 

In these cases the number of control filters is 
reduced by a factor of 2, or 4. by combining 
2. or 4. PSD frequency* components into one 
filter . The result is that the bandwidth of the 
control filters is effectively doubled or quad- 
rupled . 

Csing Equation 2 for 125- line operaiion 
of the system, a bandwidth of 4 x noise band- 
width of the 500- line ana*ysis is used since aver- 
aging 4 of the 500- line niters gives the reduced 
resolution of 125 lines. The effective noise band- 
width of the 500- line analysis is approximately 
equal to 1.5 X the filter spacing. On the 2000 Hz 
range . the filter spacing is 

= 4 Hz. 

500 

BW = 1.5x4 = 6 Hz Noise B€uiuwidth. 

Thus, when four filters are a. aged, 
a BW for 125-line opeiotion equal to 24 Hz is 
obtained. 

Using the criu • 12% accuracy' of the 
PSD estimate. Equation . used to determine 
the averaging time. 

Solving for 

T = ^ x lO"* = 1 X 10^ -2.9 sec. 

^2q 3456 

Since redundant averaging w*as used . 
or a spectrum is averaged every* 0.050 seconds, 
the closest value to 2.9 sec that was a binary* 
multiple of 0.05; i.e.. 2 x 0.05. 4 x 0.05, etc . . 
is chosen (64 x 0.05 = 3.2 seconds). 
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The reason for a binary multiple is for 
softw'are simplicity in dividing Che resultant 
averaged values by the number of averages . 

Thus . the resulting accuracy of the PSD estimate 
can be computed using Equation 2. 

T = 64 avg. X 0.05 s:>c = 3.2 sec, 

t c dB =20 log I 1 ^ - J = to . 8 dB. 

\ s'lTxXa / 

Another expression of statistical confidence 
is the number of statistical degrees of freedom 
based on a Chi Square distribution; or. 

X = 2 BT = (3) 

2 X 2^ X 3.2 = 154 statistical degrees 
of freedom . 

It has been generally accepted in the field 
that 128 degrees of freedom will usually provide 
accuracy* of t I dB in a PSD estimate. 

Thus for 125-line operation, the PSD esti- 
mate takes 3.2 seconds. The PSD estimate takes 
6.4 seconds for 250 lines and for 500 lines takes 
12.8 seconds , maintaining in all cases the same 
PSD statistical accuracy . 

The loop time, or the period of one equali- 
zation , is defined as the total time that the system 
takes to: 

a) Make the PSD measurement; 

b) Compare against the reference spec- 
trum with limit checks; 


c) Generate a new drive spectrum; and 

d) Inverse FFT to create new* drive time 
function . 

This is kept as short as possible so that 
the PSD measurement will be the major factor 
controlling speed. 

Another desigpn decision made was that 
once equalization is achieved within 1 1 dB a 
new drive spectrum or output record is no longer 
created; thus . removing a large portion of the 
major overhead in computation and approaching 
more closely the averaging time as the controlling 
loop time consideration. 

Also decided was that each equalization 
attempted would correct 100% of the error between 
the measured response and the reference response. 
Thus , the following equation is used: 

New Drive Spectrum = 

Reference Spectrum ^ Spectrum (4) 

Response Spectrum 

This computation is carried out point by 
point on the spectra involved. Figure 7. With 
a perfect PSD measuremer’ and a linear mechani- 
cal system , this would reduce the error betive^n 
measured spectra and desired spectra to zero 
after only one loop time . 

FINAL SYSTEM CO.\ FIG L RATION 

Figure 8 shows the final system configura- 
tion with the display oscilloscope , the control 
panel, the Real Time Analyzer and the computer. 
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Figure 7. F^rror Correction for Measured Re.spon.se 
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Figure 10, Esiientiul Controls to Run a Teat 


Of significance is Figure 11 showing the 
display modes of the Real Time Ana'y^er with 
its cursor* dual or single cL ^play. Having the 
curarr eliminates the problem of putting cali- 
brated gra icuLes on the X axis of the scope and 
p\x>vides a ropid meinod to read the frequency 
of some point of interest , 


Q Figure 12 shows a typical example of how 
a test spccificntion can Ije presented to a system 
operator* The tolorance/abon and standard 
spectrum can be entered using discrete frequency 
and PSD points, or by initial slope and final 
slope plus discrete frequency and PSD points* 
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Referring to Uie lest dinlog shown in Figure 
13, ihe first line of the dialog is the question. 
'^Output limit dB relative to 1 Volt?” This Is 
a question that ivhen answered, restricts the 
level of the output voltage to the power amplifier 
to a safe level in the event of no feedback from 
the response pickup . 


During the startup phase of the system 
the system looks for one of two possible con- 
ditions to exist: 


1) The system output will Increase until 
an overload on the Real Time Analyier is delected 
indicating sufficient signal ivith ivhich to operate; 
or. 


Figure 1 1 , RTA Display \lodes 


2} It will simply Increase the output 
level to a limit that is specified in this first ques 
tion of the dialog. 


Figure 13 .shows the dialog required to 
enter a test. The example shown is for the speci 
hcation given in Figure 12 . 


lines No. 2, 3. and 4 are self-explanatory 
Line 3 asks Automatic mode ”Yes” or ’’No"? In 
the Automatic mode, once the start button has 
been depressed , the system will automatically 
progress through the specified pretest levels. 

It will equalize at each level, failing to equalize 
at each le%*^el , it will make a preselected nun. her 


Once a test is entered , it can be v'erified 
by ^listing” as shovrn in Figure 14. Additionally 
if a mistake is noted or a modiil cation is desii^d 
the test can be edited as shown in Figure IS. A 
final listing is shown in Figure 16. 


test profile 


1,5 dB TOLER A^J^. 


frequency (H/1 


Figure 12. Typical Teat Specification 
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of attempts equalization at each level, then 
automatically progress '.o the next level. This 
process is continued until the full test level 
is reached. The test is then run until the timer 
has been satisfied, or the operator terminates 
the test . 


TEST N^TSERt 0 
DEFINET (Y/ll)t Y 

I. OtnPUT LIMIT (D8 RE 1 VOLT>t 20 
2* TEST OURATION (SEC) (fCOirrillUOUS>t 120 
3. CONTROL RESOLUTION (LINES) i 125 
4 * ACCELEROMETIH SENSITIVITY (HV/6>t 100 
5« AUTOMATIC MODEt (Y/N>l Y 
6« ENABLE LOOP OPEN FUNCTlONt (Y/N)t Y 
ARE YOU SURE! (Y/N)t Y 
7* NUMBER OF PRETEST LEVELSi 2 
PRETEST LEVEL (D3)i 6 
PRETEST LEVEL (DB)t 3 
6* RMS ABORT (DB)l 2 

9. TOLERANCE « ABORT PROFILE 
NUMB01 OF BREAKPCINTSs 2 

BREAKPOINT (HE)t 16 
UPPER TOLERANCE (DB)t 3 
LOVER TOLERANCE (DB)t 3 
ABORT (DB)i 6 

BREAKPOINT (HE) I 2000 
UPPER TOLERANCE (D6)t 1.5 
LOVER TOLERANCE (DB)t 1*5 
ABORT (DB)t 3 

10. STANDARD SPECTRUM 
INITIAL SLOPE (DB/0CT)t 0 
NUMBER OF BREAKPOINTS! 6 

BREAKPOINT (HT)t 8 
AMPLITUDE (G«6/NE)i 0 

BREAKPOINT (Mt>l 16 
AMPLITUDE (6«6/HE>t .0005 

BREAKPOINT (KZ>t 150 
AMPLITUDE (6«6/HZ)l .0005 

BREAKPOINT (HZ)t 300 
amplitude (6«6/MZ)i .015 

BREAKPOINT (HZ)t 700 
AMPLITUDE (G«6/HZ)l ,015 

BREAKPOINT (HZ) I 1000 
AMPLITUDE (6«6/HZ)t ,005 
FINAL SLOPE (DB/0CT)t 0 

Figure 13. Initial Dialog 


Line 6 asks whether or not to enable the 
Loop Open function . If the response to that is 
"Yes”, it asks the question, "Are you sure?" 

At full test level, by pressing the Test button 
an additional time , the feedback loop can be 
opened and the system continue to output the 
last equalized drive spectrum into the package. 
This is particularly useful for doing fixture 
evaluation and looking at responses of other 
points on a package . The only danger in this 
is that all the software limit checking and pro- 
tection IS disabled; thus, the system has no 


software protection and if no hardware protector 
is present , then the system would be totally with 
out protection . 


The Loop Open function is provided as 
a convenience for those who wish to do evalua- 
tions and look at responses from other points 
of a package. In essence, it completely frees 
the Real Time Analyzer for various other func- 
tions, even looking at responses on different 
frequency ranges . As a matter of fact durin,*^ 
operation of the system . the Real Time Analyzer 
and the Averager have been disconnected and 
removed from the rack . 


Line 7 , Number of Pretest Levels , is used 
to define up to four different pretest levels , which 
can be defined at any level relative to full test 
level . 


TEST NUMBERt 0 
DEFINET (Y/N)l 
LIST? (Y/M)f Y 
LIST ALL? (Y/M)l Y 


TEST NUMBERt 0 

1, OUTPUT LIMIT (08 RE 1 VOLT)i ' 20 

2, TEST DURATION (SEC) (0-CONTlNUOUS) t 120 

3, CONTROL RESOLUTION (LINES)! 125 

4, accelerometer SENSITIVITY (HV/G>t 100,00 

5, AUTOMAT tC MODEt YES 

6, ENABLE LOOP OPEN FUNCTION! YES 

7, PRETEST levels (DB)t 

- 6 
- 3 

0. RMS ABORT (DB)t 2,0 
9. TOLERANCE A ABORT PROFILE 


BRKPT 

♦TOLER 

-TOLER 

♦ABORT 

16 

3.0 

-3,0 

6,0 

2000 

1,5 

-!*5 

3,0 

lb, STANDARD SPECTRUM 


INITIAL 

5iL0PE 

(DB/OCT)! 

6,00 

FINAL 

SLOPE 

(DB/OCT)t 

0,00 


BRXPT (HZ) AHPLT 

6 0,0000 

)6 0,0005 

150 0,0005 

300 0,0150 

700 0,0150 

1000 0*0050 


Figure 14. Initial Listing 


Line 8, Ask for an RMS Abort in dP^. 
if the test level is exceeded by x dB, the sys- 
tem will shut the test down. Undertest, how 
ever, is not protected except for loss of feed 
back . 
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Line 9 is fairly self-explanatory. The 
tolerance and abort profiles are entered, the 
same as a standard spectrum is entered. The 
upper and lower tolerance in dB is referred to 
the level of the standard spectrum . The abort 
level is an ima^nary line that rides on the top 
of the upper tolerance profile and if any parti- 
cular spectral line exceeds the specified dB 
relative to the standard, it wiil shut the lest 
down. 

Line 10 is the dialog used to enter the 
standard spectrum . It , once again . is self- 
explanatory . 

TEST NUHB0tt § 

DEFlNEt CY/mi 
LIST? <T,3i)i 
IDIT? CY/IOt Y 
LINE NUHBERt 3 

3. CONTROL RESOLtrriON CLINES) i £S« 
LINE NUHBERi 7 

7. NtmBER OP PRETEST LEVELSt 3 
PRETEST LEVEL CDB)l 9 
PRETEST LEVEL (DB)l 6 
PRETEST LEVEL (DB)I 3 

LINE NUHBERt 

Figure 15. Editing 


TE:>'^ NURBIPi • 
UEPIUr* .Y/N)l 
LIST? <Y/N)I Y 
LIST ALL? CY/W)i Y 


TEST NtJHBfitt f 

1. OUTPUT LIHIT CD8 RE 1 WOLTji - 2 § 
e. TEST duration CSEC) (0«CONTINUOUS)t 120 
3* CONTROL RESOLUTION (LIN£S>i ^S0 
A* ACCELEROHETER SENSITIVITY (NV/6)l 100.00 

5. AUTOMATIC MOrEl YES 

6. ENABLE LOOP OPfeV PUNCTlONt YES 

7. PRETEST LEVILS <DB)i 

- 9 
• 6 

- 3 

8. RMS ABORT CDB)t 2*0 

9* TOLERANCE A ABORT PROFILE 


BRKPT 

♦TOLER 

-TOLER 

♦ABORT 

16 

3*0 

-3.0 

6*0 

2000 

1*5 

-1*5 

3.0 

10. STANDARD SPECTRIW 


INITIAL 

SLOPE 

(DB/0CT)l 

0*00 

FINAL 

SLOPE 

(DB/0CT)l 

0.00 


BRKPT fHZ) AMPLT 

8 0*0800 

16 0*0005 

1^0 0*0005 

300 0*0150 

700 0*0150 

1«00 0*0050 

Figure 16. Final Listing 


Put yourself in the place of an operator 
and take a look at Figure 17. How do you run 
a tesi? 



C PLOTTER CAL 
D F/SPSDCAL 
E TEST TIME AT ABORT 


Figure 17. 7 ypical 'I'est Sequenc e 
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First, assume that a test profile has been 
defined and assig^ned the test number 5. Dial 
”5” into the TEST NUMBER switch and press 
START. If a system calibration is necessary , 
the CAL light turns on, and the system does 
a fall calibration and functional check. When 
this is complete, the STOP light comes on. 

Now , pressing the START button turns on the 
START lamp and starts the full test sequence. 

Note that the system will not allow the test to 
proceed until a proper calibration is performed. 

If the same test is run repeatedly , a CAL is neces- 
sary only for the first run. This is the type 
of operational "interlocks" built into the system 
software to insure safe operation. 

If tlie Automatic mode in the TTY dialog 
is specified for test 5 , the operator simply sits 
back and watches. The system now equalizes 
at each of the pre-defined pretest levels . As 
equalization is accomplished, or with a preselec 
ted number of tries to equalize (if equalization 
is not achieved within the test tolerances) , the 
system steps to the next pretest level and again 
attempts to equalize. Following the last pretest 
level, the system increases the drive level to 
full test level, and turns on the TEST lamp. 

From the operator’s standpoint, he must 
keep in mind only two rules: 

1 . He can never start a test without inten 
tionally pressing START; 

2. He can stop a test at any time by press- 
ing STOP. 

If the test *s to be timed, the system auto- 
matically shuts off at the end of the test time and 
indicates an 05 in the numeric indicator. If the 
system should abort, the numeric indicator shows 
the operator the reason: 

01 External Abort #1 

02 External Abort #2 

03 External Abort #3 

04 External Abort #4 

The External Abort signals are simply contact 
closures which are connected to the system 
through the control panel. Each pair of contact 
lines may be jumpered .or Normally Open or 
Normally Closed operation. These signals would 
typically indicate system "status" lines such as 
shaker overtravel, power amplifier overheating, 
or hand-held push button switches held by test 
observers . 


06 RMS Abort (the rms level has 
exceeded the abort limit) 

07 Line Abort (one or more spectrum 
amplitudes have exceeded the 
abort limit) 

08 Loss of input signal 

During the test , the operator is informed 
about what the system is doing through the Real 
Time Spectrum display and the system front 
panel. The oociHoscope display gives a con 
tinuous 'uok at the actual power spectrum 
achieved ru the control accelerometer. The LED 
numeric readout on the system panel shows the 
grms accelerometer. If the operator 

wishes, he can press the TIME button, and mo- 
mentarily override the gj.j^g display and show 
the elapsed time for the test . Releasing this 
button again shows the iovel. 

Three indicators tell the operator how 
accurately the system is achieving equalization. 

If the yellow GAIN light is on, the system (includ- 
ing the system power amplifier) needs more 
gain to reach full test level . This light is active 
during all startup modes (START and PRETEST) 
as well as in the TEST mode , to alert the operator 
to increase the system gain , either with the control 
pemel OUTPUT GAIN control , or with the power 
amplifier geun control. 

The green TOLERANCE light tells the opera- 
tor that the system is not equali zed to within 
the specifications called out in the input dialog, 
^he red ABORT light signals any premature 
interruptioii the test. 

During the test , the operator may choose 
to store PSD spectra at any lime lor plotting at 
the conclusion of the test. He does this by press- 
ing the POST TEST DATA button . The spectrum 
thus stored is assigned a number which appears 
momentarily in the numeric indicator. The POST 
TEST rotary switch can then be used ai the test 
conclusion to reference and plot these spectra. 

Since the test verification is really the 
"final product" of any vibration test, con.siderable 
attention has been paid to producing calibrated 
hard copy plots of tho accelerometer PSD. Since 
the computer is in control of the system gain 
settings nt all times, the calioration of each plot 
is done automatically , and the operator is alerted 
as to tho lull scale value of each plot via the 
numeric indicator on the control panel. 
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EXPANSION CAPABILITIES 
Shock Synthesis And Analysis 

Expansion of the system for shock work on 
shakers iz accomplished by simply substituting^ 
either the standard Model SD320 or SD32x Extended 
Memory Shock Spectrum Analyzer for the Real 
Time Anal> zer in the Random System and loading 
new software into the PDP-11 Computer. The 
system has three decades of frequency coverage , 

10 to 10»000 Hz, 1 to 1,000 Hz and by using the 
Extended Memory , 0 . 1 to 100 Hz for seismic and 
transportation work. Specially designed phase 
equalized aliasing filters are usea to preserve 
integrity of the input transient. A sampling 
rate of four times the upper frequency limit is 
used to permit optimum rolloff and rejection by 
the aliasing filter. 

Block Diagram Figure 6 illustrates the 
simplicity of the interchange of the independent 
analyzers in the system. 

Figure 18 shows the Extended Memory 
Shock Spectrum Analyzer mounted in a com- 
panion rack for mating with the Random System, 

Analysis 

The Shock Spectrum Analyzer consists 
of a digital input memory , an analog section for 
high speed analysis and both digital and analog 
outputs. When not in use for closed loop system 
operation, the analyzer can be used as a stand- 
alone unit . 

Analysis speed depends upon the amount 
of date stored in input memory . Functionally , 
data is analyzed at a rate of 4,000 words in two 
seconds. Primarily for seismic work, the E tended 
Memory version of the analyzer permits selection 
of up to 40,000 words of data in 4,000 word incre- 
ments. Maximum analysis time under this con- 
dition would be 20 seconds . 

The analyzer's performance in the system 
provides extremely fast feedback information 
for the computer to permit rapid re- definition 
of an input time history to the vibration shock 
testing system. The transient stored in ivlemorj’ 
or the shock response spectrum may be viewed 
on an oscilloscope. See Figure 18. 

Synthesis 

The Wavelet Amplitude Equalization (WAE) 
technique is employed for shock synthesis in 
the system. The technique has been well des- 
cribed in technical literature and will not be 


expanded upon in this paper. Shock synthesis 
can be performed on either 1/3, or 1/6, octave 
spacing as defined by the operator. A simple 
dialog sequence is performed by the operator 
to select resolution on either a 1/ .i, or 1/6, octave 
spacing together with the number (up to 0 of 
pre-test excitation levels to be used in coming 
up to full test. The operator selects the frequency 
and amplitude breakpoints of the desired response 
spectrum profile , or individually defines each 
wavelet's amplitude and delay. When the spec- 
trum orofile approach is used , the computer 
calcuiates the least energy required to meet 
the desired spectrum amplitude and automatically 
adjusts each wavelet accordingly. Ab in the 
Random System , 10 test profiles may be stored 
and recalled by front panel control selection . 

Response data fron. the package is fed 
into the input of the analyzer wiiere analysis 
is performed on 1/3, 1/6 or 1/12 octave basis 
at whatever damping nas been seJected for the 
test . Data can be directly output to aii X-Y plotter 
at this point . The digitally stored response 
spectra information is automatically fed back 
to the computer where a new output waveform 
series can be generated for output to the system 
if desired. 

Sine Testing 

Originally , it was inconceivable to SDC 
des7t;nors that anyone would use an expensive 
system of tne type described above to do sine 
testing. We already make several relatively 
inexpensive analog sine controllers that are 
simple to operate. Lately, we are finding that 
more of our customers are intt rested in sine 
control expansion. We have plans to add this 
capability in the future. 

FINAL SYSTEM PERFORMANCE SUMMARY 

Frequency Ranges: 

Standaid - 2 kHz 
Highest Range - 5 kHz 
Other Ranges - Optional 

Control Resoh’tion: 

125, 250, or 500 filter elements 

Coi'e Requirements: 16k per 125 lines 
operation, 24k for 250 and 500 lin.^s. 

Dynamic Range: 

Standard Spec, cannot exceed 30 dB. 
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Output Drive Spectrum: Greater thcLi 80 dli . 

Equalization Speed. 

125 lines - 4.65 seconds 
250 lines - 8.6 seconds 
500 lines - 16.5 seconds 

Max. Output Level: 20 V peak-to-peak 

Statistical Confidence of PSD estimates: 

Program adjustable •• ± 0.8 dB standard. 

Output Waveform: Gaussian Random Noise or 
Psudeo Random Noise with Three Sigma peak 
limiting. 

Dialog: Entering a test rec[uires answering ten 
questions - Standard anc tolerance /abort 
spectra can be entered a.*; discrete frequency 
and amplitude points , or by ii.ltial and final 
slor >s plus discrete frequency and ampli' 
tU' ts . Points are joined together so 

as u jti'aight lines on a log frequency , 
log ampJtur.e scale. 

Automatic Mode. Brings system up automati- 
cally and proceeds automatically through the 
pretest levels up to full test level . 

Loop Open: System continues to output last 
equalized random spectrum. Anah ** can 
be used independently in Loop Ope., .iiode . 


Response Records: At any time during test 
or pretest operation , a record of the response 
accelerometer can be saved. Three records 
can be saved in the 16k system . Up to seven 
records can be made in the 24k system except 
when usiiig 500 -line resolution , only five can 
be saved . 

Limit Checks and Aborts: Stand '^rd spectrum 
is compared with the response spectrum; then 
the profile tolerance comparison is done on a 
relative ’’shape” basis. Abort levels are 
checked on an "absolute level” basis. Noinpu. 
signal, or loss of input signal is an abort check . 

Software: System comes loaded with FOCAT ‘ 
interpreter and all software requireo for 
Random Control Diagnostic software is also 
supplied . 

CONCLUSIONS 

By oeginnii/g with specific guidelines 
for a system configuration, the evolution of a 
simple- to-operaie automatic control «*ystenj nuh 
built-in safety protection has been snown. Sam- 
ple software routines have been given to -?how 
the with which test routines can be estab- 
lished and modified. By using the Real Time 
Analyzer for basic PSD computation, not only has 
test visibility been greatly increased, but analy- 
sis flexibility has been increased through over- 
lapped processing, and system costs kept to a 
minimum using proven, off-the-shelf hardware. 
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TECHNIQUES FOR NARROW R\KD RANDOM OR SINE ON 
WTDERAND RANDOM VIBRATION TFSTING WITH A 
DIGITAL CONTROL SYSTEM 

Michael K. ^lauffer 
Time/Data Corporation 
I^alo Alto, California 
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Raudom Sine, and Shock vibration control and testing programs were previously 
implemented in Digital Control Systems ustn^ the same basic hardware, to take 
full ad\*antage oi this existing digital control hardware, new programs ^re being 
enveloped for testing other environments v -h have been limited to au .og hard- 
ware. Two of these programs now available are *Sv.ept Narr.#wband on Broad- 
band Random Vibration Control*’ and ''Swept Sine on Random Vibration Cmicrol”. 
Although these programs use the same hardware as the basic IXgital Control 
Systems, their uuique environments »*equired special techniques to achieve the 
test criteria The tei^r»niques for implementing these new programs are de- 
scribed here in two norts. 


PART t. DIGITAL SWEPT NARROW R\ND ON 
BROADBXND RANDOM VIBRATION 
CONTROL SYSTEM 


INTRODUCTIOK TO PART I 

I^rt I of this pr er describes a digital narrowband on 
broadband random v;i. ation control system. The sys-. 
tern i'.vs been tmplt.tnented by modifying a standard 
digital vibration control system such that the refer- 
ence s|)ectrum is the superposition of a static broad- 
band spectrum and a sweeinng narrowband Spectrum- 
Depending upon test philo$oph>, superposition is de- 
fined in at least two wa;.*s: 

1. The sum of the narrow and broad spectrum 
amplitudes at each frequency. 

2. The maximum amplitude of the narrow or oroad- 
band spec ' : am (that is, whichever is greater) 

at each frequency. 


in the case of summing the two jpectrums, overall 
GRMS level is constant as the narrowband spectrum 
sweeps, the case of controlling to the maximum, 
the GR\U' .evel varies as the narrowband spectrum 
sweeps. C ontrolling to the sex is obviously a more 
severe test than controlling to the maximum. The 
system described here allows for control to either 
tne sum or the maximum r>f the narrowba'id and broad- 
band reference spectrums. 


GENERAL DESCRIPTION 

The control system i shown in Fig. 1, is b*iilt around 
a standard minicomputer and a 200 kHz data acquisi- 
tion system that serve as the heart of a family of sig- 
nal analysis and shodc and vibration products (1-5). 

The hardware configuntt}*^ is identical to that of a 
staixiard random vibration control vm as show'* 
in Fig. 2. NanxrwbaikJ on hroadban Yration cou 
trol is obtained by executing the appropriate softwai^'" ' 
in the minicomputer. A single output channel of 
shaped random Gaussian noise is generated to drive 
a shake :* system, and a single channel is sampled and 
controlled. 

Opeiator interaction with Che system occurs in two 
ways. A switch control panel, shown in Fig. .1, is 
used to select system functions and to execute control 
commands, while an alphanumeric keyboani is used 
to enter test parameters. 


IM PLEMENTATION 

A functional block liiagram of the control strategy is 
shown in Fig. 4, Basically, the system attempts to 
control the amplitude of a designated control signal 
to match the amplitude of an operator-s|)ecified swept 
narrow'band on broadband reference spectrum by co*'- 
tinuously adjusting the random noise drive signal at 
each frequency. This drive adjustment compensates 
for nulls apd resonances in the shaker system trans- 
fer function. The coi:trol strateg>- is identical to that 
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tn^u-:id of a single n^ference spectrum, three refer- 
ence ^ipectnitti? are maiatiiJned in computer memon* 
Tht^se are the hroatibund, tutrroHK-^imi, and sup-r- O 
imposed spectra. The hroadbood anrl niirrowiand 

are sjwfcifled h% a con^ersaiicsttal setufi rou* 
tine and, once getieniUsI, remain static thrcKighout 
3 test. The suiierlfnp^iei spectmm is the sum or 
the maximum t^lue of the broad and sweeping narrow - 
iKind spei‘tra. The narrowband spectrum is s»e|vt 
:icrc^5 the 1iroaLdha«l sfiecirum by u|xiatiag an index 
i-ariable (INDX) which itenotes the rebti%c position 
of the narrowband and broadl*;iiid spectra origins in 
terms of ^}JCclT3tl lines. The data ealues represent- 
ing the narrow bam I spectrum never actually shift 
position in memor> , Whenever IKDX changes a new 
total reference siiectruin is jtenerated l*y com tuning 
the broadbanti am! the 'shifted narrowband sjw^tra. 
Thi‘ narrowband sj»ectrum Is extended into an extra 
buffer to maintain continuit>* in the combiner I spectrum 
w hen the narrowbawl spectra sweefis lownrtl lower 
frequencieti . 
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used In the stiiodam -- 
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SWEEP TIMING AND FHKtfl flNO CPDATE 


Sweep timing is obtained from the end -of -fra me inier- 
ruj4s gener*UHi the* data acqutsttion and i -uttut 
hardware. Given the desired sweep rale, Utc sam- 
pling nilc, and tlK- numiwr of samples jicr output 
frame, a f resiliency change i«*r interrupt is calcutatcd. 
This frcipjcncv change is positive- for increasing fre- 
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queney swie(JS and neg.ilive for frtntvtene*. 

sweeps. The sign of this fmpiunc\ ebaitse is comple- 
mented each lime a sweep limit is reached. 

The flow chart in Fig* S shows how the swvcf^ iqx-r- 
ates, InterrujJis are aecumt :lale<l in a counter vnri* 
able* vbice each loop time a fn"'|ueacv act'umulaior 
preseni frequency of the right most 
narrow liaiid spectrum breakpoint is u|jitite"i bv a val- 
ue equal to the numtier of inlerruptf: times ibe fre- 
quenc\ change per Interrupt, The interrupt accumu- 
latf^r is I’ en reset to /en>* An integer \nriable 

OCFWTKlyi is updated hfvt 


:zmisu 


Right-most narrow breakpoint fregi^encv 
Frequenct change p»*r s|w-ctml line 


If NfTVDCU changei^, so must thcsuperlm|<ose*l spec- 
im, aiHt the shifted narrow s^iectrum Is combineti 
with the iinv'idlKiti^l s|iecirum to genent^y a new super- 
imposet! reference spec^a. The va,Tnble LVt,rA <tl»- 
cusscti previouslv is alST uplated whenever X fAV INlJ 
changes, ^ 

When sw'ecinng townirfb higher fret;uenclef?* the fre^ 
i|uenc% at tltt* highest narmwk.twl sjwclrum bn^ak- 
point iS com|cir *<1 to the ojwrHitor-s|<?ci fil'd upi^r 
weep limit, 11 it ci^b or excec^ds this iimiH 
sweep direction ts reversed. When swccfiing titwarrl 
Iwer frequet'cies, the frequenev of the lowesi narrow' 
binl sriertrum breakpoint Is compared to thi- h:mer 
sw'tH^^p limit. If it is less likin or c<|ual to this limit, 
sweep dircctieu is reverst?tl* 
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TYPICAL SYSTFM OPFRATION 

Parameters of a typical swept narrowlxind on broad- 
band random vibration control u*si are sho%*n in 
Fig. R. Xoie that specification of the narrowband 
spectrum is similar to broadband spectrum setup 
with the exception of alarm limits. Note also the 
superposition option of SCM or XL\X and the LOG/ 
LINEAR sweep and D’^fLAL DlItKCTION options 
under SWEEP PARAMETERS. 


Displays of the combined reference spectrum :** var- 
ious times thioughout a test are shown m Fig. 7, 

The control spectrum foUous this reference sfH>ctrum 
rather close.v during the test. How closeb' depends 
on the sfiectrum shapes anti the sweep rate. 


SUMMARY OF PART 1 


PART n. DIGITAL SWEPT SINE ON 

RANDOM VIBRATION CONTROL 
SYSTEM 

INTRODUCTION TO PART U 

Fart II of this paper e.xplores sc .'oral techniqi'as for 
the implementation of a digital swept sine on ra^om 
vibration control system. These techniques differ 
primarily in t^> cost-performance tradeoffs associ- 
ated with the sine frequence* generation metiiod and 
the sine control strategy. Discussion begins wiui the 
least expensive and lowest performance metnod and 
proceeds through the increasingly expensive and high- 
er performance methods. The limitations and reasons 
foK the limitations of each method are also discussed. 


METHOD I 

This method is the least expensii’e in terms of add-on 
hardware to the basic s\-stem hardware required for 
a flexible vibratimi control system. The hardware 
configuration* shown in Fig. I* is that of a standard 
\-ibration control system (l-G). 

The block diagram in Fig. ^ outlines the basic control 
strategy'. .An estimate of the shaker system transfer 
function is obtained from the input and output power 
spectral densities using the method described by 
Norin^^^. Tnis transfer functicn estimate is used in 
the random conti*ol algorithm exactlv as explained by 
Norin^^^. The sine control algorithm uses this same 
transfer function estimate in a similar manner* but 
it is otherwise independent of the random control 
strategy. 

An estimate of the sine control \*a!je is obtained as 
follows* 

SINE CONTROL H{fs) • SINE DRfA’E 

where SINE ORfYE was the amplitude of the sine 
drive output while Hffg), the transfer function esti- 
mate at the sine frequency, was calculated. 

The sine control estimate is oomixsred to a sine ref- 
erence value which remains constant in amplitude but 
\*aries in frcqueiicy according to the t\ 7 .e and rate of 
sweep. The error between the control and reference 
indicates how to u()date the sine drive value to com- 
{>ensate for the shaker trar*sfer function. 


l^rt I of this paper discussed one implementation of a 
digital swept narrowband on broadband random \ib- 
ration control system. By utilizing the same hard- 
ware as a series of othi r basic vibration testing pro- 
grams (random, shock, sine), this software is a 
flexible id economic supplement to the basic capa- 
bilities of the testing laboratory. 


This technique has the advantage that no Ixmdpass 
filter is required to extrrict the sine energy from the 
sine plus random control signal. The reason for this 
is that the sine control is calculated from the product 
of the known sine drive and the transfer function esti- 
mate. Because of the way it is calculated, the trans- 
fer function estimate is the same v • . ther c r not a 
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multiple of the frequency Increment. 
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E^AHHOVBAIO HEFEREliCE SPECTHUtU 
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alar:: limit DD : 5.23C 

9 FPEOUEI-CY HZ.: 1376. 

LEVEL GSOP/HZ.: .IC33 
ALAm: LIMIT DB : 6.200 

10 FHECLLIjCY HZ.: 1600. 

LEVEL GSen/HZ.: I .030 

alarm limit DD : 6. "'CO 

l\ rnEOU£i:CV HZ.: 1792. 

LEVrL GSOn/HZ.; 1.000 
alarm limit D2 : 6.000 

12 FI Mai. 5L0PE# DB/OCT: -43.00 


15 FHEQUEIiCY HZ.: 992.0 
LEVEL GSQH/HZ . : I .030 


16 FPEQUEliCY HZ.s 1200- 
LEVEL GSOP/KZ.: 1.330 


17 FltiAL FLOPE, DS/OCT: -64.00 


I.ARPOW GRi:S s 17.67 
TOTAL Gm:s : 32.63 

18 SWEEP PARAMETERS: 

NODE !=L0G#3=L1U: 0 
RATE, HZ. /SEC: 6.000 

tLlTlAL OIPECTIOL 1=U?, 


3=DCVM: 1 


19 LOWER SWEEP LIMIT, HZ - ; 43C-3 

20 UPP£^ SWEEP LIMIT, HZ.: 15C0. 

21 LOW Lr^EL, -DB: -23.00 

22 LEVEL IMCREMEIiT, DB: 2.000 

23 start-;? TIME SEC: .6251 

24 SHUT-DOWL time SEC: .9377 

25 TEST TIME HRS, MIM, SEC: 3,4,3 

26 COfiTP.OL CHAIiMELS: 1 

27 ACCEL SEI-S MV/G: 10. DC 

28 DRIVE CLIPPIMC I=YES,0=WO: 3 

29 ABORT LEVEL GRHS: 63.00 

30 ALARI-i LU.ES: 1,3,5 

31 ABORT LIMES: 14,16 


CORRECT! OMS 1=YES 0=MO: 3 
LIST l=Y£S 0=M0: 0 
SAVE I==YES 0=MO: 3 


BROAD GRi:S : 27.43 

Fig. 6. I.isUng of Setup Parameters for Narrow-band on Broadband 


Sint is .'.u|xrrim{K>sed on the random signal. This, of 
course, also means that the sine energN* does not af- 
fect U*e random control signal since it is calculated 
in a similar manner. 

In Method I the sine drive amplitude is adjusted on 
frequency data stored as a buffer of discrete, evenly 
s^xiced sine frequencies. The number of discrete 
sine frequencies is equal to the number of spectral 
lines in the random sjiectrum. The frequency differ- 
ence, :if, between adjacent sines is: 

^ random spectrum bandwidth 

resolution (number of spectral lines) 


The inverse Fourier transform of this sine amplitude 
Ijuffer IS calculated to generate a buffer of time do- 
main data representing the actual sinusoidal time 
wa»*eform. This buffer is added to the random data 
buffer, and the sum is output through the digital-to- 
analog converter. 

SWEEP MECHANISM 

Sweep timing is obtained from end-of-frame interrupts 
generated by the data acquisition and ouijxit hardware. 
Given the desired sweep rate, the sampling rate, and 
the number of samples per frame, a frequenev change 
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Fig. 7. Swept Narrowband on Broadband Reference Spectra 


per interrupt is calculated. This frequency change is 
positive for increasing frequency sweeps and negative 
for decreasing frequency sweeps. The sign of this 
frequency change is complemented each time a sweep 
limit is reached. The flow chart in Fig. 9 shows 
how the sweep operates. 


Interrupts are accumulated in a counter variable. 

Once each loop time a frequency accumulator for each 
sine is updated by a value equal to the number of inter- 
rupts times the frequency change per interrupt. The 
number-of- interrupts counter is then set to zero. 

Next a frequency index variable for each sine is cal- 
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culured from: 


1 HEQ INDEX = 


FREQUENCY ACCUMULATOR 
FREQUENCY CHANGE PER 
SPECTRAL Line (^f) 


rate and (banch^’idth/resolution). The relationship 
ts: 

^UXIMUM SWEEP RATE FOR A SINE 
RESOLUTION OF if = if LOOPTiME 


This index is the new sine frequency (it may be un- 
w anged). 


VThen full test level is reached, the sine waves sw*eep 
synchronously up or down at a user-progi^mmed rate 
(linear or logarithmic) until a limit frequency is 
reached by one of the sine wav*es. The sine waves 
then reverse direction, sweeping at the same rate 
until the opposite limit frequency is reached by one of 
the sine waves. The sine waves again reverse direc- 
tion and continue sweeping back and forth in this man- 
ner until test time expires. 


LIMITATIONS 

Method I has some fundamental limitations on the 
basic sine resolution and the sweep resolution. As 
mentioned previously, the fundamental sine resolution 
is limited to A f as calculated above because the time 
domain sine waveform is generated by an inverse 
Fourier transform. Also, since the sweeping sine 
frequencies are updated once each loop time, the 
sweep resolution (increment) is limited by the sweep 



Fig. *J. Swept Sine on Random Flow Chart of 
Sw'eep Frequency Updating 
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If 1 $ 16 Hz and the loop time is two seconds* then 
the maximum accurate sweep rate - 16 Hz/2 sec - 
8 Hz/sec. That is, the sines will increment by one 
spectral line (16 Hz) each loop time. If a sweep rate 
of 32 Hz/sec were specified with a loop time of two 
seconds, the discrete sine frequencies would step in 
64 Hz increments. Thus, the program provides bet- 
ter swept sine resolution at lower sweep rates. 


The sine frequency is not updated more often than 
once per loop time because of the pmcessing time re- 
quired to perform the inverse Fourier transform. 


HOW TO IMPROVE SINE RESOLUTION 

At least two different sine generation meutods can be 
used to obtain better sine frequency resolution; a soft- 
ware sine table lookup or a hardware sine synthesizer 
(Method n and Method III, respectively). 


METHOD U 

To obtiiin finer sine resolution without additional hard- 
ware, the sine values can be generated by a table 
look-up procedure, as is done in some hardware dig- 
ital sine svnthesizers^^^ Amplitude control is ac- 
complished by multiplication of each sine point by the 
amplitude factor produced by a control strategy sim- 
ilar to that used in Method I. Sine values obtained 
from the look-up table are multiplied by the amplitude 
factor and added directly to the tapered random output 
buffer. This sum buffer is then output through the 
digital-to-analog converter. The sine generation pro- 
cedure is continuous and oblivious to output buffer 
boundaries. Figure 10 shows how this procedure 
works. 

In addition to having increased resolution, this tech- 
nique doeo not require an inverse Fou'»**3r ti'ansfor.a. 
However, a multiplication is required for each output 
point. The number of multiplications increases with 
the numlxjr of sines, thus pi: cing limitations on band- 
width (bandwidth is related to how fast points can be 
output). 

Control strategies for random and sine can be the 
same as in Method I, or modified as explained below. 



NOTES. 

A = Inverse Fourier Tr aosf Of in (IFT) Generation 

B = Table Lookup 

B has twice the resolution of A and can get even better resolution 
Fig. 10. Table- Lookup Buffer Overlap 

since the sine is superimposed on the random in the 
analog domain, no digital informatiuu exists for this 
sum which is the input to the shaker system. There- 
fore, a second analog- to-digital conversion channel 
can be used to sample this analog sum simultaneously 
with the output of the shaker system. This sampled 
data is equi\*alent to the digital sum of the s .le and 
random data output through the digital-to-analog con- 
verter in Method I. 

Therefore, the techniques discussed in Method I can 
be applied to control the sine and random signals. A 
sine synthesizer is required for each s»ne to be gen- 
erated and controlled. 


CONTROL STR.ATEGIES AND SWEEP RESOLUTION 
IMPROVEMENT 

The finer sine resolution of Methods II and Til may 
necessitate changes in the sine control strategy. 

Sine control can be accomplished by several tech- 
niques. One method is to use the transfer lunction 
estimate calculated for random control to generate a 
control estimate as in Method I. This control esti- 
mate is then com^ured to the sine reference, and the 
sine drive is compensated accordingly through adjust- 
ment of the sine amplitude multiplier. As mentioned 
before, the ad\*antage ot this is that no tracking filter 
is required to extract the sine information. The dis- 
adv’untage is that the transfer function estimate has 
much coarser frequency resolution than the actual 
sine. Thus, narrow resonances and nulls may not be 
discriminated very accurately. 


METHOD HI 

Bandwidth limitati'ms imposed by Method II can be 
(j.'jroomc by adding n hardware sine synthesizer to 
the sy.«5tem described in Method I. Figure 11 is a 
block diagram of a single sine system of this type. 


Another method of control uses a tracking filter to 
extract the sine amplitude and phase information. 
This method is explained in greater detail by 
Norir/^K 

Now that two options for improving the basic sine 
resolution and t\\o basic control strategies have been 
presented, several methods for improving the sweep 
resolution will be discussed. Consider Hie sine con- 


It is desirable to use strategics similar to the above 
mothovls for random and sine control. However, 
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NOTES: 

T H = Y Buffer Power Spectral Density (PSD) 

X Buffer Power Spectral Density (PSD) 

2. Broken lines indicate ofserations which could 
be performed by a parallel microprocessor. 

Fi^. 11. Swept Sine on Random with Hardware Sine Synthesizer Block Diagram 


trol strateg>’ using the transfer function estimate 
technique. If the sine frequency is only updated once 
per random loop time (e. g. , when the transfer func- 
tion estimate has been updated), then the same sweep 
resolution exists as in Method I. 

It is desirable to update the sine frequency as often 
as possible. Once e\^er>' interrupt is the maximum 
rate. Updating the frequency at this rate is not dif- 
ficult. The problem lies in the control strategy, 
w'hose calculations can take considerable time and, 
thus, increase the overall loop time. Several strat* 
gies can be used: 

1. The control is updated each time the frequency 
is updated. This increases overall loop time 
and limits bandwidth, bc't it provides the best 
control possible with the most recent transfer 
function estimate. Note that the transfer func- 
tion estimate is updated only once per loop cime. 

2. Control update occurs less frequently than fre- 
quency update. In effect, the sine sweeps part 
of the time without being tightly controlled. 

Thus, the sine control is not as good, but more 
time is available for random control which is 
consequently better. For examnle, the control 
could be uixiated once per loop time, but the fre- 
quency updated once per interrupt. 

3. A mix of 1 and 2 is to ’’look ahead” at the trans- 
fer function where the sine w' 11 be sweeping to 


determine what sort of compensation will be re- 
quired. If the ti*ansfer function is relatively flat, 
sine control can be done less frequently than if 
the transfer function contains resor inces and/or 
nulls. 


PKRFORMANCE EXTENSIONS USING PARALLEL- 
PROCESSORS 

Since much of the sine control computation is inde- 
pendent of the random control computation, these two 
control calculations can proceed in parallel. The 
w’ide variety of low-cost microprocessors available 
makes multiprocessor systems economically viable. 
The sine control computations can be implemented 
with one or more microprocessors. 

For example, consider the block diagram of Fig. 11. 
The jjarallel microprocessor would perform (he pro- 
cessing functions enclosed by the dotted line. It 
would interrupt the main processor to get the trans- 
fer function, compute the control signal, compensate 
the drive, and update the sine synthesizer amplitude. 
This microprocessor would also update the sine syn- 
thesizer frequency using a real-time clock as a time 
reference. Memory requirements for the parallel 
processor would be minimal. One extra fxirallel pro- 
cessor should be able to handle four sine synthesiz- 
ers efficiently. 
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Another variation on this approach is to have the 
microprocessor generate the sine from a lookup 
table as well as compute the control algorithm. This 
sine data could be transferred to the main processor 
and added to the random data before output to the 
DAC, or the sine data could be output throu^ a sep- 
arate DAC and summed with the random data in the 
analog domain. Reasonable performance might be 
possible with a single microprocessor generating and 
controlling four sines. 

Use of a parallel processor makes control strategy 
number me viable. Since no random processing need 
be done, sine frequency update, control, and out,>ut 
are the only functions in the loop. This promises to 
provide the best control for reasonably fast, high- 
resolution sweep rates. The sweep resolution is. 
limited by the sine control loop time. This time iS 
greater if the microprocessor also generates the 
sine waveform as opposed to simply controlling a 
hardware sine synthesizer. Control strategies two 
and three can also be used if even faster sweep rates 
are desired. Note again that the transfer function 
estimate used in the sine control strategy is updated 
only once per random loop time. 


1 TEST IDi SR31 RCT 

2 HEADIMGt RECEIVING CHECKOUT TEST 

3 BANDVIOTHi 2048* 

4 RESOLUTION 64/126/2 56/ 512 t 128.0 
FREQUENCY INCR£MB«Ts 16.00 

REFERENCE SPECTRUHi 

5 INITIAL SLOPE# 03/OCTs 12*00 

ALARM LIMIT D3 t 6*900 

6 FREQUENCY HZ* t 16*00 

LEVEL GSQR/HZ.s 1*000 
ALARM LIMIT DB f 6*000 

7 FREQUENCY HZ*: 2000* 

LEVEL GSQR/HZ.s 1*000 
ALARM LIMIT DB t 6*000 

6 FINAL SLOPE# DB/OCTi *48*00 

RANDOM GRMSt 45*01 


9 SINE FR£QUENCY#HZt 500*0 
SINE REFERENCE LEVEL# 6*Si 5*000 


10 SINE FHEQUENCY#HZt 800*0 
SINE REFERENCE LEVEL# G*St 1*000 


IMPLEMENTATION OF METHOD I 

The techniques described in Method I have been imple- 
mented into a working digital swept sine on broadband 
random vibration control system. A typical setup 
procedure is shown in Fig. 12, and results of a test 
with four swept sines are shown in Fig. 13. This 
implementation has all the features and limitations 
described in Method I. 


SUMMARY OF PArV II 

Various techniques for implementing digital swept 
sine on broadband random vibration control systems 
have been presented. Several sine generation methods 
and sine oontrcl strategies have been proposed with 
indications of the expected relative cost- performance 
aspects of each. In particular, the use of parallel 
processors to improve system performance was dis- 
cussed. Feasibility of a digital swept sine on random 
vibration control system has been demonstrated by 
an actual implementation. 


11 SINE FR£QUENCY#aZ: 1000* 

SINE REFERENCE LEVEL# G*S: 10*00 


12 SINE FREQUENCY# HZt 1500* 

SINE REFERENCE LEVEL# G*St 5*000 

TOTAL GRMSs 45*84 

13 SWEEP PARAMETERSt 
MODE* :-LO6#0»LlNt 0 
RATE# H.VSECt 3*000 

INITIAL DIRECTION l»UP#0«DOVNt 1 

14 LOVER SWEEP LIMIT# HZt 400*0 

15 UPPER SYEEP LIMIT# HZt 1600* 

16 LOV LEVEL# -D3t *20.00 

17 LEVEL INCREMENT# DBt 2.000 

18 START-UP TIME SECt 1*875 

19 SHUT-DOWN TIME SECt *9377 

20 TEST TIME HRS# MIN# SECt 0#4#0 

21 CONTROL JHANNELSt 1 

22 ACCEL ScNS HV/Gt 10*00 

23 DRIVE CLIPPING l»Y£S#0-NOt 0 

24 ABORT LEVEL GRMSt 60*00 

25 ALARM LlNESt 

26 ABORT LINESt 

NOTE : Frequencies were automatically ad- 
justed by program to an integral 
multiple of the frequency increment. 


Fig. 12, Listing of Seum Parameters for Sine on Random 
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UHBPR LINOT 

<n> SZNE SKTRUH nFTER TBT <B> fmXVI CONTROL SPECTRUH 


Fig. 13. Test Results with Four Swept Sines 


REFERENCES 

1. Heizman, C. L. , ’’A High Perform, nee 
Digital Vibration Concrol and Analysis 
System, " Institute of Environmental 
Sciences, 1972 Proceetlings , pp. 309-315. 

2. Nelson, D. B. , ' Performance and 
Methodol Dg>’ of a Digital Random Vibration 
Control System,*' Ii^stitute of Environmental 
Sciences, 1973 Proceedings , pp. 187-192. 

3. Norin, R. S. , "Pseudo Random and Random 
Testing, " Seminar on Understanding Digital 
Control and Analysis in Vibration Test 
Systems, 1975 Proceedings. 


4. Barthmaier, J. P. , "Shock Testing Under 
Minicomputer Contrcl," Institute of 
Environmental Sciences, 1974 Proceedings. 

5. Darthmaier, J. P. , 'Two Approaches to 
Computer- Controlled Shock Synthesis, " 

The 20th International Instrumentation 
Symposium, 1974 Proceedings. 

6. Cooper, H. W. , ”\Vhy Complicate Frequency 
Synthesis Electronic Design . July 19, 1974, 
pp. 80-84. 

7. Norin, R. S. , "A Multi-Channel, Multi- 
Strategy Sinusoidal \Tbration Control 
System, *' Institute of Environmental 
Sciences, 1974 Proceedings. 



















A REVIEW OP 

ENVIRONMENTAL TEST INNOVATIONS 
PERMITTED BY DIGITAL CONTROL SYSTEMS 


W. Brian Keegan 

NASA-Goddatd Space Plight Center 
Greenbeltf Maryland 


This paper reviews the innovations of digital control 
systems from the perspective of one who establishes 
enviroranenta" tesx. requirements for aerospace hardware 
to the dyna-^Lc environments, of vibration and mechanical 
shock. It eviews the^^pe of options that specifica- 
tions might permit in those laboratories that possess 
digital control systems. The innovative tests mentioned 
consider only the capability of digitaiL control systems 
and are not presented as formally approved chaj\ges to 
recommended test philosophy. They are intended simply 
to highlight the kinds of tests that# in the future, 
might be considered to constitute acceptable portions of 
the environmental test sequence of aerospace hardware. 

Each of the environments of sinusoidal vibration, random 
vibration, and mechanical shock are addressed separately. 
In each area, potentially innovative test methods are 
discussed. The paper concludes with a warning to avoid 
the temptation to perform unnecessarily complex environ- 
mental tests simply because the capabilities of digital 
control permit them to be run. 


INTRODUCTION 

With the advent of digital control 
systems for performing environmental 
tests, sim ilation techniques previously 
impossible or at best tediously diffi- 
cult with analog control systems have 
been made far mor- straightforward. 
Digital systrms thus facilitate the 
use of more accurate simulation methods, 
enhancing the realism of tha environ- 
m.^nts to which hardware is subjected 
prior to acceptance for se.rvice. This 
paper views the advantages of using 
digital control not from the perspec- 
tive of the test equipment designer 
or operator, but rather from that of 
an envirorjnentalist responsible for 
establishing test criteria for aero- 
space equipment which must survive the 
rigors of vibration, acoustic noise, 
and mechanical shock upon Initiation of 
its service life. Simultaneously, it 
views the audience as being astutely 
aware of environmental facility design 
and operational capabilities, but not 
necessarily well grounded in environ- 


mental test philosopli' for aerospace 
equipment. 

When originally invited to ^resent 
this paper, the suggested topic ^as 
"Test Specification Changes Required 
for Environmental Test Implementation 
Using Digital Test Systems.” Upon 
consideration of the kind of material 
that could be presented, however, the 
word "required" seemed to be unneces- 
sarily restrictive, thus sparking the 
change to the current title. In fact, 
no changes are necessitated by the 
change to digital control systems with 
the exception of some minor ones to 
test tolerances, as will be discussed 
later in this paper. Since analog 
control systems pre-existed digital 
ones by several years, the original 
design of digital systems of necessity 
had to, possess all the capabilities of 
analog systems and then some in order 
to make them marketable items. Addi- 
tionally, from the standpoint of 
maintaining specifications that are 
universally applicable to all 
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environmental test laboratories# it - 
will be necessary for some time -to' 
come to consider those iabs; that still 
utilise analog control systems. Most 
tests on aerospace hardwere are com- 
ponent level tests wtich are often 
performed by smaller labs that can 
logically be expected to be among the 
latter companies purchas:^ng digital 
control systems. 

Wilting this paper was viewed more 
as an opportunity to describe how one 
might take advantage of, the increased 
capabilities offered by digital con- 
trol equipment. Thus the paper re- 
views the type of options that might 
oe permitted for those laboratories 
that possess digital control systems 
in order to perform more realistic 
simulations of the environments of 
sinusoidal and random vibration and 
mechanical shock. For each of these 
three ©reas the problems that may be 
encountered with existing test toler- 
ances brought about by the change to 
digital control will be discussed. 
Additionally# the increased test 
complexity facilitated by# as well as 
the innovative typos, of tescs oermitted 
by digital systems# will be reviewed. 

As such# the ensuing discussions are 
intended to highlight the kinus of 
tests that might in the future con- 
stitute acceptable portions of the 
environmental test sequence of aero- 
space hardware, 

SIIJLSOIDAL VIBRATION 

As implied in the previous section, 
it is not the purpose of this '^er to 
address the problem of whether ^ ■ nu- 
soidal vibration testing should be 
performed. }?ather# it stai:u3 ^Jth the 
ground rule that sine testing must be 
performed and considers how it might 
be done better using digital contiol 
equipr*.-*nt. 

First# consider a typical sinu- 
soidal test profile such as the one 
illustrated by the solid line in Figure 
1# -hich has been extrac " -id fiom l:he 
NASA- Goddard test spec*-’ -ation for 
spac icrafL. Such a i': of constant 

acc i iraoion segments a si:aplirl:ic 
env* lope of the illust- -.ced flight 
data that has been analyzed by either 
narrov:-bancI analysis c • shock response 
spectrum teenniques. in almost all 
instances such flight aata could have 
been more closely enveloped had capa- 
bility existed in the -est equipment 
to easily run non-com.tant acceleration 
versus frequency prof* les. It i3 known 
that capability has 1 .ng existed in 
analog systems for co '.trolling at the 


fixed positive slope values of +3 da 
per octave and dB per octave by 
controlling respectively to constant 
velocity or constant displacement butL. 
these arc limited in frequency since 
the once or twice differentiated 
acceleration feedback signal rapidly 
becomes too small to accurately con- 
t ol. Moreover# it is known that 
analog curve followers exist that per- 
mit test profiles to contain any 
acceleration slope value at any fre- 
quency. These# however# are bulky 
peripheral units that are time consum- 
ing and tedious to set up and use. 

Digital systems on the other hand 
are programmed such that the accelera- 
tion ver^BUS frequency profile can be 
subdivided into any reasonable number 
of straight-line segments simply by 
specifying the breakpoints. This 
drastically simplifies the process of 
running sinusoidal tests whose profiles 
are far more realistic such as that 
envisioned by the heavy dashed line 
in Figure 1, 

There do exist times when a step 
change in acceleration is the desired 
way to define c rlnusold^l test spe- 
cification. One example of such an 
instance would be the simulation of 
the Pogo event of the Delta Launch 
Vehicle. Pogo# briefly described, is 
an oscillation of the launch vehicle 
and payload at its first longitudinal, 
”accordian*' , resonant mode that is 
induced near the end of ‘ f irst-scage 
burn by coupling between the launch 
vehicle primary structure and pressure 
oscillations in the propellant feed 
line. The amplitude of this oscilla- 
tion has been measured on enough 
flights that it is amenable to statis- 
tical analysis. The frequency of 
occurrence is predictable to within a 
frequt 'ey band equal to the change in 
the first mode frequency over a few 
seconds flight time. As a result# 
there exists a narrow frequency band 
over which a relatively high sinusoidal 
acceleration test level is required. 

At frequencies outside this band, the 
flight data snows significantly lower 
levels and the data envelope is gener- 
ated by events that occur at othar 
flight times. Thus# a step increase 
and a step decrease, as illustrated in 
the lower f requenc:^ portion of r igure 
I, are a desired portion of the test 
profile. 

In such instances, digital control 
systems may not be able zo complete 
the step change ac cel ere* Mon within 
the traditional frequency tolerances 
of two percent. The customary procedure 
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for performing such a change with an 
analog system is to terminate che 
sweep, change the input vihrrrtion 
level, and when control has stabilized 
at the new level, le-initiate the 
sweep. This process undoubtedly ax- 
poses the test item to more high- 
amplitude stress reversals than is 
realistically required. It is accept- 
ed practice though because it is 5elt 
to be the optimum trade-off between 
maintaining control and completing 
the level change within an acceptable 
frequency band. Digital systems, 
being more highly automated, do not 
stop sweeping but rather change level 
as they sweep. Since control systems 
have finite compressor speeds, they 
require a finite length of time to 
change levels without losing control 
of the process. The bandwidth in 
which this stabilization to the new 
level can be completed is a function 
of the percentage change in level, 
the compressor speed, the sweep rate, 
and the frequency at which the change 
was initiated, as is defined in Ref. 1. 
There it is stated that: 



where 

Af = the bandwidth needed to accomplish 
the level change 

f = the frequency at which the level 
change is initiated 

gi/g 2 the amplitude ratio of the 
^tep change 

N = the sweep rate, and 

S = the compressor speed of the con- 
trol system at f 

The approximation error in the above 
equation is less than one percent. 

As a specific example of the ibove, 
consider the Delta Pogo specif icatior 
where at 17 Hz there is a step change 
from 1.5g to 4.5g while sweeping at 
4 octaves per minute. Using, from 
Ref. 1, an average compressor speed 
of 4 dB per se:;ond for a digital 
control system dt 17 Hz, it can be 
seen that it would take a 2 Hz band to 
complete the step change, ixiis is con- 
siderably greater than the currently 
acceptable 2 percent frequency band in 
which such step changes in level must 
be accomplished, 

Thu digital control algorithms 
must, of course, be ./ritten vo run 


safe tests and the permissible" test 
tolerances must be broadened to permit, 
operation ^f the test equipment that 
does not jeopardize the integrity of 
the test item. Of more potential sig- 
nificance, however, is, the fact that 
the test engineer must be made aware 
of how the control algorithm is 
written so that he can anticipate the 
response of- the item under test^ 
Suppose, for exan»ple, a spacecraft had 
a fundamental resonance just bclnw or 
just above the frequency band of 
Figure 1 in whic^ the input levels are 
the highest. If, in order to meet the 
inpx’u amplitude requirements in that 
band, the test level began to change 
well below the beginning of the band 
-ad was not fully reduced until well 
after the band, the possibility would 
exist of exciting the structure at its 
resonant mode to tinrealist ically high 
levels, significantly high enough such 
that special provisions should have 
been made for limiting the response of 
the structure. Thus, it can be seen 
that the problem goes well beyond the 
operational one of meeting specifica- 
tion tolerances and becomes one of 
assuring that percinent information 
regarding system perf ormar ::e capa- 
bilities has been disseminated to the 
proper individuals who mu >t utilize 
it for proper test planning. 

Thi^' brings us logically then to 
pernaps the most significant advance- 
ment in using digital control systems 
to f'erform sinusoidal vibration tests: 
the ability to automatically limit the 
response of many points to some pre- 
determined response value, thereby 
precluding inadvertent cvertests due 
to unanticipated high amplifications. 
For several years .i-.w, analog systems 
have had limited capability to trans- 
fer control to some response point 
wheJi a predetermined limit was reached 
or to initiate a controlled shutdown 
if some abort limit was exceeded. 

These teatao.es are mandatory on 
systems that are required to contro] 
tua various environments to which 
today's ccmple::, multi-million dollar 
space^ra^c are tested. In analog 
■'stems, ♦‘hese features were cor tained 
in peripheral hardware that took up 
precious srece in the control room and 
was often (iiificalt and time consumir;g 
to check cut, '"^librate and set up. 

More important!'', limitations on equip- 
ment availabi Ix-ty usually placed severe 
constraint on the number of respv^r.'^o 
points that could baI*SO monitored, 
often necoss it a. ,■* nrj a cr'mpromise be- 
tween the level of sophi.>ticc tiCi. oriv- 
desired to have in a test ^ncl the level 
one was .able to have. 
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With digital control systems, how- 
ever, such limiting and control of 
responses other than the primary con- 
trol point is contained in the control 
algorithm where it requires no addi- 
tional hardware space and can be set 
easily through the .^letype inter- 
face* Additionally, such capability 
exists in great amounts allowing 
many response points to be limited, 
thereby permitting far more complex 
tests to be run than are currently 
possible in laboratories with analog 
control equipment. A word of caution 
may be in order, however, in that one 
should not r- ' unnecessarily complex 
ests siip^ • / 'ecause the control sys- 
. m has un^‘*^^d capability. Sound 
erqT neering judg«nent must still 
r mpcr the complexity of environinental 
\.sists since the likelihood of a setup 
error increases with test complexity 
no matter how simple the cntrol 
system may be to use. 

The final point of discussion re- 
garding sinus cidal testing concerns 
the use of variable sweep rates during 
the vibration profile in the control 
algorithms of some manufacturers of 
digital equi-tmient. This feature re- 
duces the sweep rate at frequencies 
at which the system transfer function 
is rapidly changing in order to more 
accurately control the applied vibra- 
tion level to that desired. This 
approach is felt to have some distinct 
disadvantages and the potential pro- 
blems associated with its logic will be 
drawn out in the following example 
cited from a test on the Orbiting 
Astionomical Observatory (OAO) at the 
Goddard Space Flight Center. 

Using an analog control system and 
a fixed sweep rate it was desired to 
perform a base excited sinusoidal 
vibration test of the 5000 pound Space- 
craft at a constant acceleration level 
of 0.25g while sweeping at a rate of 
4 octaves per minute from 5 to 15 Hz. 
The actual control acceleration pro- 
file measured at the base of the space- 
craft was as shown in Figure 2a. As 
can be readily seen there, the input 
varied by a considerable amount from 
the desiied in the vicinity of the 9 
Hertz first cantilevered bending mode 
of the spacecraft. The explanation for 
this occurrence is that as one ap- 
proaches the fixed-base resonant fre- 
quency of such a large test item, the 
apparent mass, and hence the driving 
force required to be supplied by the 
shaker armature to maintain a constant 
acceleration level, changes more 
rapidly than the compressor of the 
control system can compensate for it 


(See Figure 2b) . Thus, the input 
acceleration falls below the desired 
level and continues to fall further 
below until the slope of the apparent 
mass decreases to the point that the 
compressor can begin to compensate and 
adjust the input to the desired level. 
After the resonant frequency is passed, 
the inverse action occurs in that as 
the required force decreases more 
rapidly than the compressor can react, 
the input rises above the desired 
level until such time as the slope of 
the apparent mass levels off and \e 
compressor is again able to prov^. e 
acceptable control. 

Although the above problem was 
described in the frequency domain, it 
is in fact a problem of the time domain, 
in tliat the main cause is the relative- 
ly slow compressor speed. Increasing 
the compressor speed is an unacceptable 
solution to the problem since it may 
cause servo loop instabilities thereby 
losing control completely as it attempts 
to correct too rapidly for changes in 
the system transfer function. De- 
creasing the sweep rate, however, 
effectively increases the compressor 
speed without creating potential in- 
stabilities by reducing the time rate 
of change of the system transfer 
function thereby allowing the compres- 
sor more time to respond and enabling 
it to significantly lessen the devia- 
tions of the actual vibration level 
from the desired. 

From a test operations viewpoint 
then, this automatic reduction of the 
sweep rate below the specified level 
at frequencies at which the system 
transfer function is changing rapidly 
is an ideal solution. It permits 
running more precise tests and is most 
likely the reason why some manufac- 
turers have included it in their 
sinusoidal vibration control algorithm. 
There are, however, considerations 
that should be made when performing 
tei»ts over and above that of maintain- 
ing precise control of the input level. 

The prime consideration during vibra- 
tion testing is that the proper load be 
applied for the desired duration. 
Variations from the desired level as 
severe as those illustrated in Figure 
2 occur only at the fundamental canti- 
levered resonances of very massive test 
items. It is at these frequencies that 
the Goddard test philosophy permits a 
reduction of the input vibration level 
below the specification value so as to 
limit the response of the test article 
to the maximum predicted by analysis 
to occur in flight. Thus, while in the 
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cit>^d axarople the input acceleration 
was tne control parameter, it was not 
a specified acceleration level but 
rather some nominal acceleration pre- 
dicted by the results of a low-level 
survey to be the one needed to develop 
the desired response load in the 
critical areas of the spacecraft 
structure. So long as the percentage 
deviation of the control acceleration 
from the desired was the same for both 
the full level run and the low level 
run that had been used as the basis 
for the extrapolation, then the full 
level run would develop the desired 
responses. By slowing down the sweep 
rate in the vicinity of test item 
resonances, one would still need to 
employ some technique that limited the 
responses of the test to those desired. 
Thus, the same peak load would be 
developed. But by slowing down the 
sweep rate, the number of cycles for 
which the peak load was applied could 
increase drastically. When sweeping 
at the constant rate of 4 octaves per 
minute, the OAO Spacecraft saw a load 
greater than 90 percent of maximum 
for 19 cycles in the vicinity of the 
9 Hertz resonance. This was already 
considerably in excess of the number 
of cycles for which the load was pre- 
dicted to last by the flight dynamics 
loads analysis. To significantly in- 
crease this number of applied load 
cycles even further in order to provide 
improved control accuracy by slowing 
down the sv/eep rate seems unreasonable. 

The fact that fatigue damage is not 
usually a problem with aerospace hard- 
ware is not an acceptable argument 
since one of the basic tenets of envi- 
ronmental testing as envisioned by the 
Goddard Space Flight Center is to 
expose hardware to as realistic a 
simulation of the service environment 
as is practicable. It is felt that 
the reduced sweep rate feature loses 
sight of this objective and for this 
reason it is recommended that it be 
excluded from the sinusoidal vibration 
control algorithm of digital control 
equipment. 

RANDOM VIBRATION 

For the control of random vibration 
inputs, digital control systems offer 
significant improvements over analog 
systems for reasons of both shorter 
equalization time and finer spectral 
resolution. Using a digital system, 
the bandwidth of the vibration profile 
can be subdivided into 400 or more 
filters. Assuming 400 filters are 
used, then a profile over the tradi- 
tional bandwidth of JO to 2000 Hertz 


could be equalized with an approxi- 
mately 5 Hertz resolution. This 
compares with existing analog systems 
that have 80 filters varying from 12 
to 50 Hertz. 

It is this improved resolution that 
is r.elt to be a potential problem area 
not from the standpoint of being not 
as good a test but from the standpoint 
of meeting existing specification 
tolerances. These traditional toler- 
ances require that in any single 
filter bandwidth the power spectral 
density (PSD) be controlled to within 
3 dB of the desired value. This is 
not a serious problem with wider 
filters because it is really the aver- 
age PSD value in the band that is 
being controlled within the specified 
tolerances. As the filters become 
very narrow, the average PSD value in 
a single filter bandwidth is radically 
affected by narrow band, high ampli- 
fication resonances in the test item. 
Although no actual data was available 
to the author, the case can be en- 
visioned in which the digital control 
system was unable to compensate 
sufficiently for such a resonance be- 
cause of its finite dynamic range. If 
such cases were encountered regularly, 
then specification tolerances would 
have to be broadened to permit greater 
than 3 dB variations from the desired 
level over a certain percentage of the 
narrow band filters. It must be empha- 
sized, however, that this is not a 
case of digital systems being unable 
to control as well as existing analog 
systems. It is rather a case of digital 
control systems being sensitive enough 
Lo detect far more sophisticated pro- 
blems than analog systems can and thus 
requiring a change to test specification 
tolerances in order to take into ac- 
count this increased level of sophis- 
tication. 

As far as random vibration control 
techniques are concerned, there is 
one area which the control algorithms 
do not currently address that is 
particulary amenable to the capabilities 
cf a digital system. That is the abil- 
ity to control vibration responses 
generated in one frequency band by 
energy input in another frequency band. 
Or.e example of such a case would be a 
resonance excited by input energy in 
th^ skirts of the filters adjacent to 
01 . in which the resonance lies. 

Another would be the excitation of re- 
sponse at one frequency generated by 
energy input at eigher a harmonic or 
sub-harmonic of that frequency. To 
the author’s knowledge, existing con- 
trol techniques utilize the measured 
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response of le control point in the 
filter band to regulate the input 
force in only that same band so that 
the measured response aqua’s that de- 
sired. Tn the cases mentioned above, 
such a cc^ntrol algorithm would have 
no eflec on the equalization of the 
PSD a . Cl e control point. 

When such a case is encountered 
using eu: ilog control systems, the 
operator, begins a process of manual 
equalization by adjusting the various 
sliders to see if equalization can be 
accomplished. This can be a time 
consuming process and often results in 
a decision to run the test with one 
filter knowingly outside the tolerance 
limit simply because one does not wish 
to take the extensive time required to 
adjust each filter until proper equal- 
ization is achieved. 

Although such instances may not be 
encountered frequently, it would seem 
worthwhile to have provisions for them 
incorporated into the random vibration 
control algorithms of digital systems. 
Such provisions would consist of se- 
quentially changing the Input energy 
in each filter band and checking to 
see if the response in the unequalized 
band was affected. In such a manner 
the filter generating the response 
could be located and a secondary con- 
trol loop established whereby its 
output energy was controlled by the 
response in the unequalized band. 

When such an equalization problem was 
encountered, this supplemental portion 
of the control algorithm could be 
activated by a simple teletype input 
in order to provide capability of con- 
ducting a good environmental test even 
when the test item possesses such 
unusual dynamic characteristics. 

The final topic for discuf.sion 
under this heading of random vibration 
concerns extension of the existing con- 
trol algorithms to a new area. Current 
practice calls for random vibration 
tests to be run for same specified rime 
at Some overall RMS acceleration Jevel. 
When one speculates as to how they 
might improve the realism of random 
simulation, the idea of a time varying 
RMS acceleration level, similar to the 
profile illustrated in Figure 3, comes 
to mind. This profile shows the RMS 
acceleration level measured ac the 
interface of a spacecraft with the 
Delta Launch Vehicle from lift-off 
through the regions of transonic 
flight and maximum aerodynamic pressure. 
While no firm test requirements for 
such a profile yet exist, it seems 
reasonable to think that such a profile 


may some day be a requirement. More- 
over, it seems that implementation of 
such a requirement could be easily 
added to the random vibration control 
algorithm of digital systems. When 
room for potential improvenents exists, 
requirements have the demonstrated 
ability of growing to meet such poten- 
tial. The writer feels that such will 
be the case for the time varying RMS 
random vibration profile discussed 
here. 

MECHANICAL SHOCK 

The area of mechanical shock is the 
one in which digital control systems 
have provided the greatest improvement. 
They offer two distinct methods of 
environmental simulation: shock re- 

sponse spectrum matching, in which the 
potential simulation accuracy has been 
significantly improved over that of 
analog systems through increased re- 
solution, and transient waveform 
control, a completely new test capa- 
bility that is not environmental 
simulation at all but rather environ- 
mental duplication. Prior to a 
discussion about either of these 
methods, let us depict the environment 
that we are trying to simulate. De- 
picted in Figure 4a is a shock tran- 
sient measured on the primary structure 
of a typical spacecraft during a 
pyrotechnic event such as separation 
from the launch vehicle. It is a 
high-G multi-frequency oscillation 
characteristic of such events. Figure 
4b meanwhile presents the acceleration 
sho*'k response spectrum of the tran- 
sient for a critical dampina ratio of 
0.05. Both of these are representa- 
tions of the same environment, one 
being a description of structural 
motion and the other a prediction of 
the amount that single degree of 
freedom systems with the defined 
damping would respond to that structur- 
al motion. If one wished to simulate 
the event, then he has the freedom to 
work in either the time domain or the 
frequency domain. 

Recent pas:t practice with analog 
systems has been to simulate the de- 
sired shock spectrum by synthesizing a 
time history whose shock spectrum 
matched, for the specified damping 
value and within the allowable toler- 
ances, the desired shock spectrum. 

This was done by summing several 
discrete frequency wavelets into a 
single complex time history with the 
amplitude of each frequency component 
individually adjusted to control the 
resultant shock spectrum at its own 
frequency. The frequency spacing, or 
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resolution, of these frequency compo- 
nents was usually one-third octave, but 
it must be understood that such shock 
spectrum synthesis techniques do not 
utilize filters that control the spec- 
trum over the entire one- third octave 
band* True control of che shock 
spectrum magnitude exists only at the 
one-third octave center frequencies. 
Because of the shape of the shock 
spectrum of each wavelet the shock 
spectrum magnitude fal3s below the 
desired value at all frequencies between 
adjacent one- third octave center fre- 
quencies* As a result, for test items 
whose predominant response frequencies 
lie between one-third octave center 
frequencies, an undertest may result. 

The concep;: of shock spectrum 
synthesis is the same for both analog 
and digital systems. The use of digital 
systems does offer improved simulation, 
however, by providing increased resolu- 
tion, su h as one-sixth or one- twelfth 
octave, thereby providing exact match- 
ing of the desired spectrum magnitude 
at a greater number of frequencies that 
are also closer together. Thereby, the 
amount of inaccuracy in spectral values 
at frequencies between two adjacent 
control frequencies is also signifi- 
cantly reduced. 

One serious deficiency of any 
method of shock spectrum matching is 
that it is good for only one value of 
damping. If the damping coefficient 
of the particular item under test does 
not equal that for which the shock 
spectrum was matched or if the item is 
sufficiently complex that its behavior 
is not at all predictable by the single- 
degree-of- freedom analogy upon which the 
shock spectrum concept is based, then 
the response induced in the test item 
by the test environment will not be the 
same as that induced by the service 
environment and the test will not have 
accomplished its objective. This is so 
for a given time history because the 
variation in shock spectrijm magnitude as 
a function of the analysis damping coef- 
ficient is dependent upon the number of 
cycles for which the time history lasts. 
This is illustrated in Figure 5 where the 
variation of peak shock spectrum magni- 
tude with damping has been plotted for 
several types of time history. As can 
be seen, if the test item has a lower 
damping value than that for which the 
shock spectrum was computed, as is often 
the case because many subsystems have 
high-Q resonances and most specifications 
are defined for Q equal to 10, and a 
typical shock spectrum synthesis technique 
is employed, then the simulated environ- 
ment will produce an overtest. 


One method of controlling the 
amount of over test or under test so in- 
duced is to control the number of 
cycles in each individual wavelet of 
the simulated environment such that 
the variation in spectral magnitude as 
a function of damping approximates that 
seem in the real environment. It is 
understood that certain digital shock 
spectrum synthesis control algorithms 
possess such capability, and it would 
seem advisable for test laboratories 
to have such capability if possible. 
This does not, however, provide a 
universal solution because the ratio 
between peak shock spectrum magnitude 
and peak magnitude of the input time 
history is reduced as the number of 
cycles in each wavelet is reduced 
thereby lowering the upper magnitude 
of ^ -e shock spectrum that can be so 
synthesized due to shaker and amplifier 
force limitations. 

One final shortcoming of the shock 
spectrum matching technique is that, 
even if excited to the proper response 
levels, the test item experiences peak 
stresses for a greater number of cycles 
because the simulated environment 
usually lasts longer than the service 
environment. The way of avoiding the 
problems of spectral matching complete- 
ly is to utilize the transient wave- 
form control capability of digital 
systems. With this technique, accurate 
duplication of the measured service 
environment can be achieved. Evalu- 
ating the responses induced in the 
test item by this method in terms of 
the foregoing shock spectrum matching 
discussion, then, there are no fre- 
quencies at which the spectrum magni- 
tude of the test environment varies 
from the desired spectrum value. In 
addition, the variation in shock spec- 
trum magnitude as a function of the 
analysis damping coefficient is pre- 
cisely the same for both the test 
environment and the service environ- 
ment. Thus the necessity of knowing 
the damping coefficient of the test 
item in advance in order to assure a 
realistic test is eliminated. Finally, 
the number of peak stress cycles is 
precisely duplicated since the duration 
of both the test and service environ- 
ments are exactly the same. Transient 
waveform control therefore is a signi- 
ficantly improved shock test method 
offered by digital control systems that 
should be adopted whenever possible. 

Prior to implementation of such a 
test method for other than development 
testing, however, an approved specifi- 
cation must exist and a method for 
generating a specification that provides 
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a qualification time history is not 
straightforward. In the case of devel- 
oping a shock spectrxim specification, 
one needs to do a shock spectrum analy- 
sis tor some defined damping coefficient 
over the frequency range of interest 
of all flight data deemed appropriate 
for inclusion into the particular 
specification being developed. A test 
specific?: cion can then be constructed 
by either enveloping the composite 
spectra encountered in service or if 
sufficient data samples exist, by per- 
forming a statistical analysis to 
determine some desired environmental 
probability level of occurrence. Com- 
bination of the various data samples 
is thus easily achieved when developing 
a shock spectrum specification. 

Still in terms of the shock spec- 
trum, let us now consider the develop- 
ment of a qualification time history. 
Seldom if ever does a single time 
history generate the worst case shock 
spectrum at all frequencies. Thus, 
there will probably never exist a single 
time history that can be uniquely de- 
fined as the worst case for a particular 
specification. Moreover, because the 
shock spectrum is not a reversible pro- 
cess, one cannot take a composite 
envelope shock spectrum of several time 
histories and generate a unique time 
history for definition in a test spe- 
cification. The most likely time- 
history specification that would be 
developed then would be a set of such 
time histories, each of which would 
generate the worst case test item re- 
sponses over a portion of the frequency 
range for which the specification was 
defined. Such specif ications should 
begin to be defined in the near future. 

One of the expected advantages of 
performing mechanical shock tests 
under the control of digital systems is 
the ability to improve the accuracy to 
which the desired responses can be 
controlled. Current tolerances on 
shock spectrum tests are +50%, -10%, a 
large band in which test amplitudes may 
fall and still be termed acceptable. 

The usually quoted accuracies on tran- 
sient waveform control are il0%. This, 
however, assumes that no non-linearities 
are encountered between the equalization 
level and the full level test. When 
equalization is conducted 6 dB or more 
below full level, as is usually re- 
commended, non-linearities well in 
excess of 10% can be expected, as is 
illustrated in Ref. 3. When this 
occurs, even perfect equalization will 
not result in tests that meet the 
stat-:d capabilities of ±10%. 


The writer feels that significantly 
improved mechanical shock test accuracy 
could be achieved if both the shock 
spectrum and transient waveform digital 
control algorithms were expanded to 
include measurement of and compensation 
for the non-linearities inherent in the 
test item. One suggested method for 
accomplishing this would be to first 
equalize initially at the -12 dB level 
(25%) , and next perform a run at the 
-6 dB level (50%) . The differences 
between the predicted inputs for the 
-6 dB run, based on a linear extra- 
polation of the - 12 dB run, and the 
actual inputs could serve as a measure 
of the non-linearities encountered in 
the transfer function of the test item 
when subjected to a 6 dB increase in 
excitation level. This non-linearity 
measurement could then serve not on3.y 
to re-equalize at the -6 dB level, but 
more importan 3y to anticipate the 
additional non-linearities now expected 
to be encountered between the -6 dB 
and the 0 dB (100%) runs. Thus, there 
would exist a technique for at least 
attempting to compensate for test item 
non-linearities. The Goddard Shock 
Test Handbook (Ref. 4) currently re- 
commends such a procedure be followed 
for manually equalizing shock spectrum 
tests using the analog synthesizer. 
Incorporation of this technique into 
the digital control algorithms would 
improve an already good test method. 

Even with the vast improvements in 
shaker shock test methods facilitated 
by digital control systems, there still 
exists room for argioment that any form 
of on - the - shaker test induces 
much greater responses throughout the 
test item than are seen if one exposes 
the same test item to the actual pyro- 
technic event. One possible explanation 
of this is that on the shaker the input 
is coherent at all mounting points of 
the test item, whereas if one measures 
the time history at the various mount- 
ing points of a test item during an 
actual pyrotechnic event, he will see 
that they are not exactly the same. 

This topic, however, does not warrant 
further discussion in this paper, but 
is undoubtly one that will receive much 
attention in the near future. 

SUMMATION 

Thus, the capabilities of digital 
control systems have been reviewed as 
they interface with environmental test 
specifications. They provide vastly 
increased capabilities in some areas, 
and in others offer test methods pre- 
viously impossible. They do not, how- 
ever, eliminate the need for constant 
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care on the part of the test engineer 
to verify that the objective of realis- 
tic environmental simulation is being 
met. 

Finally, then, it is felt that a 
warning must be sounded. Avoid the 
temptation to perform unnecessarily 
complex environmental tests simply be- 
cause the extensive capabilities of 
digital control systems to autcxna tic ally 
monitor large numbers of parameters 
permit -hem to be run. Sound engineer- 
ing judgment must continue to be the 
hallmark that identifies cost effective 
test programs that have properly bal- 
anced te&c complexity against the 
ramificationf. of inducing an unrealistic 
failure because of an improperly con- 
ducted test. 
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Figure 5: Variation of Peak Shock Spectrum Magnitude 

with Damping for Various Waveforms 
(partially extracted from Ref. 2) 
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SAFETY PROTECTION OF TEST ARTICLES 
USING DIGITAL CONTROL SYSTEMS 


R* A. Dorian 

NASA, Goddard Space Flight Center 
Gr^enbelt^ Maryland 


Safety protection of test articles is a primary concern in 
aerospace vibration testing. Digital control systems have 
greatly enhanced the protection of test articles because of 
the speed and repetition with which they can perform mono- 
tonous tasks. Because of this, the greatest source of 
error, operator error, has been extensively reduced. Alarm 
and abort tolerances of individual spectral lines are moni- 
tored and compared to the operator's predefined values, and 
an orderly shutdown of the system is initiated if the alarm 
and abort strategy is violated. The systems also initiate 
corrective action, if required, by monitoring primary power 
source levels, control signal integrity and peripheral equip- 
ment environmental requirements. However, "exciter dumps," 
which are due to unwanted transients generated at the output 
stages of the control system or the power amplifier, are 
corrected by a separate system. In this case, a system ex- 
ternal to the control system must monitor displacement, 
velocity, or acceleration of the exciter system and generate 
a controlled deceleration shutdown which prevents damage to 
the test article caused by the transients. 


INTRODUCTION 

One of the most undesirable 
occurrences during a vibration test 
is a shutdown due to the firing of 
an armature protector because of 
operator error, equipment malfunction, 
or power outage. 

Present armature protectors are 
designed to protect the moving element 
(armature) of the shaker by preventing 
it from hitting the mechanical stops. 
Unfortunately, the sudden dynamic 
braking generates large accelerations 
which have caused the spacecraft or 
test article to be subjected to ex- 
cessivi levels and possible damage. 

The firing of an armature protector 
or "exciter dump" is particularly 
important since in most cases these 
exciter dumps are due to operator 
error rather than equipment malfunc- 
tion or power outage. 

To this end equipment manufactur- 
er's have dealt effectively iii sol^Kng 
"specimen protection" problems. The 
computer has been especially helpful 


since its speed and versatility allow 
it to perform a number of operations 
such as: 

I. Checking operator error during 
the test set-up phase. 

II. Monitoring test levels through 
various control strategies and 
warning the operator of impend- 
ing problems. 

III. Housekeeping operations - moni- 
toring pertinent peripheral 
equipment parameters during a 
test. 

It would seem logical to draw the 
conclusion that with all of the above 
capability a digital control system 
would circumvent any potential "exciter 
dumps." However, this is not necessar- 
ily true. A transient occurring within 
the power amplifier circuitry or at the 
output stages of the control system 
could not necessarily be controlled by 
the computer and therefore would pass 
to rhe exciter and test article. There 
are protective circuits, however, which 
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can decelerate the exciter in an order- 
ly fashion without harming the test 
article* 

DETAILED DISCUSSION 

I DIGITAL CONTROLLER PROTECTION 

A. Checking Operator Errors During 

Test Setup Phase 

All digital control systems uti- 
lize what is commonly called a 
"conversational language" between the 
operator and the computer. During 
this conversation, the communication 
between the two takes place via tele- 
typewriter or a CRT/keyborad opera- 
tion. The operator responds to the 
question by typing an answer which is 
examined by the computer and either 
accepted if correct, or restated, if 
incorrect. Many tedious operations, 
which are handled by the computer, are 
certainly monotonous and consequenty 
prone to error, if performed by the 
operator. For example, the computer 
checks gain switches, performs cali- 
bration calculations, reference ampli- 
tude spectrum, values, error specriuiu, 
calculations, etc. and at a much 
faster rate than an operator. 

A typical test set-up "conversa- 
tional routine" is shown in Table 1, 
"Random Test Setup." As one can see, 
a number of parameters must be satis- 
fied in order to ensure a successful 
test. 

The operator specifies such items 
as transducer sensitivities, G^/Hz at 
a specific frequency, ov€»rall G rms 
level, etc. The system automatically 
controls the sic^nal conditioners, ana- 
log to digital converter, switches, 
etc., thus relieving the operator of 
these tasks. After the setup is com- 
pleted and before any full level vibra- 
tion test proceeds, the system provides 
a self-checking operation of the entire 
’dbration system, including power 
arplifier, exciter and transducer 
channels, to ensure proper system con- 
nections, calibrations, and switch 
settings. As an example, the system 
may discover that the transducer sen- 
sitivity is too high for the test 
sp'jctrum requested, in which case, the 
computer notifies the operator of such 
a condition. 

B. Monitoring Test Levels Through 

Various Control Strategies And 

Warning the Operator of Impending 

Problems 

Besides the self-checking procedure 


used by digital control systems prior 
to the test, these systems use various 
other techniques to protect the test 
specimen. 

1. Control Strategies 

Digital control can be accomplished 
through various techniques in order to 
protect the test specimen as much as 
possible. The following is an example 
of control techniques: 

a. Control on the average of multi- 
ple selected channels. 

b. Control against specified limits 
on specified signals. 

c. Control on the highest or lowest 
of selected signals. (sine sweep) 

d. Control on filtered or unfilter- 
ed (peak or rms) signals, (sine 
sweep) 

The operator selects the desired 
conc.rol strategy and thereafter the 
system performs all the monitoring and 
corrective action. There are limita- 
tions as pointed out in an article by 
Dr. A. G. Ratz (2) in which he states 
the following: 

If a system has N inputs, each of 
v;hich should be monitored and used in 
the control scheme, a problem can arise 
if only M of these inputs (N>M) can be 
analyzed each control cycle. For a 
typical gain control iteration, the CPU 
samples as many of the inputs as is 
practicable. This is done during the 
gathering time (t^ to (t^+Tf))i if 
there are too many inputs to be handled 
in an’ one gathering time, the CPU 
commutates through the inputs, handling 
a different set of inputs each control 
iteration until all inputs are covered. 
Thus, if these arc IJ inputs, only M 
of which can be handled during any one 
control iteration, it takes (N/M; con- 
trol iterations to update all inputs. 

Data for non-sampled inputs must be 
stored and used from previous control 
iterations, to compute averages, ex- 
ternal selections, etc. 

The above obviously leads one to under- 
stand that even digital control systems 
have limited corrective capability when 
the number of response channel measure- 
ments exceed the monitoring capability 
of the control system for each control 
cycle. 

2. Pretest Levels 

In Table 1, line number 9, notice 
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TABLE 1 Random Test Setup 





RANDOM EXECUTION READY 

1 

REi^D TEST TAPE ? 

2 

ENTER 

TEST IDENTITY (1/3 CHARACTERS) 

3 

ENTER 

TEST DATE (8 CHARACTERS) 

4 

ENTER 

OPERATOR IDENTITY (1^ CHARACTERS) 

5 

ENTER 

OTHER DESCRIPTORS IF ANY (1/3 CHARACTERS) 

6 

ENTER 

TRANSDUCER SENSITIVITY MV/G 

7 

ENTER 

NO. OF CONTROL CHANNELS 


ENTER 

NO. OF SCANS 


ENTER 

1 ST CONTROL CHANNEL 


ENTER 

MONITOR CHANNEL 


ENTER 

TYPE OF MULTICHANNEL CONTROL 


1. AVERAGE 2, PEAK 3. EXTREMAL 


IS EXTERNAL REFERENCE DESIRED ? ? 

8 

ENTER 

TOTAL TEST TIME IN MINUTES 

9 

PRETEST LEVELS AVAILABLE = 


1, 

-0 dB 


2 . 

-6 dB 


3. 

-12 dB 


4. 

-13 dB 


5. 

-28 dB 


6. 

-30 dB 


ENVER 1., 2., 3., 4., 5., OR 6. FOR SELECTION 

10 

CHANGE REFERENCE DESIRED ? ? 


BANDWIDTHS AVAILABLE = 


1* 

100 HZ 


2, 

200 HZ 


3. 

300 HZ 


4. 

500 HZ 


5. 

1000 HZ 


6, 

2000 HZ 


7. 

3000 HZ 


8. 

5000 HZ 


9. 

10000 HZ 


ENTER 1., 2., 3., . . . 9. FOR SELECTION 


RESOLUTIONS AVAILABLE 


1, 

W/200 


2. 

W/100 


3. 

W/50 


4. 

W/25 


ENTER 1., 2 ., 3., OR 4 . FOR SELECTION 


ENTER 

NUMBER OF SIECTRUM BREAKPOINTS 


ENTER 

FREQS. AND PSDS. ALTERNATELY 


ENTER NEXT FREQ. 


ENTER NEXT PSD 


PROGRAMMED RMS LEVEL = x,xxr 


CHANGE IN RMS DESIRED ? ’ 

11 

ENTER 

ALARM LEVEL IN DB 

12 

ENTER 

ABORT LEVEL IN DB 

13 

WRITE 

TEST TAPE ? 

14 

EXECUTE ? 
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that the control system allows the 
operator to specify to what level (in 
d&) less than full scale the operator 
would like to initiate a test. During 
this time the control system does a 
self checking operation to ensure that 
the entire system is within desired 
tolerance levels* After the system has 
accepted the level, it then increments 
itself in approximately +6 dB steps 
until it comes up to f u] 1 test level* 

In the event any abnormal responses 
occur, the system will warn the oper- 
ator of such conditions and shut down 
if necessary* 

3. Warning and Abort Valve Selection 
for Each Spectral L i ne of the 
Reference SpectrunP 

Generation of the reference spec- 
trum (Random Vibration Test) for a 
digital control system is similar to an 
analog system in that amplitude values 
are chosen as a function of frequency* 
However, the analog values are selected 
through a slide wire mechanism, whereas 
the values are typed in for a digital 
system. One major advantage wnich a 
digital system has over an analog system 
is that at each spectral line, AK\rm 
and Abort values can be specified and 
the system will constantly monitor 
these values and control to them. 

The alarm check makes the operator 
aware of test conditions which are out- 
side of the acceptable levels* The 
abort checks will cause a system to shut 
down. iSee figure 1), Each segment of 
the spectral profile is automatically 
monitored to the specified tolerances 
w’.ere the tolerances are programmed in 
dB. Alarm tolerances generally cause 
a data report to be printed if the con- 
trol amplitude exceeds the specified 
limit. The abort tolerance will cause 
a controlled shutdown and a message 
will be printed out. 

Alarm and Abort limits are gener- 
ally provided from approximately ±2 dB 
to ±®dB in adjustable steps. These 
limits are applied to the smooth PSD 
and not the "instantaneous" PSD, X©* 

To eliminate false alarms or false shut- 
downs, the alarm and abort signals are 
applied as follows: 

a. A running count n, is to be 
kept for each line (0<N), 
separate for alarm, and for 
abort • 


c. Each iteration of he control 
when the relevant limit level 
is not exceeded, n is decre- 
mented by one count. 

d. n can never go negative. 

e. When n reaches the value n =Nq# 
the relevant level is considered 
to be exce "lad, and the alaxm or 
abort action is initiated* 

f. A suitable value for Nq is three* 
No should be adjustable for spe- 
cial considerations, however, 
(Ref. (3) Application B’ lletin 
440 by Dr. A. G, Ratz; 1/23/74). 

From the above one can see that a fair 
amount of sophistication has been in- 
corporated into digital alam and abort 
specimen protection. 

4 . Gross Level Abort 

The gross level abort is an opera- 
tor selectable parcimeter that sets an 
overall RMS acceleration level which 
cannot te exceeded without terminating 
the test. Gross level aborts are used to 
protect both the system and test speci- 
men against runaway conditions. A 
loose accelerometer or inadvertently 
increasing a manual gain control causes 
the gross level abort to shut down the 
test in an orderly manner. 

5. Loop Integrity 

In the event an open loop occvirs in 
the response channel (s), catastrophic 
results could obviously come about. 

For single channel control, the system 
automatically detects the loss of con- 
trol signal at any level and initiates 
an abort. When multiple channels are 
used for control, generally each 
channel is checked for loss of signal 
at -3 dB of full scale and the system 
aborts if any loss of signal is detect- 
ed (Random Testing) . 

This check is made on control 
channel once every control loop cycle. 
These checks are made on input data 
prior to including the data in the con- 
trol calculations. Therefore, the out- 
put level will not change (increase) 
when a control signal is lost; only 
the abort action is allowed to change 
the output level to zero in a smooth 
manner. 


6. Power Fail Restart 

b. Each time (each iteration of 
the control) the alarm/abort 
level is exceeded, n is in- 
creased by one count. 
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Most control systems have an auto- 
matic power fail system. The system 
constantly monitors the line voltage- 












and when the line voltage goes out of a 
preset tolerance band, the system 
systematically and smoothly shuts down. 
The more sophisticated systems continue 
to monitor the voltage and bring the 
system back on line once the problem 
clears up. 

C. Housekeeping Operations - Monitor- 
ing Pertinent Equipment Parameters 
During a Test 

The items discussed consume, with 
drive signal updating, much of the con- 
trol system's time. However, there are 
still functions that a control system 
monitors. These functions or para- 
meters are near static or DC values 
and c^*nsequently can easily be monitor- 
ed by the system. A typical list 
might include the following: 

Amplifier door status 
Exciter overtravel status 
Exciter coolant flow status 
Exciter armature temperature 
Exciter field temperature 
Power amplifier coolant flow status 
Power amplifier heat sink temperature 
Field supply 
Amplifier overload 

1 . overcurrent 
2« overvoltage 
3. overpower 
Loss of power supplies 
Power amplifier low water warning 
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this device is to protec u the exciter 
and not the test article, in fact, 
wnen the protector is initiated , the 
deceleration levels generated are 
significantly high to cause possible 
damage to the test article. Goddard 
has consequently had two manufacturer's 
develop two Controlled Deceleration 
Specimen Protection Systems. 

The existing systems effect con- 
trolled deceleration by automatically 
sensing displacement, acceleration, or 
velocity for excessive levels, and when 
triggered, insert variable damping re- 
sistors across the armature circuit to 
effect the desired deceleration shut- 
down. 

CONCLUSION 

The safety protection of a test 
specimen has certainly been enhanced 
through the use of digital control 
systems. The primary proble... areas 
have been due to operator errors. To 
this end, the digital control system 
has been most successful, however, there 
are still areas of concern, namely, the 
"exciter dump" caused by transients 
out of the power amplifier or control 
system. This problem requires a system 
which can guarantee a controlled de- 
celeration shutdown to prevent damage 
to expensive test articles when trigger- 
ed by such transients. 


If any of these parameters go out of 
tolerance, a printout occurs stating 
the problo.a and shuts the system down 
if designated by che operator through 
software. 

II CONTROLLED DECELERATION SPECIMEN 
PROTECTION SYSTEMS 

As previously stated in the intro- 
duction, there are conditions in which 
the digital control system has limited 
ability to protect a test specimen. A 
condition of this nature occurs when a 
transient is generated within the power 
amplifier or output stages of the con- 
trol system and is directly transmitted 
to the exciter and the test article. 

To this author's knowledge there is no 
digital control scheme which can sense 
amplifier transients fast enough to 
initiate positive action to prevent the 
transient from being transmitted to the 
test article. The standard method used 
by amplifier/exciter manufacturer's to 
handle such transients has been to use 
the following technique. A concept 
known as ‘Armature Protection" is used 
in which the energy stored in the 
amplifier is discharged and the armature 
is short circuited. The purpose of 
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